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1.0 INTRODUCTORY COMMENTS 

In this part of the report, a first-principies based model is developed to predict the 
overpressure and temperature effects of a perforating orbital debris particle impact within a 
pressurized habitable module. While the effects of a perforating debris particles on crew and 
equipment can be severe, only a limited number of empirical studies focusing on space vehicles 
have been performed to date (see, e.g. [1-3]). Traditionally, crew loss or incapacitation due to a 
perforating impact has primarily been of interest to military organizations and as such have 
focused on military vehicles and systems (see, e g. [4]). 

The module wall considered in this study is initially assumed to be a standard Whipple- 
type dual-wall system in which the outer wall protects the module and its inhabitants by disrupting 
impacting panicles. The model is developed in a way such that it sequentially characterizes the 
phenomena comprising the impact event, including the initial impact, the creation and motion of a 
debris cloud within the dual-wall system, the impact of the debris cloud on the inner wall, the 
creation and motion of the debris cloud that enters the module interior, and the effects of the 
debns cloud within the module on module pressure and temperature levels. This is accomplished 
through the application of elementary shock physics and thermodynamic theory. 

Consider the normal hypervelocity impact of a projectile on the outer wall or ‘bumper’ of 
a dual-wall system as shown in Figure 1 . Upon impact, shock waves are set up in the projectile 
and bumper materials (Figure la). The pressures associated with these shocks typically exceed 
the strengths of the materials by several orders of magnitude. For example, in an 8 km/sec 
aluminum-on-aluminum impact, the ratio of the impact pressure (1 16.5 GPa=l. 15 MBar) to the 
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strength of the material (310 MPa for aluminum 6061-T6) is approximately 375, or roughly 2.5 
orders of magnitude. As the shock waves propagate, the projectile and bumper materials are 
heated adiabatically and non-isentropically. The release of the shock pressures occurs isen- 
tropically through the action of rarefaction waves that are generated as the shock waves interact 
with the free surfaces of the projectile and bumper (Figure lb). This process leaves the materials 
in high energy states and can cause either or both to fragment, melt or vaporize, depending on the 
material properties, geometric parameters, and the velocity of impact. 

In this manner, a debris cloud consisting of fragmented, melted, or vaporized projectile 
and bumper material (referred to from here on as the ‘primary debris cloud’ ) is created. This 
debris cloud travels towards, impacts, and perforates an inner bumper. For this study, this 
bumper consists of either a blanket of multi-layer thermal insulation (MLI) or a multi-material 
composite wall (‘Type A’ and ‘Type B’ inner bumpers, respectively). In addition to providing 
thermal insulation, the inner bumper also serves to trap some or most of the primary debris cloud 
particles. The debris cloud exiting the inner bumper is referred to from here on as the ‘secondary 
debris cloud’. 

Eventually, the secondary debris cloud impacts the inner or ‘pressure’ wall of the dual- 
wall system. The material in the secondary debris cloud and the portion of the pressure wall 
impacted by the secondary debris cloud are shocked and released in a manner similar to that of the 
initial impact on the bumper. In some instances this causes a perforation of the pressure wall. If 
this were to occur in an orbiting habitable module, then the debris cloud created as a result of the 
pressure wall perforation (hereafter referred to as the ‘internal debris cloud’) would be ejected 
into the module interior. In addition to the motion of the internal debris cloud within the module, 
a shock wave would also be created in the module due to the presence of the module air. 
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The focus of the model developed herein is on the motion of the internal debris cloud 
within the module and its effects on pressure and temperature within a pressurized habitable 
module. Thus, the model to be developed will consider: 

1) the shock loading and release of projectile and bumper material due to the initial impact 
on the bumper; 

2) the creation and motion of the primary debris cloud created by the initial impact on the 
bumper; 

3) the impact of the primary debris cloud on the inner bumper; 

4) the composition and motion of the secondary debris cloud; 

5) the shock loading and release of the secondary debris cloud and pressure wall material 
caused by the impact of the secondary debris cloud on the pressure wall; 

6) the module internal pressure and temperature increase at a specific location due to the 
passage of the shock wave induced by the motion of the module pressure wall; 

7) the module internal pressure and temperature increase at a specific location due to the 
passage of an internal debris cloud; and, 

8) the module internal pressure and temperature increase at a specific location due to the 
combined presence of an internal debris cloud and a shock wave. 

Section 2.0 presents the shock loading and release analysis that will be applied to the 
initial impact of the projectile on the bumper. The motion of the primary debris cloud is discussed 
in Section 3.0 with an emphasis on obtaining an appropriate velocity value to characterize the 
forward motion of the primary debris cloud. This is a critical value because it is eventually used 
as input for the pressure wall shock loading and release analysis, the results of which are in turn 
used directly in the modeling of the module internal effects due to a perforation (i.e. the pressure 
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and temperature increases). Thus, an incorrect primary debris cloud velocity value could have a 
disastrous effect on all subsequent calculations. Section 3.0 concludes with a discussion of the 
impact of the primary debris cloud on the inner bumper and the creation and motion of the 
secondary debris cloud. In Section 4.0, the shock loading and release analysis developed for the 
initial impact in Section 2.0 are applied to the impact of the secondary debris cloud on the 
pressure wall. Sections 5.0 and 6.0 present the modeling of the pressure and temperature 
increases due to shock wave motion and internal debris cloud passage, respectively. Section 7.0 
presents the manner in which the individual effects models developed in Sections 5.0 and 6.0 are 
combined to give the cumulative effect of shock wave motion and internal debris cloud passage. 
Finally, Section 8.0 presents a comparison of the predictions of the internal effects model and 
experimental data. Modifications to the model that are required to bring its predictions in closer 
agreement with the experimental results are also presented and discussed. 



2.0 INITIAL IMPACT SHOCK LOADING AND RELEASE MODELING 
2.1 Shock Loading Due to High Speed Impact 

In calculating the shock loading and subsequent release of the projectile and bumper 
materials, the shock waves are considered to be initially planar. This simplification allows one- 
dimensional relationships to be used for analyzing the creation and release of shock pressures. In 
this manner, the shock pressures, energies, etc., in the projectile and bumper materials are 
calculated using the three 1-D shock-jump conditions, a linear relationship between the shock 
wave velocity and particle velocity in each material, and continuity of pressure and velocity at the 
projectile/target interface. 

Specifically, if we consider the 1-D impact of a projectile with velocity v 0 on a sta ti onary 
target plate, conservation of mass, momentum, and energy across the shock fronts in the 
projectile and in the target yields 

Projectile Target 

Usp^op = ( u sp* u pp)^Hp u s/Vot = ( u st _u pt) /v Ht (la,b) 

P Hp = P op + u sp u p{Aop PHt = p ot + “stUpAot (2a,b) 

EHp = E op + (PHp+Pop)(Vo p -V Hp )/2 E Ht = E ot + (PHt + Pot)(V 0 t-V H t)/2 (3a,b) 

where V— 1/p is specific volume, u s and Up are shock and particle velocity, respectively; 

V 0 ,P 0 ,E 0 are the density, pressure and energy states associated with the shocked 
and initial material states. In equations (1-3), the subscripts 'p f , and *t' refer to projectile and target 
quantities, respectively. In the development of equations (1-3), the shock velocity in the projectile 
is taken relative to a 'stationary' projectile. 
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The linear shock velocity-particle velocity relationships for the projectile and bumper 
materials are taken to be in the form 

u s = c Q + ku p (4) 

where c 0 =V(KV 0 ) is the material bulk speed of sound, K=E/3(l-2v) is the adiabatic bulk 
modulus, E and v are Young's modulus and Poisson's ratio, respectively, and k is an empirically- 
derived constant. The general equations (1-4) are applied to the initial impact on the bumper of a 
dual-wall system in the following manner. Upon impact, pressure equilibrium at the projec- 
tile/target interface implies that 

p Hp = p Ht (5) 

while material continuity at the interface implies that 

v o = u pp + u pt (6) 

Because the individual bumper plates in a dual-wall system are free from any initial mechanical 
stress (they are merely supported at their four comers a fixed distance away from the inner 
pressure wall), the initial conditions ahead of the projectile and target shock waves are taken to be 
zero (with the exception, of course, of the initial material densities). Solving equations (1-6) 
simultaneously yields expressions for projectile and target particle velocities which can then be 
used to calculate shock velocities, pressures, internal energies, and material densities after the 
passage of a shock wave. For example, using this procedure to solve for initially for Upt yields 
Upt = (b-VA)/2a (7) 

where 

a = kp - k{(p 0 ^/p 0 p) (8a) 

b — 2kpV 0 + c 0 p + c ot(Pot^Pop) (8b) 
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(8c) 


A - b 2 - 4a(c 0p v 0 +k p v 0 2 ) 
Then it follows that 


u pp = v o ' u pt 

(9a) 

U St = C Ot + kt u pt 

(9b) 

u sp = c op + kp u pp 

(9c) 


The shocked densities of the projectile and target materials are found by substituting equations 
(7,9a-c) into equations (la,b) to yield 


PHp = 1 / V Hp 


u sp ! V 0 p 
u sp ' u pp 


(10a) 


PHt =1/v Ht = 


u st j Vpt 


u st -u pt 

Finally, equations (2a,b) and (3a,b) are then used to define the pressure and energy in the 
projectile and target materials, respectively, associated with the passage of the shock waves 
created by the initial impact. This completely defines the shocked states of the projectile and 
target materials due to the initial impact. 


(10b) 


While the shock loading of a material is an irreversible process that results in an increase 
of the internal energy of the shocked material, the release of a shocked material occurs 
isentropically along an 'isentrope' or 'release adiabat'. The difference between the area under the 
isentrope and the energy of the shocked state is the amount of residual energy that remains in the 
material and can cause the material to melt or even vaporize. In order to calculate the release of 
the projectile and target materials from their respective shocked states (each characterized by Pjj, 
^H> an< ^ Vh), an appropriate equation-of-state is needed for each material. To keep the analysis 
relatively simple, the Mie-Gruneisen equation-of-state [5] was used in this study. 
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2.2 Shock Release Using the Mie-Gruneisen Equation-of-State 

The Mie-Gruneisen equation-of-state (EOS) is an accurate thermodynamic description of 
most metais in the solid regime and is relatively easy to use. It has the form 

P = p H + P r (H-Ep!) (11) 

where the time-dependent Gruneisen coefficient T is given for most metals as 

r = r oPo /p (12) 

where r o =Kp/p 0 Cp is the ambient Gruneisen coefficient, K is the adiabatic bulk modulus, 

P=3a is the volumetric coefficient of thermal expansion, and Cp is specific heat at constant 
pressure. Invoking the Second Law of Thermodynamics 

dE = TdS - PdV (13) 

along with the isentropic constraint dS=0 for the release process allows us to construct the release 
isentrope in P-V space for a material referenced to the material Hugoniot in P-V space and a 
given initial shocked state defined by Pjj, Vjj, Ejj. Using the procedure outlined in Reference 
[5], the pressure P[ and internal energy Ej at a specific position *i* along the isentrope can be 
shown to be given by 

Pi = [PHi + (r/V)i(Ei_i - Pj.jAV/2 - EHi)]/[l-Kr/V)iAV/2] (14) 

where AV is tiie incremental change in volume used to create the release isentrope, and Pjjj and 
Effi are the pressure and energy along the Hugoniot corresponding to the i-th position in the 
release process. The release process is continued using equation (15) until the release isentrope 
so determined crosses the V-axis. 

It should be noted that based on its thermodynamic origins, the Mie-Gruneisen EOS 
cannot be expected to give accurate results in a highly expanded liquid regime or in a vapor re- 
gime. This is because as impact energy increases, the assumption that the Gruneisen coefficient is 
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a function of density alone is no longer valid. At high impact energies, the Gruneisen coefficient 
is a function of internal energy as well as density. Experience has shown, however, that it does 

yield fairly accurate end-state results even when there is a small percentage of molten material 
present [6], 
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Figure 1. Hypervelocity Impact of a Module Wall System 
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3.0 MODELING OF PRIMARY AND SECONDARY DEBRIS CLOUDS 

3.1 Primary Debris Cloud Mass 

The mass of the primary debris cioud consists of the mass of the original impacting 
projectile plus the mass of the removed bumper material. No mass is considered lost to 
backsplash of the bumper and projectile materials; hence, since all of the mass is presumed to be 
directed in towards the pressure wall, the model to be developed should yield conservative results. 
The contribution of the bumper material to the primary debris cloud mass can be determined once 
the bumper hole diameter is known. This diameter can be calculated using any one of a number of 
empirical equations for hole diameter in a thin plate due to a high speed impact (see, e g. [7]). 

The particular equation used in this study is given as follows [8]; 

D h /d p = 3.4(t b /dp)°-333( Vp / Cb) 0.333( 1 o-0.0308 Pb /p p ) (15) 

where cjj is the speed of sound in the bumper material. This equation was chosen because it is 
applicable over a wide range of projectile/bumper material combinations, impact velocities, and 
tt/dp ratios. 

3.2 Primary Debris Cloud Characteristic Velocities 

The characteristic velocities of interest for the debris clouds considered in this model are 
the axial and expansion velocities. These values can be determined in any number of different 
ways. Naturally, the various methods provide slightly different results. The ‘correct’ or most 
appropriate method is that which will provide characteristic velocity values that when used in 
subsequent calculations will provide pressure and temperature increases that are most consistent 
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with experimental results. The following sections present three methods to calculate 
characteristic velocities specifically for the primary debris cloud. 

3.2.1 Conservation of Momentum and Energy 

Because the initial impact occurs at such a high velocity, momentum transfer to the 
bumper itself can be ignored. As a result, the momentum balance along the axial direction before 
and after the initial impact is written as follows: 

m p v o = (m p +m b ) vjJJJ. 111 (16) 

where vj^. m is the momentum conservation based axial velocity for the primary debris cloud. In 


fact, it is also the debris cloud center-of-mass velocity. Solving equation (16) for v£[J. m yields 


v mc m " m p v o / ( m p +m b) 


(17) 


We note that v{J[J. m can be viewed as the average axial velocity of the entire debris cloud mass. 

Furthermore, because of the rapidity of the impact event, the only significant energy losses 
are to the accompanying light flash and the shock heating of the bumper and projectile mate rials 
If we neglect the energy associated with the light flash, then the following energy balance can be 
written before and after the initial impact. 

p. _ psh&rl , psh&rl mi c prim, ax , _ prim, exp /lox 

impact - E proj m p + E bmpr m b + Ej^ + E d ebc l d ( 18 > 

where Ej mpact =m p v 0 2 /2 is the energy of the impacting projectile, m p is the projectile mass, m b is 

the bumper hole out mass, E p ^ ^ and E^^j? are the residual energies per unit mass in the 

projectile and bumper materials following shock heating and release; and e£"™|® xp 


'debcld 
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are the energies associated with the axial motion and radial expansion, respectively, of the primaiy 

debris cloud. By neglecting energy losses such as those due to light flash, the results obtained 
herein should also be conservative in nature. 

If we let vP? p m denote the expansion velocity of the primary debris cloud, then we can 

write 


v 


prim 

exp 


E debdd" P =('”p + "'bXv^ p m ) 2 /2 

E d2dd X -(-p + »bX»fi") 2 /2 

Substituting equations (19a,b) into equation (18) and solving for v cxp yields 




[E p^ rim P + E bmpr m b + < m p + m b)(vEr c m > 2 l2 U ' K m p + m b )/2] 


= T K 


^o 2 /2- 


(19a) 

(19b) 


( 20 ) 


We note that the expression under the radical in equation (20) is the difference of two 
terms. Because this difference is in an expression that is under a radical, there are three 
possibilities of which the user of the model being developed herein must be aware. 

First, the difference is positive. In this case, the radical yields a real value that corresponds 
to the expansion velocity of the primary debris cloud. 

Second, the difference is zero. In this case, the kinetic energy of the impacting projectile 
IS not high enough so that there is no energy available for debris cloud expansion. This would 
result if the impact were at a velocity low enough so that the projectile passed through the bumper 
relatively unscathed (i.e. without any significant lateral dispersion). 

Third, the difference is negative. In this case, the radical yields an imaginary value; 
however, the expansion velocity must be a real number. This implies that there is an error in the 
portion of the model that has been developed thus far. For example, depending on the impact 
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velocity, it is conceivable that only a portion of the projectile mass and only some fraction of the 
bumper hole-out mass are shocked and released. Implementing some adjustment factors to 
account for partial shock heating and release for m p and m^ could significantly reduce the amount 
of energy being subtracted in equation (20). Or, for example, the hole diameter equation may be 
producing a diameter value that is too large. Correcting the hole diameter equation would also 
alter the amount of energy being subtracted. Finally, it is possible that the debris cloud axial 
velocity is too large. Including the momentum due to backsplash and the momentum transferred 
to the bumper plate would certainly alter the value of the debris cloud axial velocity. 

Finally, we note that the spread of the primary debris cloud (in this case the half-angle 
defining the spread) can be estimated using the following, relationship between the calculated 
primary debris cloud expansion and axial velocities: 


eP rim =.an-l(vPj D m /vS,r) 


3.2.2 Elementary Shock Physics 

As the shock wave created by the initial impact strikes the rear surface of the bumper, it 
creates a rarefaction wave that travels back into the bumper and eventually in some form back 
into the projectile. This action and interaction of the shock wave and the free surface impacts a 
velocity uf s ^ to the target rear surface equal to the sum of the particle velocity in the target 
material due to the shock wave u p t and the particle velocity due to the rarefaction wave u^: 

u fst = u pt + u rt (22) 

Since u rt *u pt [5], an alternative form for equation (18) is 

ufst = 2u pt (23) 

As in Reference [9], because the material along the leading edge of the primary debris 
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cloud is in fact the material at and directly behind the bumper rear free surface, the velocity of the 
leading edge of the primary debris cloud vP nm is approximated with u fst , that is. 


v le =u fst = 2u pt (24) 

Unfortunately, while the leading edge velocity may be used to characterize the axial velocity of 
the primary debris cloud, there does not exist a means by which the debris cloud expansion 
velocity can be calculated using the equations of elementary shock physics. 

3.2.3 Effective Applied Pressure Wall Velocity 

In a manner similar to Reference [10], the momentum of the primary debris cloud material 
before impact with the pressure wall is balanced with the momentum of the debris cloud material 
and the portion of the pressure wall affected by the debris cloud impact after the debris cloud 
impact. In this calculation, the presence and effect of an inner bumper are temporarily ignored. 

In this manner, the following expression is obtained: 


(m p +m b ) v P n c m - (m p +m b +m w ’)v pw (25) 

where m w ’ is the pressure wall mass affected by the debris cloud impact and v pw is the resultant 

velocity applied to the wall. Substituting for vj}[J. m according to equation (17) and solving for 
Vpw yields 

v pw = ropVo/Crop+mb+mw’) (26) 

The portion of pressure wall material affected by the debris cloud impact is taken to be the mass 
of the wail material directly beneath the impact footprint of the debris cloud on the wall, that is, 
m w = *( r w’) 2 t w Pw 

where 


(27a) 
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r w’ = Stan0P" m (27b) 

In equation (27b), S is the stand-off distance between the outer bumper and pressure wall. 

3.2.4 Observations and Comments 

A parametric study was performed using a 1.0 cm diameter aluminum projectile impacting 
an ail-aluminum dual-wall system with a 0. 16 cm thick bumper, a 0.48 cm thick pressure wall, and 
a stand-off distance of 1 1.43 cm at velocities between 1 and 16 km/sec. Table 1 provides the 

values of v , v j^p 1 , v|^ nm Vp W for the velocity regime considered. 

Table 1. Comparison of Primary Debris Cloud Characteristic Velocities 


v o 

(km/s) 

prim 

v mc 

(km/s) 

prim 

v exp 

(km/s) 

(km/s) 

v pw 

(km/s) 

1 

0.896 

0.276 

1.0 

0.111 

2 

1.689 

0.644 

2.0 

0.160 

3 

2.428 

1.048 

3.0 

0.198 

4 

3.096 

1.474 

4.0 

0.230 

5 

3.735 

1.916 

5.0 

0.258 

6 

4.340 

2.372 

6.0 

0.282 

7 

4.917 

2.840 

7.0 

0.305 

8 

5.467 

3.317 

8.0 

0.326 

9 

5.996 

3.802 

9.0 

0.345 

10 

6.504 

4.296 

10.0 

0.363 

11 

6.994 

4.797 

11.0 

0.379 

12 

7.468 

5.303 

12.0 

0395 

13 

7.926 

5.815 

13.0 

0.410 

14 

8.370 

6.330 

14.0 

0.424 

15 

8.801 

6.849 

15.0 

0.438 

16 

j 9.220 

7.370 

16.0 

0.451 


As can be seen in Table 1, Vj^ nm =v 0 for like-on-like impacts so its use as the debris 
cloud impact velocity would most likely yield unreasonably high pressure wall impact stresses. 
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Conversely, v pw is very low and would yield pressure wall impact stresses that would probably be 
too low to transfer significant amounts of energy into a module following a pressure wall 
perforation. The characteristic velocities produced by the momentum and energy conservation 

based calculation appear to be more physically reasonable and will be used to characterize the 
impact velocity of the primary debris cloud. 

3.3 Secondary Debris Cloud Mass 

The mass of the secondary debris cloud consists of a portion of the primary debris cloud 
mass plus the mass of the removed inner bumper material. No mass is considered lost to 
backsplash of the tnner bumper and primary debris cloud material. As discussed in Section 1.0, 
the inner bumper serves to trap a certain portion of the primaty debris cloud particles. This 
fraction of trapped particles can be significant, especially for Type B inner bumpers. The 
reduction in mass of the primary debris cloud as it transforms into the secondly debris cloud is 

obtained using a simple multiplicative reduction factor. The nantre of this factor, denoted by 
is discussed in Section 3.3.1. 

The contribution of the inner bumper material to the secondary debris cloud mass can be 

determined once the inner bumper hole diameter is known. This diameter is calculated using the 

following fundamental empirical equation for hole diameter in a thin plate due to a high speed 
impact [11]; 

df/dp = 0.45V(t s /dp)° 666 + o.9o 

where V is the impact velocity and , s is the thtckness of the impacted plate. This equation was 

chosen because it is applicable over a wide range of projectile/bumper material combinations, 
impact velocities, and t s /d p ratios. 
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To render equation (28) applicable to the impact of the primary debris cloud on the inner 
bumper, some simplifying assumptions must be made. First, the velocity V in equation (28) is 
taken to be the characteristic velocity of the primary debris cloud as given by equation (17). 
Second, because the inner bumper is placed relatively close to the outer bumper, the primary 
debris cloud is in a relatively compact form when it impacts the inner bumper. This can give rise 
to a fair amount of “thickness” in the leading edge of the primary debris cloud as short-time X-ray 
photographs of expanding debris clouds have previously shown (see, e.g., [12]). Thus, for the 
purposes of this model, it is assumed that the thickness ratio is equation (28) is unity, or 
t s / d p =1.0. Finally, it is assumed that the perforation of the inner bumper occurs primarily as a 
result of the impact of the primary debris cloud leading edge material that is concentrated near the 
centerline. For the purposes of this study, this material is presumed to be contained within an area 
whose diameter is approximately equal to that of the hole in the bumper plate. In this manner, the 
projectile diameter” in equation (28) is taken to be equal to Djj as calculated by equation (15), or 
d p =D h 

Making use of the terminology and simplifying assumptions discussed in the preceding 
paragraph reduces equation (28) to the following form: 

dh = (0 45vP n c m +0.90)D h (29) 

Thus, the total mass of the secondary debris cloud is given by 

m dScfd =K l( m P +m b) + *Aib d h 74 ( 3 °) 

where is the areal mass density of the inner bumper material. The component of the 
secondary debris cloud whose origin can be traced back to the primary debris cloud is hereafter 
referred to as the first component while the inner bumper material contribution to the secondary 
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debris cloud is hereafter referred to as “the second component”. 

3.3.1 The Mass Factor kj 

From the discussion of the experimental resuits in Volume I of this report, it is apparent 

a Type A inner wai, traps only a small amount of the pnmary debns cloud material, while 

Type B traps a more significant quantity of matenai. One of the most distinguishing differences 

between the two inner bumper types (and a direct consequence of their respective composition) is 

.heir area, densities: A A =0.072 gm/cm 2 for the Type A inner bumper while A B =,.02 gm/cm 2 for 

■he Type B inner bumper. Thus, i, appears that the trapprng capability of the inner bumpem is 

directly related to their areal densities. The following expression for x, appears to capture this 
feature: 


kj= 1- 


^ib 


( 31 ) 


where n is a positive rational number, A ib is the areal mass density of the inner bumper, and 
Aw=Pw«w is .he areal mass density of the pressure wall. For example, if we consider an 
aluminum pressure wall with < w -0.48 cm and n-1, then if A ib =A A , i, follows that x,=0.946, 
which IS only a slight reduction in debns cloud mass (approximately 5%). Alternatively, if 

A, b Afi ' the " lq_0 ' 237 ’ w h'eh is a more significant reduction in debris cloud mass 

(approximately 76%). Ultimately, the validity of this manipulation and, more specifically, the 

appropriate choice of the exponent will be ascertained when the predictions of the model are 
compared with experimental results. 

3.4 Secondary Debris Cloud Characteristic Velocities 

As in the case of the primary debris cloud, the characteristic velocities of interest for the 
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secondary debris cloud are the axial and expansion velocities. These values can be determined in 
a manner analogous to the procedure presented in Section 3.2.1. However, there is a significant 
difference between the method developed herein for the secondary debris cloud and that for the 
primary debris cloud in Section 3.2.1. This difference stems from the difference in the material 
composition of the two debris clouds. 

For the purposes of this model, the primary debris cloud was tacitly assumed to be 
homogeneous and spherical; although it could possibly contain two different materials, both the 
projectile and bumper material that comprised the primary debris cloud were presumed to be 
evenly distributed throughout the cloud. As a result, only one axial velocity and one radial 
velocity were required to characterize its motion and expansion, respectively. However, the 
secondary debris cloud is comprised of two very distinct material types: the primary debris cloud 
material that passes through the inner bumper (i.e. the first component) and the inner bumper 
material itself (i.e. the second component), which is typically some form of composite material. 
Thus, in this case, it is necessary to develop axial and radial velocities for each component of the 
secondary debris cloud. 

The procedures for calculating the axial and expansion velocities are very similar. In both 
cases, we first determine an average velocity quantity of interest (i.e. the average axial motion 
velocity or the average expansion velocity for the entire secondary debris cloud). In the case of 
the axial velocity, this is done by conserving momentum before and after the impact of the primary 
debris cloud on the inner bumper; for the expansion velocity, this is done using energy 
conservation before and after the inner bumper impact. Next, we postulate in both cases that the 
velocity quantities of interest for both components are proportional to the velocity quantity just 
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calculated and solve for the proportionality constants by conserving appropriate momentum and 
energy quantities. 


As m the case of the primary debris cloud calculations, losses due to the light flash and 
backsplash accompanying the impact of the pnmary debris cloud on the inner bumper are ignored, 
as are losses due to the shock heating (or burning) of the primary debris cloud and inner bumper 
materials. As a result, the state of the primary debris cloud material that does pass through the 
inner bumper to become par, of the secondaiy debris cloud remains unchanged by the impact on 
the second bumper. However, we do include the lost energy of the trapped primary debris cloud 
material is included in the energy balance equations. The next two sections present the details in 

this procedure for calculating the secondaty debris cloud axial and expansion velocities. 

3.4.1 Axial Velocities 

To begin, we apply momentum conservation in the axial direction before and after the 
impact of the primary debris cloud on the inner bumper: 

(mn+mut vP nm - m sec on sec 

t m p m b/ v mc - m debc|d v mc (J2) 

where m deb eld is 8>«" Ration (30) and where v is the average axial velocity for the 
secondary debris cloud. Solving equation (32) for v|^§ yields 


v mc = ( m p +m b ) 



sec on 
debcld 


(33) 


We note that momentum transfer to the inner bumper itself has again been ignored due to the 

rapidity of the impact event. We now rewrite equation (32) in terms of the axial velocities of the 
individual secondary debris cloud components: 
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(34) 


(m p + m b )vjj 1 n c m = Ki(m p +mb)v ^ 1 + mj b v ^£ 2 

where v mc and v^ c refer to the axial velocities of the secondary debris cloud components, 
respectively, and where from equation (30) 

™ib =*Aibdh /4 (35) 

Next, we write 

seel — /■> sec 

v mc M v mc (36a) 

and 


v 


sec 2 
me 


-C2 v 


sec 

me 


(36b) 


where Cj and C 2 are the proportionality constants that need to be determined. Substituting 


equations (36a,b) into equation (34) and rearranging yields 


Kj(m p + m b )C 1 +mj b C 2 =(m p 


+ m b) 


prim'i 
r mc 1 


1 V sec 1 

\ v mc J 


(37) 


A second equation for C \ and C 2 is obtained by conserving the kinetic energy associated with the 
axial motion of the primary and secondary debris clouds. This equation is written as follows: 

(m p -rm b )(vST c m ) 2 /2 = Ki(m p +m b Xvmc 1 ) 2 +mj b (v ^ 2 ) 2 +Ej ) ^ ax (38) 

where the final term on the right hand side of equation (38) refers to the axial kinetic energy lost 
due to the fact that the inner bumper stops a portion of the primary debris cloud mass. This 
energy loss is written as 

E&r = 0 - K 1 X m p + m b )(v!T ) 2 / 2 (39) 

Substituting equations (36a,b) and (39) into equation (38) and rearranging yields 
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2 


K 1 ( m p + m b )C ^ + m^C ^ = k i (m p + m^)j 


pnm i 
me 

i v sec 
V v mc J 


(40) 


Solving equations (37) and (40) yields 


c, = k - gC 2 

a 


(41a) 


c 2 = 


£a 

-±Va' 


a 

3 

GH 


(41b) 


where 


A = 




(42) 


and where 


a K^mp+mjj), P ~ m ib> 


£a~( m p + mb) 


pnm 

v mc 

v sec ; 

v mc / 


(43a,b) 


Thus, equations (36a,b), together with equations (41a,b), (42), and (43a,b,c) define the axial 
velocities of the two secondary debris cloud components. 

3.4.2 Comments on the Axial Velocity Calculations 

In solving equations (37) and (40) for the constants C { and C 2 , the manipulations result in 
a single quadratic equation for C 2 . The general solution of such an equation has two roots as 
indicated by the ± in equation (41b). Naturally, the choice of sign in equation (41b) will affect the 
outcome of Cj according to equation (41a). Since Cj must be a positive number, it follows from 
equation (41b) that C 2 < / 3 . Substituting for C 2 according to equation (41b) yields the 
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inequality Wa < 1 . Obviously, if we keep the sign, this inequality is always true. Of course, if 
we keep the *+’ sign and substitute for A according to equation (42), we will obtain a relationship 
between the mass reduction constant kj, the mass of the primary debris cloud, and the mass of the 
first component of the secondary debris cloud that will also make the inequality true. For the sake 
of simplicity, the model will consider only the sign in equation (41b) since this will also allow 
the model to operate in its most general form. 

Once it has been decided to keep the sign in equation (41b), another issue arises, 
namely, we must ensure that C 2 is, in fact, a positive number. According to equation (41b), C 2 
will be positive so long as 


£-Va>o 

a 



Substituting for A according to equation (42) and simplifying, the inequality in equation (44) 
reduces to simply 1 > iq2 Since the definition of >q is based on the premise that iq < 1, the 
inequality in equation (44) is always true, which in turn implies that C 2 will always be a positive 
number if we keep the sign in equation (41b). 

The final issue that needs to be discussed is raised by the presence of the radical in 
equation (41b), namely, we must be sure that A > 0. According to equation (42), this will be true 
only if 



1 



a 


( 45 ) 


Substituting for a and p in equation (45) according to equations (43a,b) yields the following 
quadratic inequality: 
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where 


(46a) 


K1 2 


p + yKi - 1 > o 


y = m ib 


m p + m b 

Therefore, in order for the inequality in equation (45) to be true, we must have 


(46b) 


M > 


-y + y/r 2 + 4 


2 (47) 

We note that, since x, > 0 by definition, the root with the V sign in front of the radical has been 
discarded. In addition, the discriminant in equation (47) is always positive (so there is no need to 
be concerned about complex roots) and the numerator is always positive (so there is no need to 
be concerned about negative iq values. Thus, if the choice of the exponent V in equation (31) 
results in a mass reduction factor k, that violates equation (47) (i.e. the exponent is too large and 

the resulting reduction factor is too small), then the exponent needs to be reduced in order to 
increase the reduction factor value to an acceptable level. 

3.4.3 Expansion Velocities 

To begin, we conserve energy before and after the impact of the primary debris cloud on 
the inner bumper. This energy balance is written as follows: 


E. pnm = g sec, ax , p sec, exp prim 
impact ^debcld + fc debcld + E lost 


(48) 


where E j mpact is the total kinetic energy of the impacting primary debris cloud (i.e. axial motion 

and expansion components combined), while E**“ and Eg** refer to the kinetic energies 
associated with the axial motion and radial expansion, respectively, of the secondary debris cloud. 
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The lost energy of the trapped primary debris cloud particles, E , consists of both 
components of primary debris cloud motion, and is given by 

e 1 P oT = 0 - K 1 )(m p + mb XviT ) 2 / 2 + (1 - K 1 X m p + m b XvS?™ ) 2 / 2 (49) 


If we let v eX p denote the average expansion velocity of the secondary debris cloud, then 


we can write 


P sec, exp 
^debcld 


_ sec on 
m debcld 


fv sec 

t v exp 


) 2 / 2 


(50a) 


where m debcld is given by equation (30). Since 


P sec, ax 
aebcld 


(50b) 


we substitute equations (49) and (50a,b) into equation (48) and solve for v to obtain 

v Sp=^{M[(mp +m bX v mc m ) 2 +(m p +m b Xvgj n ) 2 ]-m^ l d (v^) 2 }/ mjg^ (51) 

We again note that the expression under the radical in equation (51) is the difference of two 
terms. Because this difference is in an expression that is under a radical, care must again 
be taken to ensure that the difference is in fact positive. 

To determine the expansion velocities of the two secondary debris cloud components, we 
again write 

v ixp 1= Dj v| X p (52a) 

and 

v exp 2= E > 2 v Ixp (52b) 
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where 1 and v|®£ 2 are the expansion velocities of the two secondary debris cloud 

components, and Dj and D 2 are the corresponding respective proportionality constants that need 
to be determined. 


Next, we conserve the kinetic energy associated with the expansion of the two secondary 
debris cloud components: 

(m p +m b )(v^) 2 /2 = jq(m p +m b )(v|^ 1 ) 2 /2 + mj b (v|^ 2 ) 2 /2 + E 1 pr ^ exp (53) 

where the final term on the right hand side of equation (53) refers to the expansion kinetic energy 

lost due to the fact that the inner bumper traps a portion of the primary debris cloud mass. This 
energy loss is written as 

E & eXP =( 1 -MXmp + m b )(v e P x 7) 2 / 2 (54) 

Substituting equations (52a,b) and equation (54) into equation (53), yields 

K l( m p +m b )D 2 +mj b D 2 = k j(m p +m b ) 


prim 

'exp 


sec 


(55) 


A second equation for Dj and D 2 is obtained by conserving the momentum associate with 
the expansion of the primary and secondary debris clouds: 

(m p +m b )v p 7 =K!(m p +m b )v|^l +mibv |©c2 (56) 

Substituting equations (51a,b) into equation (56) yields 


K l( m p +m b )Dj +mj b D 2 =(m p +m b ) 


( prim^ 
v exp 


sec 

V v exp j 


(57) 


Thus, the constants Dj and D 2 are found by solving equations (55) and (57). However, if we 


27 



compare the equation pair (55,57) to the equation pair (37,40), we find them to be identical with 
the exception that equations (55) and (57) are written in terms of expansion velocities while 
equations (37) and (40) are in terms of axial motion velocities. Other than that difference, their 
basic forms are identical. Therefore, the solution for the constants Dj and D 2 is written down 
directly from equations (41) through (43) as follows: 

D, ^e-PP2 


D 2 = 


a 


a 


(58a) 


0l| +1 


(58b) 


where 


A = 




K 


2 

1 



and where 


a = Ki(m p +mb), (3 = m; b , 


£e 


f 

(m p + m b ) 

V. 


prim 

'exp 

sec 

v exp j 


(59) 


(60a,b) 


Thus, equations (51a,b), together with equations (58a,b), (59), and (60a,b,c) define the axial 
velocities of the two secondary debris cloud components 

3.5 Debris Cloud Spread Angles 

Finally, the spread angles of the two secondary debris cloud components (in this case the 
half-angles defining the spreads) can be estimated using the following relationships between the 
calculated secondary debris cloud expansion and axial velocities: 
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(61a) 


9 dc Cl ~ tan ' 1 ( v exp l /^mc 1 ) 

9 dc C2 * tan ' l ( v Sp 2 /vS 2 ) (62b) 

Based on equations (27a,b), the portions of pressure wall material affected by the impacts of the 
two secondary debris cloud components are given by: 


m w 1 *i\v 1 *wPw 

m w2 ~ 7Cr w2^ t wPw 
respectively, where 

r wl = ^2 tan 9^ C * 

r w 2 = S 2 tan 6 ^ c ^ 


(63a) 

(63b) 

(64a) 

(64b) 


In equations (64a,b), S 2 is the stand-off distance between the inner bumper and pressure wall. 
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4.0 PRESSURE WALL IMPACT SHOCK LOADING AND RELEASE MODELING 


In modeling the shock loading and subsequent release of the secondary debris cloud and 
pressure wall materials, it is again presumed that the shock waves are initially planar. This 
simplification again allows one-dimensional relationships to be used for analyzing the creation and 
release of shock pressures. The process used to calculate shock pressures, energies, etc. in the 
debris cloud and pressure wall materials is identical to that used in modeling the initial impact of 
the projectile on the bumper, with the following exceptions. 

1) Only the secondary debris cloud material component consisting of primary debris cloud 
material that passed through the inner bumper contributes to the impact loading of the pressure 
wall. As such, the properties of that component of the secondary debris cloud and pressure wall 
properties are used in place of projectile and bumper properties, respectively. 

2) The initial density of the secondary debris cloud material considered in the impact 
analysis is taken as the average of the final densities after release from their respective shocked 
states of the projectile and bumper materials in the primary debris cloud. 

3) The secondary debris cloud axial velocity calculated using the momentum conservation 
method developed in Section 3.4.1 is used as the impact velocity in the impact analysis. 

4) The component of the secondary debris cloud consisting of inner bumper material 
contributes to the loading of the pressure wall in the form of an impulsively delivered velocity. 
This causes additional motion of the pressure wall and so affects only the shock wave created at 
the pressure wall/module wall interface. 
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5) Unlike the stress-free bumper plate, the pressure wall of a pressurized habitable module 
experiences an initial bi-axiai state of stress due to the internal ambient pressure. However, 
because the pressure wall is thin (i.e. on the order of 5 mm in thickness), it may be assumed that 
the stress through the thickness of the plate is nearly zero and may be ignored. Since the shock 
waves in the pressure wail propagate through the thickness as well, it appears reasonable to 
presume that the initial conditions ahead of the pressure wall shock wave are also zero. 

6) The release of the pressure wall material from its shocked state terminates when the 
pressure m the wall material reaches a level corresponding to that which induced in the wall by the 
motion of the rarefaction wave through the shocked material. Since the pressure wall is backed 
by module air at 1 atm, this value is not zero as in the case of the bumper which was surrounded 
by a vacuum. The means of obtaining this value is presented in the next section. 
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5.0 MODELING OF INDIVIDUAL INTERNAL EFFECTS: SHOCK WAVE PASSAGE 
As discussed in Chapter 4, the magnitude of the shock wave induced in the module by the 
motion of the pressure wall is the combined result of two phenomena. The first is the impulsive 
loading imparted to the pressure wall by the impact of the secondary debris cloud inner bumper 
material component. The second is the impact of the secondary debris cloud material component 
consisting of primary debris cloud material on the pressure wall. The motion of the pressure wall 
due to the impulsive loading is discussed in Section 5. 1 while the effects of the impact loading are 
discussed in Section 5.2. The application of these two phenomena to the creation, propagation, 
and attenuation of a shock wave in the module interior is presented in Sections 5.3 and 5.4. 

5.1 Pressure Wall Impulsive Loading Effect 

In a manner similar to Section 3.2.3, the effect of the impulsive loading imparted to the 
pressure wall by the impact of the second secondary debris cloud component is characterized by 
an instantaneously applied velocity. This velocity is obtained by balancing the momentum of the 
second component of the secondary debris cloud material before its impact on the pressure wall 
with the momentum of the debris cloud material and the portion of the pressure wall affected by 
the impact after the debris cloud impact. In this manner, the following expression is obtained: 

m ib v me 2 = ( m ib +m w2)up W (65) 

where m ib> an d m w2 are given by equations (36b), (35), and (63 b), respectively, and Up W 

is the resultant velocity that is applied to the wall. Solving for u pw yields 

u pw = m ib v me ^ /( rn ib +m w 2 ) (66) 
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The application of this velocity to the creation of the module interior shock wave is presented in 
Section 5.3. 

5.2 Pressure Wall Impact Loading Effect 

As the shock wave generated by the impact of the first component of the secondary debris 
cloud on the pressure wall front surface propagates through the pressure wall, it encounters the 
pressure wall rear surface. In addition to creating a rarefaction wave that is reflected back into 
the pressure wall to relieve shock pressures, a shock wave is transmitted into the air that is 
adjacent the pressure wall rear surface. The properties of the reflected and transmitted waves are 
found using a technique based on the method of impedance matching (see, e.g., [9]). In this 
technique, continuity of pressure and particle velocity are enforced at the pressure wall/module air 
interface. 

The equations governing the reflection and transmission of shock waves at material 
interfaces are presented in the following section. In this development, all materials are assumed to 
be stationary. As such, the superposition of the effect of an additional pressure wall velocity is 
deferred until Section 5.3 where the results of the development to follow are combined with the 
results of the preceding section. 

5.2.1 Impedance Matching Technique 

Figure 2a shows an incoming shock wave in material A, a reflected wave from the A/B 
interface, and a transmitted shock wave into material B. Shock wave I is shown moving into 
undisturbed material (denoted with a ’O' subscript). The shocked state of the material in its wake 
is denoted by a T subscript. The reflected wave II moves back into this shocked material and 
leaves behind it material whose state is denoted by a subscript of '2'. The transmitted shock wave 
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Ill moves into undistributed material (denoted by a '4' subscript); the condition of the shocked 
materia] behind it is labeled with a '3' subscript. Figure 2 b shows the same configuration only all 
motion is shown under steady conditions. 


Across shock front I we have 

Pl(Ui-u 1 ) = p 0 A (U 1 -uo) (67a) 

p 1 - p 0 = Po^KU i -uo)(u i -u 0 ) (67b) 

as well as the constitutive relationship between the shock wave speed Uj and the particle velocity 
u} induced in the shocked material 

u l=c 0 A + k A ui ( 68 ) 

In equations (67a,b), po has been replaced by p 0 A , the ambient density of material A. Assuming 
stationary conditions at zero pressure ahead of the shock wave (i.e. uo=Po= 0 ) and that Pj is 
known (it is the pressure associated with the impact of the initial debris cloud on the pressure wall 
as calculated in Section 4.0), equations (67-68) can be used to solve for u\ and p j as follows: 


U 1 = -(c 0 A / 2 k A ){ l-[l+4k A p 1 /p 0 A (c 0 A )2] 1 / 2 } (69a) 

PI =p 0 A Ui/(U 1 -u 1 ) (69b) 

This completely defines the state of the material behind shock wave I. For shock wave II we have 
Pl(U 2 + u i) = P 2 (U 2 +U 2 ) (70a) 

p 2 - p l = P l (U 2 +u 1 )(u 1 -u 2 ) (70b) 

and for shock wave III we have 

P 3 OJ 3 -U 3 ) = p 0 b (U 3 -u 4 ) (71a) 

p 3 - p 4 = Po B (U3 * u 4>( u 3 - u 4> (7 1 b) 

u 3 = c o B + kBu 3 (71c) 
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where P 4 has been replaced by p 0 ®, the ambient density of material B. In anticipation of the fact 
that material B will represent the module air, we assume stationary conditions but a non-zero 
pressure ahead of shock wave III. Enforcing pressure and velocity continuity at the A/B material 
interface (i.e. P 3 =P 2 a °d U 3 =U 2 ) reduces equations (71a-c) to 

p 3 (U 3 -U2) = Po B U 3 (72a) 

F2=P4 + Po B U 3 u 2 (72b) 

U 3 = c 0 B + k^u 2 (72c) 

Equations (70a, b) and equations (72a-c) form a system of 5 equations in 6 unknowns (P 2 , 
P 2 ,U 2 ,U 2 ,P 3 ,U 3 ). The elimination of one unknown is illustrated graphically in Figure 3. In 
Figure 3, Curves A and B are the Hugoniots of materials A and B, and Curve A' is the Hugoniot 
of materials A reflected about point C which denotes the initial shocked state in material A (i.e. 
prior to the passage of the reflected wave). The shocked state of material B must lie at the 
intersection of its Hugoniot (Curve B) and the reflected Hugoniot for material A (Curve A'). This 
state is denoted by point D. The particle velocity corresponding to point D is the interface 
velocity U 2 =U 3 while the pressure corresponding to point D is the interface pressure P 2 = P 3 - 

Knowing that curve A' is the reflection of Curve A, that is, it passes throughout the points 
(u p =U|,P=Pi), (up=2ui,P=0), and (up=0,P=2p o A u \ (c 0 A +2k A u 1 )), allows us to obtain the 
following functional form for Curve A' in P-Up space: 

P A' = 2 Po Au 1 (c 0 A+ 2 k A u 1 ) - P( Ac 0 A +4kAu j -k A u p )u p (73) 

Thus, when we set P^' equal to the functional form of Curve B in P-Up space, we have an 
equation for the particle velocity that corresponds to point D. Solving for this particle velocity 
yields: 
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ud = (P - VP 2 - 4ay ) / 2a 


(74) 


where 


a = p 0 AkA- Po B k B (75a) 

P = Po^o A + Po B °o B + 4p 0 A kA u j (75b) 

Y = 2p o A u 1 (c 0 A + 2k A Ul ) - P 4 (75c) 

in which uj, of course, is known. Setting U 2 =ud * n equations (70a,b) and (72a-c) allows us to 
solve for all the remaining quantities: 


U3 = c 0 B + k B U2 (76a) 

P3 =p 4 + Po B U3U2 (76b) 

P3 = po B U3/(U3-U2) (76c) 

U 2 = (P 1 -P2)/p 1 ( u 1 - u 2) - u 1 (77a) 

P 2 = Pl( u 2 +u lV( u 2 +u 2 ) (77b) 


5.3 Shock Wave Induced Pressure and Temperature Changes 

The analysis in Sections 5.1 and 5.2 is applied to the reflection and transmission of shock 
waves from the pressure wall/module air interface as follows. Clearly, material A is the pressure 
wall material (aluminum) and, as mentioned previously, material B is air. Since the axial velocities 
of the two secondary debris cloud components are not equal, the pressure wall impact event 
occurs at a time either before or after the impulsive loading event. As a result, it may be 
presumed that the impact event does not affect the impulsive loading event and vice versa. This 
implies that the total pressure wall rear surface velocity may be written as simply the sum of the 
two velocities U 3 and Up W as follows: 

u pa = u 3 +u pw (78) 
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This in turn is then also the total particle velocity in the module air immediately adjacent to the 
pressure wall rear surface. The velocity of the shock wave in the module air due to the combined 
effects of the impact and impulsive loadings is then given simply as 

u sa = c 0 B + k B u pa (79) 

while the initial pressure behind the shock front in the module air is given as 

p Ha = p 4 + Po Bu sa u pa (80) 

It is also noted that P3=P2 is the pressure level in the pressure wall material at which the release 
process for the pressure wall material terminates. For the purposes of this model, the shock wave 
is presumed to emanate spherically from the center of the pressure wall perforation. Naturally, 
the shock wave decays as it travels into the module environment. However, prior to its decay, the 
temperature increase in the immediate vicinity of the pressure wall is given by 

AT - (Epia - Eoa)/Cp a (81) 

where E oa = Cp a T 0 is the ambient module air energy (note: T 0 = 297-deg Kelvin; all temperatures 
are measured in degrees Kelvin), Cp a = 0.24 cal/gm°C is the specific heat of air, and Ejj a is given 
by 

B Ha ~ ^oa + ^ ( p Ha" l ’ p oa)( ^Poa* ^PHa) (82) 

In equation (82), p oa = 1 .22x10*3 gm/cm^ is the ambient density of air, P oa = 1 atm is the 
ambient module air pressure ahead of the shock wave, and 

PHa— u saPoa (83) 

u sa- u pa 

is the density of the shocked air. Thus, AT as calculated by equation (81) and AP=Pj^ a -P oa , 
where Pjj a is given by equation (80), define the change in pressure and temperature due to a 
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shock wave originating at the rear surface of the pressure wall in a pressurized module in the 
immediate vicinity of the pressure wall itself. The attenuation of the pressure increase as a 
function of distance away from the shock wave is discussed in the following section. 

5.4 Attenuation of Shock Wave Induced Pressures and Temperatures 

Figure 4 shows a sketch of a generic shock wave and its attenuation with distance. The 
decay with distance of the initial pressure and temperature induced by the shock wave in the 
module are obtained using the non-dimensional attenuation curve for normally incident spherical 
air blast waves [14], A closed-form analytical expression of the tabular P s ’ vs R’ data in [14] was 
obtained through a curve-fitting exercise for four regimes in R’ with the following results: 


0.01 <R’ <0.10 

P s ’ =3.803(R’)-1 2627 

0.10<R’ <1.0 

P s ’ = 0.476(R’)-22416 

1.0 <R’< 10.0 

P s ’ = 0.441 (R’)-l -3821 

10.0 < R’ < 1000.0 ... 

P s ’ = 0.241 (R’)- 11 120 


(84a) 

(84b) 

(84c) 

(84d) 


The initial pressure Ppj a is attenuated as a function of distance R away from the origination point 
of the shock wave by scaling values of P s ’ calculated using equations (84a-d) to Pjj a using the 


maximum value of P s \ i.e., the value of P s ’ at R’=0.01, as follows: 


p a( R ) = 


Ps’CR') 

P S '(R ,= o.oi) 


p Ha 


(85) 


where P a (R) refers to the pressure in the module at a distance R away from the origination point 
of the shock wave. The dimensional distance R is obtained from its non-dimensional counterpart 


R’ using the relationship in [14], that is, 


38 



R«R'(E/P oa ) l/3 (86) 

where ‘E’ refers to the amount of energy that has been transformed into a shock wave blast. This 
quantity can be easily calculated by performing an energy balance before and after the impact of 
the first secondary debris cloud component on the pressure wall. If energy losses due to light 
flash are again ignored, this energy balance is written as follows: 


F seconl _ rpsh&rl , , psh&rl v, , c initi 

E debcld _ K2 *- E pwall m wl + E debcld K l( ^ m P + m b)l + E s h w 


initial 

shwave 


(87) 


pop | 

where E<jeb C ld ' s t ^ ie tota * kinetic energy of the secondary debris cloud component considered in 


the impact analysis (i.e. it is the sum of the axial and expansion energies of the first secondary 
debris cloud component) prior to its impact on the pressure wall, E^Jy and E^j^j are the 

residua] energies per unit mass in the pressure wall and the secondary debris cloud following 
release from their respective shocked states, E^Jj^ is the initial energy of the shock wave (i.e. 
in the immediate vicinity of the pressure wall), and m w j is the mass of the pressure wall material 
subjected to shock heating and release as defined by equation (63a). The quantities E?J^y and 

| e i 

E debcld ^ ca l cu l ate d using the 1-D shock physics equations developed in Section 2.0. A 


constant K 2 < 1 has been inserted into equation (87) to allow for the possibility that not all of the 
secondary debris cloud and pressure wall materials will be loaded to a fully shocked state. The 
placement of K 2 outside the square brackets implies that the fractions of shocked and released 
material for the pressure wall and the secondary debris cloud material are assumed to be equal. 

Then, solving for yields 
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E = E 


initial _ F sec,ax P sec,exp rpS h&rl , , p sh&ri „ . mt « , Q9 . 

shwave “ E debcld E debcld K2 *- E pwall mwl E debcld Kl( - m P + ( 88 ) 


• sh& rl 


sh&rl 


We note that the quantity E as required by equation (86) and defined by equation (88) 
must be a positive number. Therefore, the fact that the right hand side is a difference between 
initial and expended energies is a signal that care must be taken to ensure that sufficient energy 
remains to allow shock wave creation and propagation. If, for example, the mass of the pressure 
wall material that is shocked and release is too large, it is conceivable that the right hand side of 
equation (88) may be a negative number. If this were the case, then the situation must be 
corrected either by decreasing K 2 or by decreasing m w [. The latter would be accomplished, for 
example, by decreasing the spread angle of the appropriate secondary debris cloud component. 

To find temperature decay with distance R, we use the following slightly modified forms 
of equations (81) and (82): 

AT = [E a (R)-E oa ]/Cp a (89) 

where P a (R) is given by equation (85), 

Ea(R) — E oa + ^ [P aO*-) + P oalt^Poa'^ Pa(®-)1 (90) 

is the energy behind a decaying shock wave as a function of the distance R, and p a (R) is the 
corresponding shocked air density. In this case, the density of the air behind a decaying shock 
wave is given as a function of R by the following modified form of equation (83): 


Pa(R) = Us(R)Poa 
u s (R)-u p (R) 


(91a) 


where 


u s (R) ~ c oa + ^a u p(R) 


(91b) 
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is the speed of the decaying shock wave as a function of R and Up is the corresponding particle 
velocity (also a function of the distance R). Finally, the particle velocity as a function of distance 
is found by substituting equation (91b) into the following alternate expression for P a (R): 


PaW Poa + Poa u sO*-) u p(R) 
Solving for Up(R) thus yields 

u n>\ - ' p oa c oa 

UP(R) ~ 2k aPoa 


(92) 


(93a) 


where 

A = (Poa c oa) 2 _ 4k a p oa [P oa*^ a(^)] (93b) 

5.4.1 Pressure and Temperature Values at a Specified Location 

Because the shock wave is presumed to propagate with spherical symmetry from the point 
on the pressure wall corresponding to the center of impact by the initial debris cloud, determining 
pressure and temperature values due to the passage of a decaying shock wave at a specific point 
of interest merely involves checking the radial distance from the decaying shock front to the point 
of interest. As the attenuation process progresses, if the point of interest is over-run by the shock 
wave, then a simple interpolation of the pressures and temperatures before and after the point of 
interest should yield a valid approximation of the pressure and temperature at that point. If the 
pressure jump across the shock wave has decayed to less than 1% of ambient and the point of 
interest has not yet been reached, then the attenuation process stops and we presume we are at 
ambient conditions. In this case, personnel at the point of interest will not feel any effect from the 
shock wave. 
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Figure 2a. Shock Wave Reflection and Transmission at a Material Interface — Unsteady Conditions 



U2=U3 

Figure 2b. Shock Wave Reflection and Transmission at a Material Interface - Steady Conditions 


■2 



Pa(R) 
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6.0 MODELING OF INDIVIDUAL INTERNAL EFFECTS: 

INTERNAL DEBRIS CLOUD PASSAGE 
6.1 Debris Cloud Induced Pressure and Temperature Changes 

In the event that module wall perforation results in a debris cloud entering the module 
interior, the temperature rise and pressure increase accompanying the passage of a shock wave 
created by the motion of such a debris cloud are determined as follows. 

First, the speed of the shock wave created by the motion of the internal debris cloud 
through the module air is taken to be equal to the leading edge velocity of the debris cloud. As 
discussed previously, this can be approximated by the particle velocity at the pressure wall/module 
air interface calculated using the impedance matching technique described previously, that is, 
u sdc =u 2 =u 3 While this may be presumed to be only valid in the immediate vicinity of the leading 
edge of the debris cloud near its centerline, applying this equality to the entire shock wave 
generated by the debris cloud should yield conservative results. Furthermore, we do not add the 
pressure wall velocity Up W to U3 to obtain the particle velocity at the pressure wall/module air 
interface. The velocity Up W is assumed to affect only the creation and movement of the internal 
shock wave that is created by pressure wall motion. 

Second, it is further assumed that the speed of the debris cloud and its accompanying 
shock wave do not decay until it reaches a certain distance away from the perforated pressure 
wall. If this distance is less than a module diameter, then the debris cloud effectively “runs out of 
steam” at this location and the velocity drops quickly to zero. It is at this distance (and point in 


45 



time) that the attenuation of the shock wave induced pressures and temperatures begin. The 
methods for calculating this distance is presented subsequently in Section 6.2.2. This assumption 
allows us to treat the shock wave velocity as essentially constant, which is consistent with the 
conservative nature of the model being developed. 

Figure 5 shows a sketch of a generic debris cloud moving through the interior of a module 
and the shock wave accompanying its motion. Following the definition of the shock wave 
velocity given in the preceding paragraphs, the particle velocity in the air induced by the shock 
wave is given by 

u pdc = ( u sdc " c oaV^a (94) 

from which the pressure behind the shock is given by 

PHa = p oa + Poa u sdc u pdc (95) 

where c oa =0.33 km/s, k a =1.07 and as before P oa =l atm. The shocked density is again given by 

PH a _ u sdcPoa (96) 

u sdc ~ u pdc 

As before, the energy associated with the passage of the shock wave induced by the internal 
debris cloud is 

E Ha " E oa + \ (PHa +p oa)( 1 /Poa* 1 /PHa) (97) 

and the associated temperature change is again 

AT = (E H a-E 0 a)/C pa (98) 

Thus, AT as calculated by equation (98) and AP=Pjj a -P oa as calculated in equation (95) define 
the change in pressure and temperature due to a internal debris cloud traveling within a 
pressurized module in the immediate vicinity of the debris cloud itself. The attenuation of the 
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pressure increase as a function of distance away from the debris cloud is discussed in the 
following section. 

6.2 Attenuation of Debris Cloud Induced Pressures 

As the debris cloud moves through the module interior, it continuously sets up a shock 
wave in the module air. Thus, attenuation of the shock wave so created can occur only when 
there is no longer any debris cloud motion. In keeping with the assumption of constant debris 
cloud velocity, it is presumed that debris cloud motion ceases either when the debris cloud 
impacts an obstacle such as the far module wall or when the drag force due to the presence of air 
within the module exhausts the axial kinetic energy of the debris cloud, whichever comes first. It 
is at this instant that the attenuation of the shock wave is begun. Figure 6 shows a sketch of this 
scenario and the geometric parameters required to perform this calculation. 

As can be seen in Figure 6, the distance parameter for the attenuation of the shock wave is 
T), not R as in Section 5.1.3. In fact, the decay of the pressure and temperature near the shock 
wave generated by the motion of an internal debris cloud in the module as a function of the 
distance r\ is actually obtained using the same procedure as that which was developed for the 
attenuation of the initial shock wave induced pressure and temperature. The only notable 
differences in the procedure are: 1) the quantity R in equation (86) is replaced by r| and refers to 
the distance away from the shock wave along the direction of travel (i.e. in the x-direction) of the 
internal debris cloud as shown in Figure 6; and, 2) the quantity E in equation (86) now refers to 
the total kinetic energy of the internal debris cloud (i.e. the sum of the axial motion and expansion 
energies). Thus, it is presumed for the purposes of this model that the decaying shock wave 
keeps its shape and direction of motion during decay; only pressure and temperature magnitudes 
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associated with the shock wave decrease with distance once attenuation has begun. 

The energy quantity required by equation (86) is found by applying momentum and energy 
conservation in a manner similar to that in Section 3.2.1. First, we conserve momentum before 
and after the impact of the first component of the secondary debris cloud on the pressure wall. 
Thus, we have 

Ki(m p +mt,)v^£ 1 =[Ki(m p +mb) + m w h]u mc (99) 

where v^c^ is the axial velocity of the secondary debris cloud first component, m w h is the mass 
of the pressure wall hole-out material, and u mc is the axial velocity of the internal debris cloud. 
We note that in equation (99) it has been assumed that the entire mass of the secondary debris 
cloud component used in the impact event analysis has been assumed to enter the module as part 
of the internal debris cloud. Second, once u mc is known the kinetic energy associated with the 
axial motion of the internal debris cloud can be written as 


E debcld ~ 1 ^ m P + ^ + m wh ] u mc / 2 (100) 

Next, energy conservation is applied before and after the impact of the secondary debris 
cloud on the pressure wall as follows: 


usee fax , pSec l,exp 
b debcld b debcld 



( 101 ) 


If we let u eX p denote the expansion velocity of the internal debris cloud, then we can write 

E deb did = t K 1 ( m p + m b ) + m w h l^mc / 2 (102) 

Substituting equations (102) into equation (101) and solving for u exp yields 
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“exp + EdS^ld^ _ ( K 2[ E ^S] 1 m w l + E debdd K l( m p + m b)l + E debdd ^ 1 ^ K l( m p + m b) + m wh ]/ 2} 

(103) 

Then, the energy quantity required by equation (86) is given by 


P_ P int,ax ,pint,exp 
il_ ^debcld + tl debcld 


(104) 


where and ^debcW are §‘ ven by equations (100) and (102), respectively. 

We note that the expression under the radical in equation (103) again is in the form of a 
difference between two quantities. Thus, there again exists a potential for spurious results if 
sufficient attention is not paid to, as before, the proper modeling of debris cloud velocities and, in 
this particular case, to the mass of the pressure wall hole-out material. The latter issue is 
addressed in the next sub-section. 

6.2.1 Pressure Wall Hole-Out Mass 

A question now arises as to what is the value of m w h, the mass of the pressure wall hole- 
out material. Clearly, its exact value is a function of material properties as well as geometry and 
the conditions of the initial projectile impact. However, previous high speed impact testing of 
aluminum dual-wall structures at impact velocities between 6 and 8 km/sec has revealed an the 
following two features of the nature of pressure wall perforation in dual-wall systems under 
hypervelocity impact. 

1) If an inner bumper is not placed between the bumper and the pressure wall, then the 
post-perforation pressure wall is relatively flat and its perforation consists of a single, well- 
defined, jagged hole in the center of the pressure wall damage area. In this case, a significant 
amount of mass is in fact removed from the pressure wall itself. 

2) If a relatively light inner bumper is placed between the bumper and the pressure wall 
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(e.g. a Type A inner bumper), then a pressure wall perforation again consists of a single, well- 
defined jagged hole in the center of the pressure wall damage area. In this case, a significant 
amount of mass is removed from the pressure wall. 

3) If a relatively heavy metallic or multi-material composite panel (e.g. a Type B inner 
bumper) is placed between the bumper and the pressure wall, then the perforated pressure wall is 
severely petalled. For impact velocities in excess of 6.5 km/s and projectile diameters greater than 
1 cm this results in huge through-holes. However, observation of the petals themselves has 
shown them to be mostly sharp, and not at all stubby or short. This implies that if the petals were 
to be bent back to their original flat positions, they would very nearly meet at the center of 
impact. In these cases, therefore, almost no pressure wall mass is lost despite the inordinately 
large holes created. 

Since the dual-wall construction to be used on the International Space Station is most 
likely to include some sort of thermal insulation and/or an inner bumper between the bumper and 
the pressure wall, it is assumed for the purposes of this study that the amount of pressure wall 
hole-out mass can be written as a fraction of the mass of pressure wall material that is subjected to 
shock heating and release by the impact of the secondary debris cloud, that is, 

m wh “ K 3( K 2 m wl) (105) 

where K 3 < 1 . For a Type A inner bumper, K 3 can be as large as 0.5 or 0.8 while for a Type B 
inner bumper, K 3 can be as small as 0. 1 or 0 . 2 . 

As mentioned in the preceding section, an improper choice of K 3 could have an adverse on 
the resulting internal debris cloud expansion velocity value. For example, if K 3 is too small (i.e. 
we are underestimating the mass of the pressure wall hole-out material), then according to 
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equation ( 99 ), u mc will be too large and it is conceivable that the radical in equation ( 99 ) may 
yield an imaginary value, which is clearly not acceptable. Naturally, there is more than one value 
of K3 that will yield expansion velocity values that are physically appropriate. However, the one 
that is ultimately most correct is the one that yields pressure and temperature values that are most 
in line with experimental results. 

6.2.2. Onset of Pressure and Temperature Attenuation 

Referring to Figure 6, the distance away from the pressure wall at which the attenuation of 
the shock wave induced pressures and temperatures begins is denoted by ‘H\ This distance is 
calculated by balancing the axial kinetic energy of the debris cloud and the work done by the drag 
force due to the presence of the module air as the debris cloud moves through the module. 
Naturally, if the value so calculated exceeds, for example, the distance from the impact site on the 
pressure wall to the nearest obstacle within the module (such as, for example, the far module 
wall), then the distance H would simply be the distance to the nearest obstacle. 

Thus, we write 

E debdd= F DH (106a) 

where is the axial kinetic energy of the internal debris cloud and Fj) is the drag force 

due to the module air. Elementary fluid dynamics provides the following expression to calculate 
the drag force: 

FD = PaCDA p v2/2 ( 106 b) 

where p a =1.22 kg/m^ is the density of air, Cj) is the drag coefficient, Ap is the presented area of 
the moving object, and V is the velocity of motion through the air. The drag coefficient is 
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typically a function of the Reynolds number Re which itself is given by 


Re p a VDp/[i a (107) 

where D p is the presented diameter and |i a =1.78xl0' 5 kg/m-sec is the viscosity of air. Assuming 
a nominal debris cloud diameter of 10 cm (i.e. Dp=0. 1 m) and a nominal debris cloud axial 
velocity of 1 km/sec (i.e. V=1000 m/sec) and substituting these values into equation (107) yields a 
Reynolds number in excess of 10^. In this Reynolds number regime, the drag coefficient curve 
for a sphere is relatively flat with a value between 0. 1 and 0.3. The specific value of Cj) used in 
this model will be discussed in a later section. 

To solve for the distance H, we note that Ap=7tR < j c 2. In addition, at the onset of 
attenuation (see again Figure 6), the radius of the debris cloud R^ c can obtained using the 
following relationships between debris cloud spread half-angle, the axial and transverse debris 
cloud velocities, and the geometry involved: 


tan^dc = 


sin<t>d c = 


Thus, we have 


4»dc = tan 


u exp 


u 


me 


K-dc 


H-R dc 


(108a) 


(108b) 


-1 


^ u exp^ 
l u mc> 


(109a) 


Rdc = Hf*, f* = Mn .*? c (109b, c) 

y 1 + sin<J>dc 

Substituting equations (106b) and (109b) into equation (106a), we obtain 

E debdd = * c DPa f ^ u mc H3 1 2 (1 10 ) 
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where V, the debris cloud characteristic velocity, has been replaced by u mc . Solving for H yields 


H = 


1 


pint, ax 
11 deb cl d 


^DPafj^mc / 2 


( 111 ) 


6.2.3 Pressure and Temperature Values at a Specified Location 

In order to determine the pressure and temperature values at a specified location within 
the module due to the passage of the internal debris cloud, several preliminary issues must be 
addressed. These issues are discussed sequentially, and are followed by a development of the 
equations required to obtain the desired information. 

First, if the point of interest is within the cone angle formed the motion and expansion of 
the internal debris cloud (region ABC in Figure 7), then personnel at the point of interest will 
probably be struck by the debris cloud itself; hence, the issue of pressure and temperature 
becomes secondary to the major concern of surviving the debris cloud impact. 

Second, if the point of interest is not within the cone angle characterizing the spread of the 
debris cloud, but is within the volume of air that is traversed by the shock wave induced by the 
motion of the debris cloud (region CAD in Figure 7), then personnel at the point of interest will 
feel the full pressure and temperature associated with the passage of the shock wave itself. 

Hence, it is not necessary to perform any attenuation of the shock wave induced pressure and 
temperature increases for such points of interest. 

Third, although the shock wave is presumed not to fan out during decay, its growth in the 
lateral direction (i.e. perpendicular to the axial direction of its motion) can still be observed as 
shown in Figure 7. If the point of interest is outside region HGEF, i.e. the cylindrical volume of 
module air whose base corresponds to the maximum extent of the attenuating shock wave, then 
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personnel at the point of interest will not feel any of the effects of the debris cloud shock wave. 

Fourth, once the pressure jump across the decaying shock wave (i.e. in the x-direction) is 
less than 1% of ambient, it is assumed that we are in fact at ambient conditions and the 
attenuation process stops. Any point of interest that lies “ahead of’ the shock wave in such a 
state will not experience any of the effects of the debris cloud or the shock wave it creates while a 
point of interest “behind” it will. Naturally, the attenuation process must first be performed to 
identify the position of this decayed shock wave, after which it can then be determined whether or 
not the point of interest lies “behind” it or “ahead of’ it. If during the attenuation process the 
point of interest is traversed by the decaying shock wave, then the attenuation process stops and 
the pressure and temperature at the point of interest are available. 

6.2.4 Definition of Shock Wave Position 

Based on the discussion in the preceding section, it is apparent that a functional form of 
the shock wave at the time attenuation begins and throughout the attenuation process is required. 
We begin by considering a shock wave in two dimensions and then extend the results to three 
dimensions. The equation itself derived from the fact that the curve describing the shock wave 
must satisfy at least the following two conditions. First, far from the leading edge of the shock 
wave, the angle defining the slope of the curve is given by 

p = sin- 1 (l/M) (112) 

where M=u sc j c /c oa is the shock wave mach number. Second, at the leading edge, the curve is 
perpendicular to the direction of travel. Although there are a number of curves that satisfy these 
conditions, the simplest curve that does so is the hyperbola. The specific form of the hyperbola 
that is appropriate for shock wave characterization is obtained as follows. 
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Consider the generic two-dimensional hyperbola shown in Figure 8. We recall the 


following properties of such a hyperbola in Cartesian coordinates: 

Equation : (x/a) 2 - (y/b) 2 = 1 

Asymptote Slope : m = b/a 

Distance From Origin to Vertex : a 

Distance From Origin to Focus : c 

Relationship Between c & a : 

Next, we note that equation (112) implies that the slope of the curve can also be written as 
1 


(113a) 

(113b) 


c 2 — a 2 + fc>2 


(113c) 


m = 


Vm 2 - 


(114) 


1 


Equating equations (1 13b) and (113d) allows us to solve for ‘b’ in terms ‘a’, to wit. 


b = 


a/m 2 - 


(115) 


Finally, as shown in Figure 8, we place the center-of-mass of the debris cloud at the focus of the 
hyperbola. In this case, we note that the leading edge of the shock wave is located at some 
distance A ahead of the debris cloud, which has some characteristic radius Rd c . This distance A is 
obtained from elementary shock physics theory (see, e.g., [15]) and is given by 

A = 0. 8 Rdc(Poa / PHa) ( 1 1 6 ) 

where pjj a is given by equation (83) and R<j c is given by equation (109b,c). 

We note in passing that Region ABC as discussed in Section 6.2 can now also be defined. 
Specifically, Region ABC in Figure 7 is defined by tan4),j c =u eX p/u mc . If we assume axial 
symmetry with respect to the x-axis (i .e. the longitudinal axis of the debris cloud), then to 
determine whether or not the point of interest, whose coordinates are given by (xp 0 j, yp 0 i> ZpoiX 
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is within Region ABC, we note that the equation of the circle in the y-z plane with its center at 
x=x poi lying on the cone subtended by the semi-vertex angle 4> dc is given by 


y 2 + z 2 = (x po jtan<t> dc ) 2 = R Xpoi 2 
If 


(117) 


ypoi 2 + ^oi 2 < R Xp 0 j 2 (118) 

then the point of interest is within Region ABC and personnel at that point will probably be struck 
by the expanding debris cloud. However, if 


ypoi 2 + z poi 2 > R x po , 2 


(119) 


then the point of interest is outside Region ABC. 

Once the center-of-mass of the internal debris cloud is placed at the focus of the 
hyperbola, it follows that the distance from the origin to the focus can also be written as 

c = a + A + R dc (120) 

Substituting for ‘b’ in terms of a according to equation (1 15) into equation (113c), solving for ‘c’, 
equating the result to equation (120), and solving for ‘a’ finally yields 


A + R dc 
E_ M 


( 121 ) 


yl M 2 -1 


- 1 


Since a and b are now both determined, the equation for the hyperbola approximating the shape of 
the shock wave created by a moving internal debris cloud is also completely determined. To be 
applicable to the two-dimensional coordinate system shown in Figure 8, the x-coordinate in 
equation (113a) must be offset by an amount h, that is, the equation for the two-dimensional 
hyperbola is now written as 
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x — h 


■> 


( 122 ) 


The quantity h is obtained by noting that the origin in Figure 7 is shifted over an amount H + A + 
a when compared to the origin in Figure 8; hence, 

h = H + A + a (123) 

In three dimensions, making use of the axial symmetry about the x-axis, equation (123) becomes 


x-h) 2 fy'i 2 (V 2 


= 1 


(124) 


Thus, the position of the shock wave at the onset of attenuation is now defined, and 
Region CAD as discussed in Section 6.2 can now also be defined. Specifically, the y- and z- 
coordinates (yp W , Zp W ) of the circle on the pressure wall with its center at x=0 and passing 
through Point D in Figure 7 is obtained simply by setting x=0 in equation (124), to wit. 


^y pw ^ 


v D J 


+ 


^pw V 


V D J 


H + A + a^ 


-1 


(125) 


To determine whether or not the point of interest, whose coordinates are given by (xp 0 j, yp 0 j, 
Zp 0 i), is within a region like the two-dimensional region CAD, we note that the equation of the 
circular cross-section of the three-dimensional hyperboloid at x= Xp Q j is given by 


0 


if*© 2 - 


x poi — (H + A + a) 


i2 


->=( R x P oi y 


(126) 


If it is true that 






z poi 
V b 




< (R x • Y 
v x poi J 


(127) 
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then the point of interest is within a region like Region CAD and the full effects of the shock wave 
are felt at the point of interest. However, if 


'y P oP 2 


V 


b J 


+ 


f z A 2 
z poi 


b J 


>(R X • y 

v x poi y 


(128) 


then the point of interest is outside a region like Region CAD and only the effects of an attenuated 
shock wave (or not at all) are experienced at the point of interest. 

To define the curve corresponding to the position of the shock wave at the completion of 
the attenuation process (i.e. pressures and temperatures have decayed to ambient conditions), we 
recall that we have assumed that the shock wave merely translates axially without fanning out 
during the attenuation process. In this manner, the x-coordinate of the origin of the hyperbola 
shifts by an amount rij, where r|j is the distance away from the full-strength shock wave at the i-th 
attenuation step. Thus, when r|j = 11^5, if we let (y pw ,amb> z pw,amb) denote the y- and z- 
coordinates of the circle on the pressure wall that has its center at x=0 and that lies on the 
“ambient” hyperboloid, then equation (125) becomes 


ypw,amb 


-pw,amb 


fH+A+a+ 


Ramb 


) 2 - 


1 


(129) 


As before, to determine whether personnel at a point of interest defined by the coordinates (xp 0 j, 
y po i, Zpoi) experience the effects of an attenuated shock wave or not at all, we first calculate the 
quantity 


< R Xpoi )2 “ 


x poi “(H + A + a + rii) 


i2 


-1 


(130) 


where r|j is the distance away from the full strength shock wave at the i-th attenuation step. If 
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then it is outside this zone of influence at the current attenuation step. If however, equation (132) 
is true for the entire attenuation procedure, then no shock effects are experienced at the point of 
interest. Likewise, if the pressure jump across the shock wave has decayed to less than 1% of 
ambient and the point of interest has not yet been reached, then the attenuation process stops and 
we again presume that we are at ambient conditions and that personnel at the point of interest do 
not experience any shock wave effects. 

Finally, we note that as the attenuation process progresses, if the point of interest is over- 
run by the shock wave, then an interpolation procedure is employed to determine the pressure and 
temperature at the point of interest. However, in the case of a shock wave due to the passage of 
an internal debris cloud, the interpolation process is somewhat more complicated than the 
elementary one employed for the spherically expanding shock wave in Section 5.3.1. In this case, 
since the attenuation proceeds in the x-direction, we must first calculate, xf s , the x-coordinate of 
the point on the hyperboloid corresponding to the y- and z-coordinates of the point of interest. 
This is done by solving for x in equation (89c). Once this value is obtained, then the interpolation 
is performed between the pressures and temperatures at the i-th and (i-l)-th attenuation steps, 
where the x-coordinate of the point of interest has been over-run in the i-th step, by referencing 
x poi t0 xfs+rii-l 
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Figure 5. Generic Internal Debris Cloud With Accompanying Shock Wave 
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Figure 6. Onset of Internal Debris Cloud Shock Wave Attenuation 
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Figure 7. Pressure and Temperature Values at a Specific Location Due to a 
Debris Cloud Induced Shock Wave 
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Figure 8. Two-Dimensional Hyperbola as a Model for a Debris Cloud Induced Shock Wave 
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7.0 SUPERPOSITION OF INDIVIDUAL INTERNAL EFFECTS 
In modeling the interaction of the phenomena discussed in Chapters 5 and 6, one of the 
key issues is the partitioning of the kinetic energy of the initial debris cloud when both a shock 
wave and an internal debris cloud are created and propagate into the module interior. Once this 
partitioning is performed, the analysis required to model the interaction of an initial shock wave 
and an internal debris cloud on the pressure and temperature at a specific point within the module 
proceeds in a relatively straightforward manner. 

7.1 Energy Partitioning 

To begin, we write the following energy balance between the energy of the system before 
and after the impact of the initial debris cloud on the pressure wall: 


psec fax 
^debcld 


+ E 


secfexp _ 
debcld “ 



A A 

pint, ax pint, exp 
11 debcld + 11 debcld 


(133) 

We note that equation (133) is very similar to equations (87) and (101) in Section 5.3 and 6.2. 
The ‘a* symbol over the initial shock wave and internal debris cloud energy terms in equation 
(133) denoted modified values calculated when a shock wave and an internal debris cloud are 
both present. 

As in Section 6.2, the kinetic energy associated with the axial motion of the debris cloud is 
written as 


A 


E d"bdd + m b) + "' w h](u mc .l 2 


(134) 
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A 

where u mc is the axial velocity of the of the internal debris cloud. If the momentum of the 
shock wave is ignored in the conservation of momentum before and after the impact of the 
secondary debris cloud first component on the pressure wall, then this quantity is found by solving 
equation (99) with the following result: 


A 

u mc 


K l(m p +m b )vmc 1 
Ki(m p +m b ) + m w h 


(135) 


where v^c ^ is the momentum conservation based characteristic velocity of the first component 
of the secondary debris cloud. We note that, under the assumptions made, the value of u mc does 
not depend on whether or not an initial shock is created in addition to the internal debris cloud. 
However, the presence of the initial shock wave does become a factor when we calculate the 
kinetic energy associated with the expansion of the internal debris cloud. As in Section 6.2, this 
energy is written as 

A 

E debdd = ^ K l< m p + m b ) + m wh ](u exp ) 2 ( 136 > 

where, in this case, u exp is still unknown. As such, equation (101) as written in Section 6.2 
would have two unknowns and could not be solved directly for u exp . Thus, another equation 
relating u mc and u exp is required. This equation is obtained by noting that the expansion velocity 
of a debris cloud, at least near the front of the debris cloud on the trajectory centerline, is also the 
difference between the leading edge velocity and the center-of-mass velocity of the debris cloud. 
If we presume this to be true for the entire debris cloud, we can then write 

u exp = u le ~ u mc (137) 
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where uj e is the leading edge velocity of the internal debris cloud. 

We note that if we simply took uj e to be equal to the interface particle velocity obtained 
from the impedance matching calculations as in Section 6. 1, i.e. u| e =U2 =u 3 ^sdc' then the 
interaction between the creation of the initial shock wave and the motion of the internal debris 
cloud would be minimal. Specifically, if we claim that ui e or u S( j c is already known, then all that 
remains to be done is to solve for the energy associated with the initial shock wave using equation 
(133). Once this value is obtained, we then obtain peak pressures and temperatures and their 
attenuated values using the techniques presented in Chapters 5 and 6. 

However, the presumption that ui e =u sc j c tacitly implies that the radial expansion as well as 
the axial motion of the internal debris cloud is independent of the presence of an initial shock 
wave in addition to the debris cloud. After all, we are using values obtained in Chapter 6 to 
characterize its motion and expansion, and those values were obtained assuming the presence of 
only an internal debris cloud. If we wish more of an interaction to exist between the axial motion 
and radial expansion of the internal debris cloud and the creation of the initial shock wave, then 
the leading edge velocity of the internal debris cloud should be modified somehow from the value 
calculated in Chapter 6. Specifically, the modified value should reflect the fact that some portion 
of the kinetic energy associated with the expansion of the internal debris cloud has been absorbed 
by the creation of the initial shock wave and vice versa. 

To this end, we propose to use the following relationship between debris cloud leading 
edge velocities and initial shock wave energies: 


( u sdc) _ (usdc) 2 

a ~ £ initial 

£ initial sh wave 

11 sh wave 


(138) 
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In equation (138), E^wave an< * u sdc re f er t0 the initial shock wave energy and internal debris 


cloud leading edge velocity that were computed independently in Chapters 5 and 6 while 
corresponding quantities with a ‘a‘ symbol denote a coupled initial shock wave energy and 
internal debris cloud leading edge velocity. Then, according to equation (138), we can write 


A 

^initial 
11 sh wave 


_ jjinitial 
shwave 


A 

u sdc 

u sdc 


(139) 


Substituting for 


A 

• initial 
'shwave 


in 


equation (133) according to equation (139) and rearranging terms 


A 

yields the following equation for u sc j c : 


u 


2 

AE 


A 2 ^ 2 

E o( u sdc) +( u sdc- u mc) 


(140) 


where 


ginitial 

shwave 

[ K l(m p + mtj) + m wh]Ug £ j c / 2 

u 2 = AE 

^ [*l( m p +mb) + m w h]/2 


(141) 


(142) 


and 

48 - Edebdd 1 + E Sd'd XP - l K 2l 8 pwaU m wl + E ffid*l< m P + m b)]+(Kl<">p + m b ) + / 2) 

(143) 

Solving equation (143) for the leading edge velocity of the internal debris cloud in the presence of 
an initial shock wave yields 
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(144) 


A 

u sdc 


u 


me 


+ -J u mc (Eq + 0( u mc u af.) 


E 0 + 1 


A 

Then, once u s d c *s known, the expansion velocity of the internal debris cloud is found using 

A 

equation (137) with ui e =u sc j c (which is now also the velocity of the shock wave created by the 

motion of the internal debris cloud), and the energy associated with the creation of the initial 
shock wave is found directly using equation (139). Once again, we note that the expression for 

A 

u sdc contains a radical under which is an expression that is in the form of a difference between 


A 

two quantities. In addition, we note that substituting the value of u S( j c as obtained using 

equation (144) into equation (137) may yield a negative expansion velocity. In both cases, 
adjusting the value of u mc (or, rather, m w j 1 ) would improve the results. 

7.2 Initial Shock Wave Pressures and Temperatures 
7.2.1 Peak Pressures and Temperatures 

In the preceding section, a procedure was presented to calculate the leading edge velocity 
of an internal debris cloud (which is then also the velocity of the shock wave created by the 
motion of this debris cloud) when an initial shock wave was also present. As was stated during 
the course of that development, the leading edge velocity so calculated is directly related to the 
particle velocity at the module wall/air interface when both an initial shock wave and an internal 
debris cloud are created. Applying this information to the concepts and equations in Chapter 5 
allows us to write the total particle velocity in the air immediately adjacent to the module wall as 

A A 

u pa =u sdc + u pw (145) 
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from which the initial shock wave velocity is given by: 


A A 

u sa = c 0 + ku 


pa 


( 146 ) 


A A 


Now that Up a and u sa are known, the peak pressure and temperature induced by the initial 
shock wave are readily obtained using the equations in Section 5.2 as follows: 


A A A 

^Ha = P oa + Poa u sa u 


pa 


(147) 


and 


AT = ( E Ha- £ oa) / C 


pa 


(148) 


where 


E Ha - E oa + 2 ( E Ha + E oa / Poa ~ 1 / PHa) 


(149) 


and 


PHa=^f- 

u sa u pa 


(150) 


7.2.2 Attenuated Pressures and Temperatures 

The attenuation of the peak pressure and temperature defined in the preceding section 
follows a procedure virtually identical to that developed in Section 5.4. The only difference is 
that we now use the a quantities where appropriate to denote the shock wave pressures, 
velocities, etc. that are calculated considering the presence of an internal debris cloud as well. 
Specifically, we now have 
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(151) 


Pa(R) = 


V(R) 

P S '(R' = 0.01) 


PHa 


A 


R = R'(E/P 0a )' /3 .E = E'"“f ve 

AT = [Ea^R)- E oa ] / C pa 

Ea(R) = Eoa + 2^[Pa(R) + PoalH f Poa - 1 ^ Pa(R)] 


(152a,b) 

(153) 

(154) 


A 


Pa(R)= “ s(R)P “ 

(155) 

u s (R)-u p (R) 


A A 


u s (R) = c oa +k a u p (R) 

(156) 

/a 

ii 'rp'i- ~Poa c oa +VD 
Up(R> - 2k aPoa 

(157) 

D = (Poa c oa) - 4k a p oa [P oa - Pa(R)] 

(158) 


7.2.3 Pressure and Temperature Values at a Specified Location 

The procedure for calculating the pressure and temperature increase at a specified location 
within the module for an initial shock wave that has been created in conjunction with an internal 
debris cloud is identical to that already presented in Section 5.4.1. 

7.3 Internal Debris Cloud Pressures and Temperatures 
7.3.1 Peak Pressures and Temperatures 

As was stated in Section 7.1, the leading edge velocity is also the velocity of the shock 
wave created by the motion of the internal debris cloud. Following the procedure in Section 6.1, 
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the particle velocity in the air induced by the motion of the shock wave accompanying the internal 
debris cloud is given by 

A A 

u pdc = ( u sdc~ c oa) / ^a (159) 

from which the peak pressure behind the shock wave and the corresponding peak temperature 
increase are given by 

A A A 

l*Ha = ^oa + Poa u sdc u pdc (160) 

A A 

AT = (Ejja - E 0 a) I Cpa (161) 

where 

^Ha = Eoa + 2 ^(^Ha + l*oa)(l ! Poa - 1 / PHa) (162) 

and 


A 

PHa = 


A 

u sdc Poa 

A A 


u sdc u pdc 


(163) 


7.3.2 Attenuated Pressures and Temperatures 

The attenuation of the peak pressure and temperature defined in the preceding section 
follows a procedure virtually identical to that developed in Sections 5.4, 6.2, and 7.2.2. Once 
again, referring to Figure 6, we use rj and not R as the distance parameter for the attenuation 
process, and, as in Section 7.2.2, we use the *a‘ quantities where appropriate to denote the shock 
wave pressures, velocities, etc. that are calculated considering the presence of an internal debris 
cloud as well as an initial shock wave. Finally, as in Section 6.2, the quantity E in equations (86) 
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and (152a) now refers to the kinetic energy associated with the motion of the internal debris cloud 
and is given by the sum of the energies in equations (134) and (136). 

7.3.3 Pressure and Temperature Values at a Specified Location 

The procedure for calculating the pressure and temperature increase at a specified location 
within the module for shock wave that has been created due to the motion of an internal debris 
cloud is nearly identical to that presented in Sections 6.2. The major difference lies in the 
definition of the internal debris cloud spread cone angle and the radius of the internal debris cloud. 
Specifically, we now have 

A 


A 

tan<|)dc = 


u sdc u mc 
u mc 


(164) 


and 


A 

x A R dc 
sm(J)dc = 

H-Rdc 


(165) 


from which 


A 


4>dc =tan 


-1 


A 

u sdc - u mc 


u 




me 


J 


(166) 


and 


A 

r A c _ Hsinfrdc 

A 

1 + sin <t>dc 

Following these revisions, the process to calculate the 
location due to the motion and attenuation of a shock 


(167) 

temperature and pressure at a specified 
wave created by an internal debris cloud 
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follows the same procedure as outlined in Sections 6.2. 1 and 6.2.2. 

7.4 Cumulative Pressure and Temperature Values at a Specified Location 

Once the pressure and temperature values at a specified location have been determined 
due to the passage of the initial shock wave and the shock created by the passage of the internal 
debris cloud, their cumulative effect at that point is obtained by merely adding the results 
together. 
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8.0 COMPARISON WITH EXPERIMENTAL RESULTS 


8.1 Test Configurations and Impact Parameters 

The predictions of the model developed herein are compared against the experimental 
results obtained from 4 high speed impact tests conducted at the University of Alabama in 
Huntsville Aerophysics Research Center. A detailed description of the test facility has been 
presented in Volume I of this two volume report. Table 2 below presents a summary of the 
geometric and impact parameters for the 4 impact tests. 

Table 2. Impact Parameters and System Geometries 


Test 

mm 

mm 

mm 

S 

MgM 

Inner 

No. 

■SSI 


RSI 

(cm) 

SI 

Bumper 

1 

1.27 

6.70 

0.16 

11.43 

0.48 


3 

1.59 

6.42 

0.16 

11.43 

0.48 

ES73X9 

5 

1.59 

6.58 

0.16 

11.43 

0.48 

ESI 

6 

1.59 

6.50 

0.16 

11.43 

0.48 

E353D 


The MLI in both cases consisted of 20 mylar layers. In the case of the Type B inner 
bumper, the multi-layer composite material panel consisted of 6 layers of Nextel AF62, 6 layers of 
Kevlar 710, and a 0.203 cm thick graphite/epoxy panel. Additional information regarding the 
composition of the Type B inner bumper can be found in Reference [16]. Finally, in all cases the 

initial ambient temperature behind the pressure wall was 24°C. 

Table 3 below presents a summary of the x-, y-, and z-coordinates of the various pressure 
(P) and temperature (T) sensors that were used to gather the pressure and temperature change 
data. It is these coordinates that are input into the model to produce the theoretical predictions. 
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The final column in Table 3 shows the distance from the origin of the coordinate system to the 


sensor location. 


Table 3. Sensor Coordinates 


Sensor 

Sensor 

Types 

Shot 

No. 

X 

(cm) 

y 

(cm) 

Z 

(cm) 

Rd 

(cm) 

IE 

P, T 

All 

18.0 

61.0 

0.0 

63.6 

2E 

P, T 

All 

119.0 

61.0 

0.0 

133.7 

3E 

P, T 

All 

203.0 

61.0 

0.0 

212.0 

1W 

P, T 

All 

18.0 

-61.0 

0.0 

63.6 

2W 

P, T 

All 

119.0 

-61.0 

0.0 

133.7 

3W 

P 

All 

203.0 

-61.0 

0.0 

212.0 

3W 

T 

1 

203.0 

-61.0 

0.0 

212.0 

3W 

T 

3,5,6 

203.0 

0.0 

0.0 

203.0 


As can be seen in Table 3, Pressure Sensors 1 W, 2W, and 3W were mirror images of 
Pressure Sensors IE, 2E, and 3E, respectively, with respect to the x-axis for all the tests while the 
symmetry of the temperature sensors was in effect for Test No. 1 only. For Test Nos. 3,5 and 6 , 
Temperature Sensor 3E was place in the same position as in Test No. 1 while Temperature 
Sensor 3W was placed directly on the x-axis at the corresponding distance. Therefore, in 
subsequent tables, when comparing model predictions to experimental results, the results for 
Sensors 1E,W, 2E,W, and 3E,W will be grouped together. 

8.2 Model Parameter Values 

As discussed in the preceding chapters, the results produced by the analytical model 
depend on the values of the empirical constants ‘n’, K 2 , and K 3 . Poorly chosen values could lead 

to either spurious results (e.g. negative or imaginary debris cloud expansion velocities) or 
unrealistic results (e g. debris clouds that area either too heavy or too light). In addition, an 
improper internal debris cloud drag coefficient value could result in debris clouds that stop either 
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too soon or too late. In the first case, the resulting pressures and temperatures at points far 
removed from the pressure wall would be far too low while in the second, the shock wave 
produced by the motion of the internal debris cloud could possibly engulf far removed points of 
interest and result in pressures and temperatures at those points that are unrealistically high. 
Naturally, a fairly substantial experimental database is required to determine specific values of the 
empirical constants and the drag coefficient as functions of projectile diameter, impact velocity, 
system geometry, etc. 

Before determining the actual values of the constants for the analytical model, a series of 
trial runs were performed to gain insight into the effects of the various constants on the 
predictions of the analytical model. These trial runs were performed at an impact velocity of 6.5 
km/sec, and a dual-wall system with a 0. 16 cm bumper, a 0.48 pressure wall, and a 1 1 .43 stand- 
off distance, and a Type A inner bumper; a drag coefficient of 0.2 was assumed for the purposes 
of calculating the stopping distance of the internal debris cloud. Results were obtained at 
coordinates corresponding to those of Sensor IE as given in Table 3 . 

Tables 4,5 and 6 present a summary of the output obtained from the trial runs that were 
performed to determine appropriate empirical constants. Table 4 addresses the effect of changing 
‘n’ while holding K 2 and K 3 constant, while Tables 5 and 6 show the effect of varying K 2 only and 

K 3 only, respectively. Because changes in K 2 and K 3 affect only the internal debris cloud, 

comparisons of secondary debris cloud characteristics are omitted in Tables 5 and 6 . Examination 
of these tables reveals the effects on model prediction of changing a given parameter while 
holding the others fixed. 



Table 4. Effect of Changing ‘n’ on Model Predictions (k 2=0.2, *3=0. 3) 


Parameter 


K 1 


e dc Cl(de *> 


m^ cl (gms) 


v mc 1 (km/s) 


vg^lcm/s) 


6^f 2 (deg) 


m“ c2 (gms) 


v^ 2 (kn,/s) 


v exp 2 (km/ s) 



Uga (km/s) 


E SW (J) 


e3?(deg) 


m dc (g ms ) 


u mc ( km/s ) 


u| e (km/s) 


Edc(J) 


^sw (psi) 


AP dc (psi) 


AP tot (psi) 


AT SW (OC) 


AT dc (OC) 


AT tot (OC) 


n=0.2 

n=1.0 

0.989 

0.945 

34.54 

34.54 

4.558 

4.355 

4.100 

4.080 

2.823 

2.810 

34.54 

34.54 

2.128 

2.128 

0.102 

0.525 


0.290 

0.289 

3.301 

3.291 

3.863 

3.852 


10.08 

11.01 

7.476 

7.13 

2.500 

2.496 

0.511 

0.503 

3.011 

2.999 

2.4x1 0 4 

2.4xl0 4 

0.90 

.86 

13.26 

9.58 

14.16 

10.44 

3.63 

3.47 

37.02 

28.42 

40.65 

31.89 



77 


























































For example, in Table 4, we see that increasing ‘n’ from 0.2 to 1.0 serves to decrease 
pressure and temperature increases due to both the shock wave and the internal debris cloud. 

This is due to the fact that an increase in ‘n’ results in a decrease in iq, which in turn implies that 

less primary debris cloud mass is passing through the inner bumper. Because of the manner in 
which the energy and momentum of the primary debris cloud are partitioned among the two 
secondary debris cloud components, this in turn results in a decrease in the characteristic 
velocities of the first component of the secondary debris cloud and an increase in the second 
component velocities. The decrease in the characteristic velocities of the first component are 
directly responsible for the drop in the shock wave induced pressure and temperatures increases. 
However, the drop in the debris cloud induced pressure and temperature increases is due to the 
fact that less mass is coming through the inner bumper and hence into the module itself. Because 
of the manner in which the creation of the module shock wave is presumed to interact with the 
creation of the internal debris cloud, the velocities of the internal debris cloud are only slightly 
lower for the larger value of ‘n\ Hence, the energy associated with the internal debris cloud in 
the second case (a larger value of ‘n’) is less than that in the first case (a smaller value of ‘n’), 
which in turn results in lower pressure and temperature increases. 

In Table 5 we see that increasing vq (i.e. increasing the amount of shocked and released 

secondary debris cloud and pressure wall material) decreases the amount of energy converted into 
shock wave formation which in turn decreases the pressure and temperature increases due to the 
module shock wave. However, increasing K 2 also serves to increase the mass of the internal 

debris cloud which decreases its axial velocity. This in turn increases the expansion velocity and 
the spread angle of the debris cloud. The final result is that the shock wave induced by the 
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internal debris cloud fans out more and affects a larger volume of module air. This could cause 
remote points to experience higher pressure and temperature increases than if the lower value of 
K2 were used. 


Table 5. Effect of Changing K2 on Model Predictions (n=0.5, K3=0.3) 


Parameter 

k 2 =0.2 

OO 

?■ 

u pa (km/s) 

3.301 

2.800 

u sa (km/s) 

2.862 

3.326 

E sw (J) 

2.8xl0 4 

1.5xl0 4 

C(deg) 

12.06 

50.39 

m dc (g ms ) 

7.399 

16.154 

u mc (^s) 

2.481 

1.137 

u exp (km/s) 

0.530 

1.374 

uj e (km/s) 

3.011 

2.510 

E dcP) 

2.4x1 0 4 

2.6xl0 4 

APsw (P s 0 

0.89 

0.45 

APdc (Psi) 

10.75 

68.31 

APtot (P si ) 

11.64 

68.76 

AT SW (°C) 

3.59 

1.86 

AT dc (°C) 

31.22 

155.40 

AT, ot (°C) 

34.81 

157.26 


In Table 6 we see that increasing K3 (i.e. the fraction of shock and released pressure wall 
material that becomes part of the internal debris cloud) naturally increases the internal debris 
cloud mass. This again decreases the internal debris cloud axial velocity, increases the expansion 
velocity and the debris cloud spread angle and has the same net result as that which was observed 
with an increase in K2- 
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Table 6. Effect of Changing K3 on Model Predictions (n=0.5, K2=0.2) 


Parameter 

k 3 =0.3 

K3=0.8 


3.301 

3.191 

u sa (km/s) 

2.862 

3.745 

E sw CD 

2 . 8 xl 0 4 

2 . 6 xl 0 4 

ejfcWeg) 

12.06 

43.16 


7.399 

12.263 

u mc Conk) 

2.481 

1.497 


0.530 

1.404 

uj e (km/s) 

3.011 

2.901 

EdcCD 

2.4xl0 4 

2 . 6 xl 0 4 

AE sw (psi) 

0.89 

0.79 

APdc (psi) 

10.75 

16.40 

AP to t (psi) 

11.64 

17.19 

AT SW (°C) 

3.59 

3.19 

AT dc (°C) 

31.22 

44.06 

> 

0 

/~N 

O 

O 

34.81 

47.25 


Based on the trends observed Tables 4, 5, and 6 , it appears that if, for example, in 
contradiction with observed pressure and temperature phenomena, the model predicts a remote 
point of interest to be engulfed by a full strength debris cloud induced shock wave, then it would 
be appropriate, for example, to decrease either K 2 or K 3 , or both. An alternative means of 

producing such a result would be to increase the drag coefficient of the internal debris cloud. This 
would cause the debris cloud to stop sooner so that remote points of interest would feel only 
attenuated pressures and temperatures. 

To determine the values of the constants ‘n’, < 2 , K 3 , and Cd for use in the initial 
comparisons between the predictions of the model with experimental results, another series of 
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runs was performed using the model with various constant combinations. Naturally, the choice of 
constant combinations was guided by the trends observed in the trial runs and the actual 
experimental values as recorded by the pressure and temperature sensors at the locations indicated 
in Table 3. As such, the runs were performed using the system geometries and impact velocities 
as specified in Table 2. Table 7 presents a summary of the initial empirical constant values as a 
function of inner bumper type and projectile diameter. These constants should produce 
reasonable results for the impact velocities, projectile diameters, materials, and system geometries 
that were tested as part of this study. 


Table 7. Initial Empirical Constant Values 


Empirical 

Constant 

Type A 
Inner Bumper 

TypeB 
Inner Bumper 

Projectile 

Diameter 

( cm ) 

1.27 

1.59 

1.59 

n 

0.50 

2.50 

1.10 

K 2 

0.19 

0.14 

0.23 

K 3 

0.30 

1.00 

0.88 

c D 

0.20 

0.30 

0.11 


Naturally, when a sufficiently large amount of test data becomes available, it should be 
possible to obtain functional expressions for each of the empirical constants in Table 7 in terms of 
system geometry, materials, impact conditions, etc. 

8.3 Comparison of Model Predictions with Experimental Results 

Tables 8 through 1 1 present a comparison of the experimental results and the predictions 
of the model developed herein for the sensor locations presented previously. It is noted that only 
initial peak experimental values (i.e. values which occur at times less than approximately 150 
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msec) are compared against the predictions of the analytical model. In its present stage of 
development, the model is simply not capable of predicting very late time effects which involve 
numerous shock wave interactions and reflections. As seen from the complete data presented in 
Volume I, the temperature at very late times often exceed by a factor of 2 or 3 the initial peak 
temperature increases. 

Tables 8 through 1 1 also show the times at which, the initial peak values occurred at the 
various sensors. These times are presented to determine whether or not complementary sensors 
(e.g. Sensors 1E,W) registered peak values at approximately the same time. Ideally, due to the 
symmetric nature of the impact event and the sensor placement (with the exception of 
Temperature Sensors 3E,W for Test Nos. 3, 5, and 6), peak values should occur (nearly) 
simultaneously at complementary sensors and should have (nearly) identical magnitudes. Any 
significant deviation from nearly identical peak values and their times of occurrence is an 
indication of an asymmetry in the propagation or composition of the internal debris cloud. 
Examination of the experimental results in Tables 8 through 1 1 reveals that the pressure and 
temperature change magnitudes recorded and their times of occurrence were fairly symmetrical. 
Naturally, some extreme deviations did exist in certain individual instances, but as a whole the 
recorded values were fairly symmetrical. 

Table 8. Pressure and Temperature Data Comparison (Basic Model) — Test No. 1 



Pressure Increase 

Temperature Increase j 

Sensor 

Time 

(msec) 

Experimental 
Result (psi) 

Model 

Prediction (psi) 

Time 

(msec) 

I 

Model 

Prediction (°C) 

IE 

1.12 

17 

15.2 

13.3 

18 

43.2 

1W 

1.07 

18 

15.2 

19.0 

27 

43.2 

2E 

1.93 

4 

0.5 

21.6 

23 

2.2 

I 2W 

2.15 

6 

0.5 

23.0 

20 

2.2 

3E 

2.33 

3 

0.2 

645.8 

25 

0.7 

3W 

2.00 

2 

0.2 

n/a 

n/a 

0.7 
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Table 9. Pressure and Temperature Data Comparison (Basic Model) -- Test No. 3 




Pressure Increase 

Temperature Increase 1 

Sensor 

Time 

(msec) 

Experimental 
Result (psi) 

Model 

Prediction (psi) 

Time 

(msec) 

Experimental 
Result (°C) 

Model 

Prediction (°C) 

IE 

0.86 

21 

24.3 

8.5 

104 

63.9 

1W 

1.00 

21 

24.3 

6.5 

118 

63.9 

2E 

1.33 

17 

0.7 

10.9 

139 

3.0 

2W 

1.50 

21 

0.7 

10.8 

94 

3.0 

3E 

2.13 

8 

0.2 

n/a 

n/a 

' 2.1 

3W 

2.83 

8 

0.2 

115.0 

40 

2.1 


Table 10. Pressure and Temperature Data Comparison (Basic Model) - Test No. 5 




Pressure Increase 

Temperature Increase 1 

Sensor 

Time 

(msec) 

Experimental 
Result (psi) 

Model 

Prediction (psi) 

Time 

(msec) 

Experimental 
Result (°C) 

Model 

Prediction (°C) 

IE 

0.95 

34 

34.2 

13.4 

64 

82.7 

1W 

0.80 

40 

34.2 

17.9 

54 

82.7 

2E 

2.13 

10 

0.0 

19.6 

34 

0.0 

2W 

1.93 

15 

0.0 

20.0 

31 

0.0 

3E 

4.07 

3 

I 0.0 

n/a 

n/a 

0.0 

3W 

3.90 

4 

0.0 

n/a 

n/a 

0.0 


Table 11. Pressure and Temperature Data Comparison (Basic Model) - Test No. 6 




Pressure Increase 

Temperature Increase 1 

Sensor 

Time (msec) 

Experimental 
Result (psi) 

Model 

Prediction (psi) 

Time (msec) 

Experimental 
Result (°C) 

Model 

Prediction (°C) 

IE 

1.33 

16 

24.6 

25.0 

20 

64.8 

1W 

1.06 

31 

24.6 

26.3 

65 

64.8 

2E 

1.71 

15 

0.8 

92.0 

52 

3.2 

2W 

1.54 

14 

0.8 

53.0 

63 

3.2 

3E 

2.65 

12 

0.2 

52.0 

44 

1.1 

3W 

2.66 

7 

0.2 

n/a 

n/a 

2.1 


An examination of the information in Tables 8 through 1 1 reveals that the model performs 
reasonably well in predicting the pressures experienced by Sensors 1E,W. However, with regard 
to the other sensors, the although the model correctly predicts substantial pressure decreases at 
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Sensors 2E,W and 3 E,W, the predicted values are significantly lower than the experimental 
values. It is almost as if the attenuation of the peak pressures occurs too rapidly with distance 
away from the shock wave. Based on the significant differences observed between the predictions 
of the analytical model and the experimental results at Sensors 2 E,W and 3E,W, it is apparent that 
some adjustment of the analytical model is warranted. The modifications made to the model and 
the constants are presented and discussed in the next section. 

8.4 Internal Effects Model Modifications 
8.4.1 Model Modifications 

As noted in the preceding section, the model performs reasonably well in predicting the 
pressures experienced by Sensors 1E,W, but at Sensors 2 E,W and 3E,W, the predicted values are 
significantly lower than the experimental values. To lessen the rapidity with which the model 
predicts the attenuation of the peak pressures and temperatures, it is proposed that equation 
(152a), which calculate distance using a non-dimensional position quantity, be modified by 
introducing an empirical constant to control the rate of attenuation as follows: 

R = K 4 R’(E/P oa ) 1/3 (168) 


where K 4 > 1 . Choosing a value of K 4>1 would serve to “stretch out” the attenuation process so 

that at points far from the shock wave, the pressures are not quite as low as they would be 
without the empirical factor. 

The following expression for K 4 appears to capture the desired attenuation retardation 


feature: 


I 1 m 

[l + m(100R’) m 


,Xpoi ^ H 
5 x poi > H 


(169) 
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where m is a positive rational number, R’ is the non-dimensional attenuation distance parameter, 
and H is, as calculated in Section 6 . 2 . 2 , the distance away from the pressure wall at which internal 
debris cloud motion stops and attenuation of peak pressures and temperatures begins. This form 
is motivated by the desire to include “stretching effects” in the attenuation of peak pressures and 
temperatures only for points of interest that are relatively far removed from the shock wave or the 
internal debris cloud. Points of interest relatively near the shock wave or debris cloud experience 
“normally” attenuated pressure and temperature magnitudes. In this manner, some late time 
effects are included in the analytical model. Additionally, the modification proposed in equations 
(168) and (169) is applied only to the attenuation of debris cloud induced pressures and 
temperatures; the attenuation of shock wave induced quantities proceeds without any adjustment 
or stretching out of the attenuation. It is also noted that when m=0 (i.e. no “stretching out”), 
equation (168) does in fact reduce to its original form given in equation (152a). 

8.4.2 Empirical Constant Adjustments 

Following the modification of the analytical model as described in the preceding section, a 
value of m was determined that would extend the effects of the shock wave induced pressure and 
temperature increases further away from the shock wave and the debris cloud. In doing so, the 
empirical constants n, K 2 , and K 3 as well as the drag coefficient C[> were kept at the values shown 
in Table 7 The final values of the constant m are presented in Table 12 . 


Table 12 . Final Empirical Constant Values 


Empirical 

Constant 

Type A 
Inner Bumper 

TypeB 
Inner Bumper 

Projectile 

Diameter 

(cm) 

1.27 

1.59 

1.59 

m 

0.60 

0.70 

1.25 
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It must be emphasized that the constants in Tables 7 and 12 will likely yield acceptable 
results only for the system geometries, projectile diameters, and impact velocities considered in 
the experimental phase of this effort, the results of which served to determine the actual constant 
values. While the analytical model will certainly produce results for other system geometries, 
projectile diameters, and impact velocities, those results may or may not be appropriate for the 
alternative parameters considered. Additional experimental evidence will be required to validate 
the analytical model developed herein for impact parameters and systems configurations other 
than those which produced the empirical data upon which the empirical constants in Tables 7 and 
12 were based. It is indeed likely that these will need to be adjusted to accommodate other 
empirical information as it becomes available. 

8.5 Comparison of Modified Model Predictions with Experimental Results 

Tables 13 through 16 present a comparison of the experimental results and the predictions 
of the model developed herein using the empirical constants given in Tables 7 and 12. As can be 
seen in Tables 13 through 16, the predictions of the model now compare much more favorably 
with the experimental results. 

Table 13. Pressure and Temperature Data Comparison (Improved Model) — Test No. 1 



Pressure Increase 

Temperature Increase 1 

Sensor 

Time 

(msec) 

Experimental 
Result (psi) 

Model 

Prediction (psi) 

Time 

(msec) 

Experimental 
Result (°C) 

Model 

Prediction (°C) 

IE 

1.12 

17 

15.2 

13.3 

18 

43.2 

1W 

1.07 

18 

15.2 

19.0 

27 

43.2 

2E 

1.93 

4 

7.2 

21.6 

23 

23.0 

2W 

2.15 

6 i 

7.2 

23.0 

20 

23.0 

3E 

2.33 

3 

3.2 

645.8 

25 

11.6 

3W 

2.00 

2 

3.2 

n/a 

n/a 

11.6 
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Table 14. Pressure and Temperature Data Comparison (Improved Model) — Test No. 3 




Pressure Increase 

\ Temperature Increase 1 

Sensor 

Time 

(msec) 

Experimental 
Result (psi) 

Model 

Prediction (psi) 

Time 

(msec) 

sfTTr 


IE 

0.86 

21 

24.3 

8.5 

104 

63.9 

1W 

1.00 

21 

24.3 

6.5 

118 

63.9 

2E 

1.33 

17 

17.2 

10.9 

139 

46.6 

2W 

1.50 

21 

17.2 

10.8 

94 j 

46.6 

3E 

2.13 

p 8 ; 

8.0 

n/a 

n/a 

25.1 

3W 

2.83 

8 

8.0 

115,0 

40 

39.0 


Table 15. Pressure and Temperature Data Comparison (Improved Model) - Test No. 5 




Pressure Increase 

Temperature Increase 1 

Sensor 

Time 

(msec) 

Experimental 
Result (psi) 

Model 

Prediction (psi) 

Time 

(msec) 

Experimental 
Result (°C) 

Model 

Prediction (°C) 

IE 

0.95 

34 

34.2 

13.4 

64 

82.7 

1W 

0.80 

40 

34.2 

17.9 

54 

I 82.7 

2E 

2.13 

10 

10.2 

19.6 

34 

30.0 

2W 

1.93 

is 

10.2 

20.0 

31 

30.0 

3E 

4.07 

3 

5.6 

n/a 

n/a 

18.2 

3W 

3.90 

4 

5.6 

n/a 

n/a 

21.5 


Table 16. Pressure and Temperature Data Comparison (Improved Model) — Test No. 6 




Pressure Increase 

! Temperature Increase 1 

Sensor 

Time (msec) 

Experimental 
Result (psi) 

Model 

Prediction (psi) 

Time (msec) 

Experimental 
Result (°C) 

Model 

Prediction (°C) 

IE 

1.33 

16 

24.6 

25.0 

20 

64.8 

1W 

1.06 

31 

24.6 

26.3 

65 

64.8 

2E 

1.71 

15 

17.8 

92.0 

52 

48.1 

2W 

1.54 

14 

17.8 

53.0 

63 

48.1 

3E 

2.65 

12 

8.3 

52.0 

44 

25.8 

3W 

2.66 

7 

8.3 

n/a 

n/a 

40.1 


The most significant difference between experimental result and model prediction 
remaining is found for the temperature increases for Test No. 3 as shown in Table 14. 
Interestingly enough, the corresponding pressure increase values predicted by the model for this 
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test do agree rather well with the experimental results. A comparison of the test conditions, etc. 
for Test Nos. 3 and 6 reveals that the two are virtually identical; the only difference is that the 
velocity for Test No. 6 was 0.8 km/sec higher than that for Test No. 3. Inspection of the 
experimental results for these two tests as shown in Tables 14 and 16, respectively, shows that the 
pressure increases sustained during the two tests are in fact very similar. However, with the 
exception of Sensors 3E,W, the temperature increases are quite different: those recorded in Test 
No. 6 are substantially lower than those in Test No. 3. Since the analytical model predictions 
agree well with the pressure readings in both tests as well as the temperature readings in Test No. 
6, the inference to be drawn is that the temperature readings in Test No. 3 must be due to some 
phenomenon not included in the analytical model developed herein. 



9.0 SUMMARY AND RECOMMENDATIONS 

9.1 Summary 

A first-principles-based model has been developed to predict the effects of a perforating 
orbital debris particle impact within a pressurized habitable module. The model was developed in 
a way such that it sequentially characterizes the phenomena associated with the impact event, 
including the initial impact on the outer bumper, the creation and motion of the debris clouds 
within the module wall system, the impact of the debris clouds on the module pressure wall, the 
creation and motion of the debris cloud that enters the module interior, and the effects of the 
debris cloud within the module on module pressure and temperature levels. 

The model development has been accomplished through the application of elementary 
shock physics and thermodynamic theory, as well as the principles of mass, momentum, and 
energy conservation. The predictions of the model were compared against experimental results 
obtained from four high-speed impact tests conducted at the University of Alabama in Huntsville. 
The basic model predictions compared favorably with the experimental data obtained at the 
sensors located nearest to the impact site; at locations further removed from the impact site, the 
comparisons were not as favorable even taking into account test-to-test variations in the data. 
Several modifications were made to the basic model to bring its predictions closer in line with the 
experimental results, including decreasing the rapidity with which the model predicts the 
attenuation of peak pressures and accompanying temperature increases. Following the adjustment 
of several empirical constants, the predictions of the modified internal effects model were in much 
closer agreement with the experimental results. 
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9.2 Recommendations 


Following a review of the methodology used to develop the internal effects model 
presented herein, the following recommendations are offered as suggestions for improving the 
robustness of the model as well as for improving its ability to model the phenomena associated 
with the high speed impact of a pressurized module. The recommendations are grouped 
according to the part of the model that would be affected by their implementation. 

Outer Bumper Impact 

1) Include the effects of backsplash due to the impact of the projectile on the outer 
bumper . This would have the effect of decreasing the mass of the primary debris cloud, 
but, due to momentum conservation, increasing its center-of-mass velocity ( the so- 
called “momentum enhancement effect”). Since the energy of the debris cloud is 
proportional to the square of the velocity, this could have a significant impact on the 
subsequent effects produced by the primary debris cloud. 

2) Include the effects of light flash in the energy balance for the initial impact on the outer 
bumper . This would have the effect of producing more accurate primary debris cloud 
velocities. 

Inner Bumper Impact 

1) Develop a more suitable equation for the diameter of the hole in the inner bumper 
produced by the impact of the primary debris cloud . The present approximation is at 
best an order-of-magnitude estimate. A more appropriate hole-out equation would 
produce a more accurate mass value for this component of the secondary debris cloud. 
This in turn would have some effect on the magnitude of the velocity imparted to the 
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pressure wall by the impact of this component of the secondary debris cloud. 

2) Include the effects of inner bumper burning and/or melting . This may be a significant 
energy absorbing mechanism and may have a significant influence on the energy balance 
that is applied to the system before and after the impact of the primary debris cloud on 
the inner bumper. 

Pressure Wall Impact 

1) Include the effects of pressure wall deformation, including hole formation, bulging. 
cracking, etc. This will affect the amount of energy that is available for the motion of 
the internal debris cloud and the creation of the accompanying shock waves within the 
module interior. 

2) Develop an improved mass model for the internal debris cloud . This includes 
developing a more rigorous means of calculating pressure wall hole-out mass as well as 
an improved estimate for the amount of secondary debris cloud mass that is trapped by 
the pressure wall and the amount that enters the module interior. 

Internal Effects Modeling 

1) Improve the shock wave attenuation scheme for the internal debris cloud . This 
includes adjusting the present model to include an internal debris cloud whose speed 
decays as it moves into the module interior as opposed to ceasing motion in a rather 
abrupt manner. 

General Comments 

1) Extend t he model to oblique impacts . This would render the model applicable to a 
much wider variety of impact scenarios. 
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2) Perform additional internal effects testing at different velocities using alternative system 
geometries and impact conditions . This would provide additional information that 
would facilitate the implementation of the recommendations presented herein as well as 
validate the model for a much larger family of impact events. 

In terms of increasing the applicability of the model to a wider class of impact scenarios, 
additional internal effects testing at different velocities using alternative system geometries and 
impact conditions is absolutely required. This new information could then be used to, for 
example, extend the applicability of the model to oblique impacts. In terms of increasing the 
accuracy of the model developed herein, it would appear that, in order of decreasing importance, 
the improvements that need to be made are: 

• include the effects of pressure wall deformation, including hole formation, 
bulging, cracking, etc.; 

• develop an improved mass model for the internal debris cloud; 

• develop a more suitable equation for the diameter of the hole in the inner 
bumper produced by the impact of the primary debris cloud; 

• include the effects of inner bumper burning and/or melting; 

• improve the shock wave attenuation scheme for the internal debris cloud; 

• include the effects of backsplash due to the impact of the projectile on the 
outer bumper; and, 

• include the effects of light flash in the energy balance for the initial impact on 
the outer bumper. 
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APPENDIX A - INTEFF.FOR 



on no no 


$ DEBUG 


C 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM EFFCALC 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DOUBLE PRECISION KSB, KSP, KST, KA, KB, KP, KDC, KT, NUB, NUP, NUT, MB, MIB 
DOUBLE PRECISION MP, MW, MW1 , MW2 , MWH, MDC, MDC1,MDCT, PRT (10, 2 ) 
CHARACTER* 2 BID, PID,TID, IB, BIDCHK, PIDCHK, TIDCHK 
CHARACTER* 10 BMAT , PMAT , TMAT , I BMAT 
INTEGER VOPT 

COMMON/HYPERB/A, B, DEL, H 

OPEN ( 1 , FILE- * IMPDAT 1 ) 

OPEN ( 2 , FILE— * IMPOUT * ) 

OPEN ( 4 , FILE- 1 GPARAM 1 ) 

PI-3.141592 

READ DEBRIS CLOUD, PRESSURE WALL, AND AIR MATERIAL PROPERTIES. 
THE PARAMETERS MUST BE IN THE FOLLOWING UNITS: 


BID,PID,TID,IB MATERIAL ID CODES 

BMAT , PMAT , TMAT , I BMAT . . . MATERIALS 

COB, COP, COT BULK SOUND SPEED, KM/S 

RB, RP, RT AMBIENT MATL DENSITY, GM/CUCM 

ARIB INNER BMPR AREAL DNSTY, GM/SQCM 

KB, KP, KT SLOPE OF US-UP LINE 

EB,EP,ET ELASTIC MODULUS, LBS/SQ.IN. 

ALFAB , ALFAP , ALFAT LINEAR COEFF OF TERMAL EXP, 1/C 

CPSB, CPSP, CPST SPECIFIC HEAT (SOLID), CAL/GM-C 

CPLB, CPLP, CPLT SPECIFIC HEAT (LIQD) , CAL/GM/C 

TMB , TMP , TMT MELT TEMPERATURE, C 

TVB , TVP , TVT VAPORIZATION TEMPERATURE, C 

HFB, HFP, HFT LATENT HEAT OF FUSION, CAL/GM 

HVB, HVP , HVT LATENT HEAT OF VPRZTN, CAL/GM 


WRITE ( * , 3 ) 

3 FORMAT ( 1 ENTER 
READ (4,5) PID 
5 FORMAT (A2) 

WRITE (*,31) 
31 FORMAT {' ENTER 
READ (4, 51) BID 
51 FORMAT (A2) 

WRITE (*, 7) 

7 FORMAT ( 1 ENTER 
READ (4,9) TID 
9 FORMAT (A2) 

WRITE (*, 11) 
11 FORMAT (' ENTER 
READ (4,111) IB 
111 FORMAT (A2) 

IF (IB. EQ. 1 AA' ) 


PROJECTILE MATERIAL ID CODE (A2) AND HIT ENTER*) 

BUMPER MATERIAL ID CODE (A2 ) AND HIT ENTER*) 

PRESSURE WALL MATERIAL ID CODE (A2) AND HIT ENTER’) 

INNER BUMPER ID CODE (A2) AND HIT ENTER') 

THEN 


I BMAT- * MLI-BLNKT ' 


RIBA-0.072 

ENDIF 

IF (IB. EQ. ' BB ' ) THEN 
IBMAT- 'MLI/K/N/GE* 
RIBA-1.02 


ENDIF 


C 
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REWIND 1 
READ (1, 4) 

4 FORMAT (////) 

C 

99 READ (1,1) PIDCHK 

1 FORMAT (A2) 

IF (PID.EQ. PIDCHK) THEN 
READ (1,10) PMAT, COP, KP, RP, GPI 
10 FORMAT (A10,4F10. 5) 

READ (1,100) EP, NUP, ALPHAP, CPSP, CPLP 
100 FORMAT (2( El 0.3, FI 0.5), FI 0.5) 

READ (1,102) TMP , TVP , HFP , HVP 
102 FORMAT (4F10. 5) 

ENDIF 

IF (PID.NE. PIDCHK) THEN 
IF (PIDCHK. EQ. 'XX' ) THEN 
WRITE (*,17) 

17 FORMAT (' PROJECTILE MATERIAL NOT FOUND IN MATERIAL LIBRARY.',/, 

$’ PLEASE CHECK DEBRIS CLOUD MATERIAL ID CODE AND BEGIN AGAIN.') 
STOP 
ENDIF 

IF (PIDCHK. NE. 'XX' ) THEN 
READ (1,2) 

2 FORMAT (///) 

GOTO 99 

ENDIF 

ENDIF 

C 

REWIND 1 
READ (1, 4) 

999 READ (1,1) BIDCHK 

IF (BID. EQ. BIDCHK) THEN 

READ (1,10) BMAT, COB, KB, RB, GBI 

READ (1, 100 ) EB, NUB, ALPHAB, CPSB, CPLB 

READ (1,102) TMB, TVB, HFB, HVB 

ENDIF 

IF (BID. NE. BIDCHK) THEN 
IF (BIDCHK. EQ. 'XX' ) THEN 
WRITE (*,117) 

117 FORMAT (' BUMPER MATERIAL NOT FOUND IN MATERIAL LIBRARY.’,/,' PLEAS 
$E CHECK BUMPER MATERIAL ID CODE AND BEGIN AGAIN.') 

STOP 

ENDIF 

IF (BIDCHK. NE. 'XX' ) THEN 
READ (1,2) 

GOTO 999 

ENDIF 

ENDIF 

REWIND 1 
READ(1, 4) 

9999 READ (1,1) TIDCHK 

IF (TID.EQ. TIDCHK) THEN 

READ (1,10) TMAT, COT, KT, RT, GTI 

READ (1,100) ET, NUT, ALPHAT, CPST, CPLT 

READ (1,102) TMT, TVT, HFT, HVT 

ENDIF 

IF (TID.NE. TIDCHK) THEN 
IF (TIDCHK. EQ. 'XX' ) THEN 
WRITE (*,1117) 
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non 


1117 FORMAT (' PRESSURE WALL MATERIAL NOT FOUND IN MATERIAL LIBRARY.',/ 
PLEASE CHECK PRESSURE WALL MATERIAL ID CODE AND BEGIN AGAIN.') 
STOP 
ENDIF 

IF (TIDCHK.NE. 'XX' ) THEN 
READ (1,2) 

GOTO 9999 

ENDIF 

ENDIF 


C 

C REQUIRED PROPERTIES OF AIR AND THEIR UNITS: 

C 

C COA ... ADIABATIC BULK SOUND SPEED, KM/S 

C CPA . . . AMBIENT SPECIFIC HEAT, CAL/GM/DEG-C 

C RP AMBIENT DENSITY, GM/ CU.CM. 

C KA SLOPE OF US-UP LINE, DIMENSIONLESS 

C TOA . . . AMBIENT TEMPERATURE, DEG-C 

C POA ... AMBEINT PRESSURE, ATM 

C CD . . . . DRAG COEFFICIENT 

C 


CPA=0 . 24 

COA=0 . 33 

KA=1. 07 

RA=1 . 22E-03 

T0A=24 . 0 

P0A=1 . 0 

TA=T0A+273 . 0 

PA=P0A*101325 . 0 

E0A= ( CPA* 418 6. 0 ) *TA 


C 

C READ GEOMETRIC PARAMETERS AND ASSUMPTIONS 

C 

C DP ... . PROJECTILE DIAMETER, CM 

C TS . . . . BUMPER THICKNESS, CM 

C TW PRESSURE WALL THICKNESS, CM 

C S BUMPER-TO-PRESSURE WALL STAND-OFF DISTANCE, CM 

C S2 INNER BUMPER-TO-PRESSURE WALL DISTANCE, CM 

C HMOD . . MODULE DIAMETER, CM 

C 


READ (4,113) DP, TS, TW, S, S2, HMOD, CD 
113 FORMAT (7F10. 5) 

C 

C READ PROJECTILE IMPACT VELOCITY IN KM/S 

C 

WRITE (*,29) 

29 FORMAT (' INPUT PROJECTILE IMPACT VELOCITY IN KM/SEC (F5.2) AND HIT 
$ ENTER') 

READ (*,30) VP 

30 FORMAT (F5. 2) 

CALCULATE PROJECTILE AND BUMPER HOLE-OUT MASSES 

MP= ( PI/6 . 0 ) * (DP/ 100. 0) * (DP/ 100.0) * (DP/100 . 0) * (RP*1000 . 0 ) 

R1=VP/ COB 
R2=TS/DP 
R3=RB/RP 

DBDP=3 . 4* (R2**0. 33333) * (R1**0 . 33333 ) * ( 1 . 0-0 . 0308*R3 ) 

DB=DP*DBDP 

MB= (PI/4.0)* (DB/100. 0) * (DB/100. 0) * (TS/ 100.0) * (RP*1000. 0) 

C 
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EB=EB*68947 . 0 
BETAB=3 . 0*ALPHAB 
IF (NUB. LT. 0.5) THEN 
KSB=EB/3 .0/(1. 0-2 . 0*NUB) 

COBC=DSQRT ( (KSB/10. 0) / (RB*1000. 0) ) / 1000.0 
ENDIF 

IF (NUB. EQ. 0.5) THEN 
KSB=-1 . 0 
COBC=-l . 0 
ENDIF 

IF (NUB. LT .0.5) GB=2 . 3885E-08*KSB*BETAB/CPSB/RB 
IF (NUB. EQ. 0.5) GB=GBI 
GRB=GB*RB*1000. 0 
C 

EP=EP*68947 . 0 
BETAP=3 . 0*ALPHAP 
IF (NUP.LT.0.5) THEN 
KSP=EP/ 3 . 0/ ( 1 . 0-2 . 0*NUP) 

COPC=DSQRT ( (KSP/10. 0) / (RP*1000 . 0) ) /1000. 0 
ENDIF 

IF (NUP.EQ.0.5) THEN 
KSP=-1. 0 
COPO-1. 0 
ENDIF 

IF (NUP.LT.0.5) GP=2 . 3885E-08*KSP*BETAP/CPSP/RP 
IF (NUP.EQ.0.5) GP=GPI 
GRP=GP*RP*1000. 0 
C 

WRITE (2,40) PMAT , BMAT , I BMAT , TMAT 

40 FORMAT (' NORMAL IMPACT OF A ' ,A10,' PROJECTILE ON A DUAL-WALL SYSTE 
$M WITH A' ,/,A10, ’ BUMPER, A ',A10,' INNER BUMPER, AND A ',A10,' PR 
SESSURE WALL’) 

WRITE (2,45) PMAT , COP, KP, RP, DP, MP* 1000 . 0, VP, BMAT, COT, KT, RT, TS, DB 
$ ,MB*1000 . 0 

45 FORMAT!/, 'PROJECTILE PROPERTIES ...',/, 3X, 'MAT = ',A10,/,3X, 

$’CO = 1 , F6 . 3, ' KM/S', /, 3X, 'K = ' , F6-3,/,3X, 'RHO = ',F6.3,' GM/C 
$U . CM. 1 , / , 3X, ' DP = 1 , F6 . 3, ' CM',/,3X,'MP = ',F6.3,' GMS',/,3X, 

$'VP = 1 , F6 . 3, ' KM/S' ,//, 'BUMPER PROPERTIES ...',/, 3X, 'MAT = ',A10 
$,/,3X,'CO = ' , F6. 3, ' KM/S' ,/,3X, 'K = ' , F6. 3, /, 3X, ' RHO = ',F6.3, 

$' GM/CU.CM. ',/,3X, 'TS = ' , F6 . 3, ' CM',/,3X,’DH = ',F6.3,' CM',/, 

$3X, 'MB = ' , F6. 3, ' GMS ' ) 


C 

C CALCULATE PARTICLE AND SHOCK WAVE VELOCITIES AND HUGONIOT 

C PRESSURE DUE TO PROJECTILE IMPACT 

C 


V=VP 

IF (BMAT. EQ. PMAT) GOTO 35 
A=KP-KB* (RB/RP) 

B=2 . 0*KP*V+COP+COB* (RB/RP) 

C=COP*V+KP*V*V 

D=B*B-4 . 0*A*C 

UBP= ( B-SQRT (D) ) / (2. 0*A) 

GOTO 38 

35 UBP=V/2.0 

38 UPP=V-UBP 
UBS=COB+KB*UBP 
UPS=COP+KP*UPP 
PP=RP*UPS*UPP 
PB=RB*UBS*UBP 
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c PROJECTILE AND BUMPER SHOCK LOADING RESPONSE AND RELEASE 

C CALCULATION PHASE 

C 

WRITE!*, 5080) 

5080 FORMAT!/,' BEGINNING PROJECTILE AND BUMPER SHOCK LOADING RESPONSE 
$AND ' , / , ' RELEASE CALCULATIONS . . . ' ) 

WRITE(2, 509) 

509 FORMAT!/,'**** PROJECTILE AND BUMPER SHOCK LOADING RESPONSE AND RE 
$ LEASE CALCULATIONS ****■) 

WRITE (2, 6011) VP,UPP,UPS, PP,UBP,UBS, PB 

6011 FORMAT (/,' PROJECTILE IMPACT VELOCITY VP = ',F7.3,' KM/S',/, 

$ ' PROJ MATL PARTICLE VELOCITY ... UP = ',F7.3,' KM/S ' , / , ' PROJ MATL 

$ SHOCK WAVE SPEED US = ',F7.3,' KM/S ’ , /, 'HUGONIOT IMPACT PRESS 

SURE PH = ' , F7 . 3 , ' GPA',/,'BMPR MATL PARTICLE VELOCITY ... U 

$P = ' , F7 . 3 , ' KM/S' ,/, 'BMPR MATL SHOCK WAVE SPEED US = ',F7.3, 

$' KM/S' ,/, 'HUGONIOT IMPACT PRESSURE PH = ',F7.3,' GPA' ) 

C 

VPO=l. 0/RP 
VP1=RP*UPS/ (UPS-UPP) 

VP1=1.0/VP1 

C 


PH=PP* 1 . 0E09 
C 

WRITE (2, 705) EP/10 . 0, NUP, KSP/10 . 0, ALPHAP, CPSP, CPLP 
705 FORMAT (/,' PARAMETERS REQUIRED FOR CALCULATING PROJECTILE MATERIAL 
$ RESPONSE AND ',/,' RELEASE FROM SHOCKED STATE USING THE MIE-GRUNEISE 

$N E-O-S: ',/,3X, 'ELASTIC MODULUS E =',E10.4,' N/SQ 

$.M.',/,3X, 'POISSON RATIO NU = ’ , F10 . 3, /, 3X, ' BULK 


$ MODULUS K = ' , E10 . 4, 1 N/SQ. M. ’,/, 3X, ' LIN. COEF. 

$ OF THERM. EXP. ... ALFA = ' , E10 . 4 , ’ /DEG-C' , /, 3X, ' SP HEAT (SOLID) 

$ CPS = ' , F10 . 3, ' CAL/ GM/ DEG-C ' , / , 3X , ' S P HEAT (LIQUID) 

$ CPL = ' , F10 . 3, ' CAL/GM/DEG-C’ ) 


PHMB=PH/100 . 0E+09 

WRITE (2, 800) PH, PHMB, VPO, VP1, GP, GPI 
800 FORMAT (3X, 'HUGON IMP PRESS ( PA, MBAR) ... PH =' , E10 . 4, ' , ' , F5 . 3, /, 

$3X, ' SP VOL AT REST V0 =',F10.3, ' CU . CM. /GM’ , / , 3X, 

$'SP VOL AT IMPACT VI =',F10.3,' CU.CM. /GM' ,/, 3X, 'AM 

$B M-GRUN COEF (CAL,INP) ... GAM0 = ' , F10 . 3, ' , ’ , F5 . 3) 

WRITE (2, 805) TMP , TVP , HFP , HVP 

805 FORMAT (3X, 'MELT TEMPERATURE TM =',F10.2,' DEG-C',/, 

$3X, 'VAPOR TEMPERATURE TV =',F10.2,’ DEG-C ',/, 3X, ' HEA 


$T OF FUSION HF =',F10.2, ' CAL/ GM' ,/, 3X, ' HEAT OF V 

$APORIZATION HV =',F10.2,' CAL/GM' ) 

C 


C CALCULATE RELEASE OF PROJECTILE MATERIAL UP UNTIL ZERO PRESSURE 

C IS REACHED 

C 


WRITE (2,1701) 

1701 FORMAT!/, 'RELEASE OF SHOCKED PROJECTILE MATERIAL ...’) 
PFIN=0 . 0 

CALL RELS ( COP , KP , RP , GRP , VPO , VP1 , PH , EXTP , UPP , PFIN , VFP ) 

C 

C CALCULATE TEMPERATURE INCREASE IN PROJECTILE MATERIAL 

C 

CALL TINC (CPSP, CPLP, TMP, TVP, HFP, HVP, EXTP) 

C 

VBO=l . 0/RB 
VB1=RB*UBS/ (UBS-UBP) 

VB1=1 . 0/VB1 
C 
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n n a o noon 


PH=PB*1. 0E+09 


WRITE (2,7051) EB/10 . 0, NUB, KSB/10 . 0, ALPHAB, CPSB, CPLB 
7051 FORMAT (/,' PARAMETERS REQUIRED FOR CALCULATING BUMPER MATERIAL RESP 
$ONSE AND' ,/, 'RELEASE FROM SHOCKED STATE USING THE MIE-GRUNEISEN E- 

$0-S : ' , / , 3X, 1 ELASTIC MODULUS E E10 . 4 , ' N/SQ . M ' 

$,/, 3X, ' POISSON RATIO NU =' , F10. 3, /, 3X, ' BULK MODU 

* LUS K = 1 , E10 . 4 , ' N/SQ. M. ' ,/,3X, ’LIN. COEF. OF 

$THERM. EXP. ... ALFA =' , E10. 4, ' /DEG-C ' , /, 3X, ’ SP HEAT (SOLID) 

f CPS = ' , F10 . 3, ' CAL/GM/DEG-C’ ,/,3X, 'SP HEAT (LIQUID) ... 

^ CPL = ' , F10 . 3 , 1 CAL/GM/DEG-C) 

PHMB=PH/100. 0E+09 

WRITE (2, 8001) PH, PHMB, VBO, VB1, GB, GBI 
8001 FORMAT (3X, ' HUGON IMP PRESS (PA, MBAR) ... PH =’,E10.4,' ' F5 . 3 / 

$3X, 'SP VOL AT REST V0 =',F10.3,' CU . CM . / GM ' , / ,3x', 

$ SP VOL AT IMPACT VI =',F10.3, ’ CU.CM. /GM' ,/, 3X, 'AM 

$B M-GRUN COEF (CAL,INP) ... GAM0 =' , F10. 3, ' , ' , F5.3) 

WRITE(2, 8051) TMB, TVB, HFB, HVB 

8051 FORMAT (3X, 'MELT TEMPERATURE TM =',F10.2,' DEG-C’ / 

$3X, 'VAPOR TEMPERATURE TV =',F10.2,’ DEG-C ',/, 3X, ' HEA 

$T OF FUSION HF =',F10.2,' CAL/GM' , /, 3X, ’ HEAT OF V 

SAPORIZATION HV =',F10.2,' CAL/GM') 


CALCULATE RELEASE OF BUMPER MATERIAL UP UNTIL ZERO PRESSURE IS 
REACHED 


WRITE (2,1703) 

1703 FORMAT (/, 'RELEASE OF SHOCKED BUMPER MATERIAL ') 
PFIN=0 .0 

CALL RELS (COB, KB, RB, GRB, VBO, VB1, PH, EXTB,UBP, PFIN, VFB) 


CALCULATE TEMPERATURE INCREASE IN BUMPER MATERIAL 
CALL TINC (CPSB, CPLB, TMB, TVB, HFB, HVB, EXTB) 


WRITE {*, 1707) 

1707 FORMAT { ' PROJECTILE AND BUMPER MATERIALS SHOCK LOADING AND RELEASE 
$ 1 , / , ' CALCULATIONS COMPLETE . ' ) 


WRITE(*, 58) 

58 FORMAT!/, ' BEGINNING PRESSURE WALL RESPONSE AND IMPEDANCE MATCH CA 
$LCULATIONS . . . ' ) 


CALCULATE PRIMARY DEBRIS CLOUD MASS AXIAL VELOCITY, AND EXP VEL 


C 

c 


MDC=MP+MB 

RATIOP=MP/MDC 

RATIOB=MB/MDC 

VCOM=MP*VP/ (MP+MB) 


EINIT=0 . 5*MP* (VP* 10 00.0) * (VP* 100 0.0) 

EDIFF=EINIT-EXTB*MB-EXTP*MP-0 . 5*MDC* (VCOM*1000. 0) * (VCOM*1000 0) 

IF { EDIFF . LE .0.0) THEN 
WRITE (2,991) 

WRITE (*,991) 

991 FORMAT ( 1 * * * FATAL ERROR 001: INSUFFICIENT ENERGY AVAILABLE FOR PRI 
$MARY DEB CLD EXPANSION *++») 

STOP 

ENDIF 


101 



c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


VEXP=DSQRT ( EDIFF/ { 0 . 5*MDC) ) / 1000 . 0 

CALCULATE PRIMARY DEBRIS CLOUD LEADING EDGE VELOCITY 

VLE=2 . 0*UBP 

CALCULATE PRIMARY DEBRIS CLOUD 1/2-ANGLE SPREAD, PRESSURE WALL 

LOADING AREA, MASS OF SHOCKED AND RELEASED PRESSURE WALL MATL, 

AND MASS OF PRESSURE WALL HOLE-OUT MATERIAL 

THDC=ATAN (VEXP/VCOM) 

THDEG=180 . 0*THDC/PI 

DW=2 . 0 * S * DTAN ( THDC ) 

RW=DW/2 . 0 

MW=< PI/4.0) * (DW/100.0) * (DW/100.0) * (TW/100.0) * (RT*1000.0) 


READ (4,49) AK2 , AK3 , AN , AMEXP 
49 FORMAT ( 4 FI 0 . 5 ) 

C 


c CALCULATE DEBRIS CLOUD IMPACT VELOCITY BASED ON MOMENTUM 

C CONSERVATION AT THE PRESSURE WALL 

C 

VPW=MDC*VLE/ (MDC+MW) 

C 

c ASSIGN DEBRIS CLOUD IMPACT VELOCITY 

C 

READ (4,1709) VOPT 


1709 FORMAT (II) 

C 

IF (VOPT. EQ. 1) V=VCOM 
IF (VOPT.EQ.2) V=VLE 
IF (VOPT.EQ. 3) V=VPW 
C 

WRITE (2,1710) VEXP, VCOM, VLE, VPW, V 

1710 FORMAT (/,' PRIMARY DEBRIS CLOUD CHARACTERISTIC VELOCITIES ...',/, 3X 
$, 'VEXP = 1 , F6 . 3, ' KM/S',/,3X, 'VCOM = ',F6.3, ' KM/S ' , /, 3X, ' VLE = ' 
$, F6. 3, ' KM/S ',/, 3X, 'VPW = ',F6.3,' KM/ S ',/,' SELECTED PRIMARY DEBR 
$IS CLOUD IMPACT VELOCITY ...’,/, 3X, 'VDC = ',F6.3,' KM/S') 

C 

ET=ET*68947 . 0 
BETAT=3 . 0*ALPHAT 
IF (NUT.LT.0.5) THEN 
KST=ET/3.0/ (1.0-2. 0*NUT) 

COTC=DSQRT( (KST/10 . 0) / (RT*1000 . 0) )/1000.0 
ENDIF 

IF (NUT. EQ. 0.5) THEN 
KST=- 1 . 0 
COTC=-l . 0 
ENDIF 

IF (NUT.LT.0.5) GT=2 . 3885E-08*KST*BETAT/CPST/RT 
IF (NUT. EQ. 0.5) GT=GTI 
GRT=GT*RT*1000 . 0 
C 

KDC=KP*RATIOP+KB*RATIOB 
RPSR=1 . 0/VFP 
RBSR=1 . 0/VFB 

RDC=RP S R* RAT I O P + RBS R* RATI OB 
CODC=CO P * RAT I O P +COB * RAT I OB 
GRDC=GRP*RATIOP+GRB*RATIOB 
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c 

WRITE (2,50) THDEG, RATIOP, RPSR, RATIOB, RBSR, CODC, KDC, RDC, MDC*1000 . 
50 FORMAT (/,' PRIMARY DEBRIS CLOUD PROPERTIES 3X, ' 1/2-ANG = 

$F6 . 3, ' DEG' , / , 3X, ' RATIO-P = ' , F6. 3, /, 3X, ' RPSR = \F6.3, 

$’ GM/CU. CM. ' , 3X, ' RATIO-B = ' , F6 . 3, /, 3X, ' RBSR = ',F6.3, 

$' GM/CU. CM. ',/,3X, 'CO = 1 , F6. 3, ' KM/S (AVG) ' , 3X, ' K 

$ ' , F6 . 3, ' (AVG) ',/,3X, ' RHO = ’,F6.3,' GM/CU. CM. (AVG)',/,3X, 

$ 'MDC = 1 , F6 . 3, 1 GMS ' ) 


C 

C CALCULATE SECONDARY DEBRIS CLOUD COMPONENT MASSES, AXIAL MOTION 

C VELOCITIES, AND EXPANSION VELOCITIES 

C 


DH1= ( 0 . 45*V+0 . 90) *DB 

MIB=PI* (RIBA* 10. 0) * (DH1/100 . 0 ) * (DH1/100. 0) /4. 0 
RWA=RT*TW 

AK1= ( 1 . 0-RIBA/RWA) **AN 
GAM=MIB/MDC 

AKMN= ( -GAM+DSQRT ( GAM*GAM+4 . 0 ) )/2.0 
IF (AKMN.GT.AK1) AK1=AKMN 
MDC1=AK1*MDC 
MDCT=MDC 1 +MI B 
C 

VMCA=MDC*V/MDCT 

C 

R2=V/VMCA 

CALL ROOTS (AK1, R2,MIB,MDC,MDC1, C1AX, C2AX) 

VMC1=C1AX*VMCA 

VMC2=C2AX*VMCA 

C 

EPA=MDC* (V*1000.0) * (V*1000 . 0 ) /2 . 0 
EPE=MDC* (VEXP*1000. 0) * (VEXP*1000. 0) /2.0 
EL1=MDCT* (VMCA*1000 . 0) * (VMCA*1000. 0) /2 . 0 
EL21=MDC* ( 1 . 0-AK1 )*(V* 1000. 0)*(V* 1000. 0)/2.0 
EL22=MDC* (1.0-AK1)* (VEXP* 1000.0) * (VEXP*1000.0)/2.0 
EDIF=EPA+EPE-EL1- (EL21+EL22) 

IF (EDIF.LE.0.0) THEN 
WRITE (2,9921) 

WRITE (*,9921) 

9921 FORMAT ('*** FATAL ERROR 002: INSUFFICIENT ENERGY AVAILABLE FOR SEC 
$ONDARY DEBRIS CLOUD EXPANSION *** ' ) 

STOP 

ENDIF 

VEXPA=DSQRT (EDIF/ (MDCT/2.0) ) / 1000.0 
C 

R2=VEXP/VEXPA 

CALL ROOTS (AK1, R2,MIB,MDC,MDC1, C1EXP, C2EXP) 

VEXP1=C1EXP*VEXPA 
VEX P2 =C2 EX P * VEX P A 
C 

c CALCULATE SECONDARY DEBRIS CLOUD SPREAD ANGLES 

C 

THDC1=ATAN (VEXP1/VMC1) 

TH1DEG=180 . 0*THDC1/PI 
THDC2=ATAN (VEXP2/VMC2 ) 

TH2 DEG= 180 . 0*THDC2/PI 
DW1=2.0*S2*THDC1 
DW2=2 . 0*S2*THDC2 
RWl=DWl/2 . 0 
RW2=DW2/2 . 0 

MW1= (PI/ 4.0)* (DW1/100. 0) * (DW1/100 . 0) * (TW/100. 0) * (RT*1000 . 0) 
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non 


MW2= (PI/4 . 0) * (DW2/100. 0) * (DW2/100. 0) * (TW/100. 0) * (RT*1000. 0) 

.... CALCULATE ADDITIONAL PRESSURE WALL VELOCITY 

MDC=MDC1 
MW=MW1 

MWH=AK3 * (MW1*AK2 ) 

V=VMC1 
VEXP=VEXP1 

UPW=MIB* (V+VEXP) / (MIB+MW2) 

WRITE (2,450) RIBA, DH1,S2 

450 FORMAT (/, 'INNER BUMPER PROPERTIES 3X , ’ RHO-A = ’,F6.3, 

$’ GM/SQ.CM. ’ ,/,3X, 'DH = \F6.3,' CM',/,3X,'S2 = \F6.3, 

$' CM’) 

WRITE (2,451) VMCA, EPA, EPE , ELI , EL2 1 , EL22 , EDI F, VEXPA, UPW 

451 FORMAT (/,' SECONDARY DEBRIS CLOUD PROPERTIES (GENERAL) 3X, 

$’SEC DEB CLD AVG AXL VEL ',F6.3,' KM/S ' , /, 3X, ' PRIM DEB CLD AX 

$L ENERGY ',E11.4,' J' , /, 3X, ' PRIM DEB CLD EXP ENERGY ', 

$E11 . 4 , ' J' ,/,3X, 'SEC DEB CLD AXL ENERGY ',E11.4,' J',/,3X, 

$ ' PRIM DEB CLD AXL EN LOST ... 1 ,E11.4,' J' , /, 3X, ' PRIM DEB CLD EXP 
$EN LOST ... * , Ell . 4 , ' J ',/, 3X, ' EN AVAIL SEC DEB CLD EXP ... ', 

$E11. 4, ' J’,/,3X, 'SEC DEB CLD AVG EXP VEL ' ,F6.3, ' KM/S',/,3X, 

$ ' ADDL PRESS WALL VEL ',F6.3,' KM/S') 

WRITE (2, 452) TH1DEG, RDC, AN, AK1, AKMN, MDC* 1000 .0, C1AX, C1EXP, V, VEXP 

452 FORMAT (/, 'SECONDARY DEBRIS CLOUD PROPERTIES (1ST COMP) ...', 

$/,3X, '1/2-ANG = ' , F6 . 3, ' DEG’ , /, 3X, ’ RHO = ',F6.3,' GM/CU.CM.', 
$/,3X,’N = 1 , F6. 3, /, 3X, ' K1 = ’,F6.3,' (>', F5 . 3, ’)',/, 3X, 

$ 'MDC1 = ' , F6. 3, ' GMS ' , / , 3X, ' C#1 = ' , F6. 3, / , 3X, ' C#2 = ', 

$F6. 3, /, 3X, 'VMC = ' , F6. 3, ’ KM/S ’,/, 3X, 'VEXP = ',F6.3,' KM/S' 

$) 

WRITE (2, 453) TH2DEG, RIBA, MIB* 1000 . 0, C2AX, C2EXP, VMC2, VEXP2 

453 FORMAT (/,' SECONDARY DEBRIS CLOUD PROPERTIES (2ND COMP) ...',/, 3X, 

$' 1/2-ANG = ' , F6 . 3, ' DEG' ,/,3X, 'RHO-A = ',F6.3,' GM/CU . CM. ' , /, 3X, 

$ 'MIB = 1 , F6. 3, ' GMS ' , /, 3X, ' C#1 = • , F6. 3, /, 3X, ' C#2 = ', 

$F6.3,/,3X, 'VMC = ' , F6. 3, ' KM/S ',/, 3X, 'VEXP = *,F6.3, 

$• KM/S') 

WRITE (2,454) TMAT, COT, KT, RT,TW, S, COA, KA, RA, CPA,T0A, P0A, HMOD 

454 FORMAT (/, 'PRESSURE WALL PROPERTIES ...',/, 3X, 'MAT = ' ,A10,/,3X, 
$'CO = 1 , F6 . 3 , ' KM/S' ,/, 3X, 'K = ' , F6 . 3, /, 3X, ' RHO = ',F6.3, 

$' GM/CU.CM. /,3X, 'TW = ',F6.3,' CM',/,3X,'S = ',F6.3,' CM’,//, 

$ 'MODULE AIR PROPERTIES ...',/, 3X, ' CO = ',F6.3,' KM/S', /,3X, 

$'K = ' , F6. 3, /, 3X, ' RHO = ’,E10.3,' GM/CU.CM. ',/, 3X, ’ CP = ',F6.3, 

$' CAL/GM/DEG-C' ,/,3X, 'TO = ',F6.3,' DEG-C* , /, 3X, ' P0 = ',F6.3, 

$' ATM',//, 'MODULE DIAMETER ...’,/, 3X, ’ HMOD = ',F7.2,' CM') 

C CALCULATE PARTICLE AND SHOCK WAVE VELOCITIES AND HUGONIOT 

C PRESSURE DUE TO SECONDARY DEBRIS CLOUD IMPACT 

C 

IF (TMAT. EQ. BMAT. AND. BMAT. EQ. PMAT) GOTO 351 
A=KDC-KT* ( RT/RDC) 

B=2.0*KDC*V+CODC+COT* (RT/RDC) 

C=CODC*V+KDC*V*V 

D=B*B-4.0*A*C 

UTP= ( B-SQRT ( D) )/(2.0*A) 

GOTO 381 
351 UTP=V/2.0 

381 UDCP=V-UTP 
UTS=COT+KT*UTP 
UDCS=CODC+KDC*UDCP 
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PDC=RDC*UDCS*UDCP 

PT=RT*UTS*UTP 

C 

c PRESSURE WALL RESPONSE AND IMPEDANCE MATCH CALCULATION PHASE 

C 

WRITE (2,59) 

59 FORMAT (/,' **** PRESSURE WALL RESPONSE AND IMPEDANCE MATCH CALCULAT 
$IONS ****') 

WRITE (2, 60) V, UTP, UTS, PT 

60 FORMAT (/, 'ASSIGNED PRESS WALL IMP VEL VW = ',F7.3,' KM/S', 

$/,' PRESSURE WALL PARTICLE VELOCITY ... UP = ’ ,F7.3,' KM/S ',/,' PRES 
$SURE WALL SHOCK WAVE SPEED .... US = ',F7.3,' KM/S ' , / , ’ HUGONIOT IM 
$PACT PRESSURE PH = ',F7.3,' GPA' ) 

C 

VTO=l . 0/RT 

VT1=RT*UTS/ (UTS-UTP) 

VT1=1 . 0/VT1 

VDCO=l . 0/RDC 

VDC1=RDC*UDCS/ (UDCS-UDCP) 

VDC1=1 . 0/VDC1 
C 

PH=PT*1 . 0E09 
C 

WRITE (2, 75) ET / 1 0 . 0 , NUT , KST/10 . 0, ALPHAT, CPST, CPLT 

75 FORMAT (/,' PARAMETERS REQUIRED FOR CALCULATING PRESSURE WALL MATERI 
$AL RESPONSE AND ',/,’ RELEASE FROM SHOCKED STATE USING THE MIE-GRUNE 

$ISEN E-O-S : ' , 3X, 'MATL ELASTIC MODULUS E =', 

$E10 . 4 , 1 N/SQ.M. ' ,/,3X, 'MATL POISSON RATIO NU =', 

$F10. 3, /,3X, 'MATL BULK MODULUS K =’,E10.4, 

$' N/SQ.M. ',/, 3X, 'MATL LIN. COEF. OF THERM. EXP. ... ALFA = ' , E10 . 4 , 

$' /DEG-C' ,/,3X, 'MATL SP HEAT (SOLID) CPS =',F10.3, 

$’ CAL/GM/DEG-C' , /, 3X, 'MATL SP HEAT (LIQUID) CPL =' 

$, F10 . 3, ' CAL/GM/DEG-C’) 

PHMB=PH/100. 0E+09 

WRITE (2,80) PH, PHMB, VTO, VT1, GT, GTI 

80 FORMAT (3X, 'MATL HUGON IMP PRESS ( PA, MBAR) ... PH = ’ , E10 . 4, ' , ' , 

$F5 . 3, / , 3X, 'MATL SP VOL AT REST V0 =',F10.3, ' CU.C 

$M./GM' , /, 3X, 'MATL SP VOL AT IMPACT VI =',F10.3, ' CU 

$.CM./GM' ,/, 3X, 'MATL AMB M-GRUN COEF (CAL,INP) ... GAM0 =',F10.3,', 
$ ' , F5 . 3) 

WRITE (2,85) TMT, TVT, HFT, HVT 

85 FORMAT ( 3X ,' MATL MELT TEMPERATURE TM =',F10.2,' DEG- 

$C*,/.3X, 'MATL VAPOR TEMPERATURE TV =',F10.2,' DEG-C' 

$,/,3X, 'MATL HEAT OF FUSION HF =',F10.2,' CAL/GM' 

$,/,3X, 'MATL HEAT OF VAPORIZATION HV =',F10.2,' CAL/GM') 


C 

C PERFORM IMPEDANCE MATCH CALCULATIONS TO DETERMINE SHOCK WAVE 

C PROPERTIES FOR WAVE TRANSMITTED INTO MODULE AND REFLECTED INTO 

C PRESSURE WALL 

C 


CALL IMPMTCH ( RT, RA, KT, KA, COT, COA, UTP, PA, UAP) 

UAS=COA+KA*UAP 
PHA=RA*UAS*UAP 
PHAB=(PHA/ 101325.0) ‘1.0E+09 
C 

WRITE (2,175) UAP, UAS , PHA, PHAB 

175 FORMAT (/, 'RESULTS OF IMPEDANCE MATCH CALCULATIONS ...',/, 3X, 

$ ' INTERFACE (REFL i TRANS) PARTICLE VEL ...',F8.5,' KM/S', /,3X, 


$' TRANSMITTED SHOCK WAVE SPEED ',F8.5,' KM/S',/,3X, 

$’ INTERFACE (REFL & TRANS) PRESSURE ' ,E12.5,’ GPA (',F5.1, 
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$ 


ATM) 1 ) 


C 

c CALCULATE RELEASE OF SECONDARY DEBRIS CLOUD MATERIAL UP UNTIL 

C ZERO PRESSURE IS REACHED 

C 

WRITE (2,171) 

171 FORMAT (/,' RELEASE OF SHOCKED SECONDARY DEBRIS CLOUD MATERIAL ') 
PFIN=0 . 0 

CALL RELS (CODC, KDC, RDC, GRDC, VDCO, VDC1, PH, EXTDC,UDCP, PFIN, VFDC) 

C 

c CALCULATE RELEASE OF PRESSURE WALL MATERIAL UP UNTIL INTERFACE 

C PRESSURE IS REACHED 

C 

WRITE (2,173) 

173 FORMAT(/, 'RELEASE OF SHOCKED PRESSURE WALL MATERIAL ...') 

PFIN=PHA* 1 . 0E+09 

CALL RELS ( COT , KT , RT , GRT , VTO , VT 1 , PH , EXTPW, UTP , PFIN , VFPW ) 

.... CALCULATE TEMPERATURE INCREASE IN PRESSURE WALL MATERIAL 
CALL TINC (CPST, CPLT, TMT, TVT, HFT, HVT, EXTPW) 

WRITE (*,177) 

177 FORMAT ( ' PRESSURE WALL & IMPEDANCE MATCH CALCULATIONS COMPLETE.') 
.... MODULE INTERNAL EFFECTS CALCULATION PHASE 
WRITE (*,181) 

181 FORMAT!/, ' BEGINNING MODULE INTERNAL EFFECTS CALCULATION PHASE 
$') 

C 

READ (4,183) IOPT 
183 FORMAT (II) 

C 

WRITE (2,186) 

186 FORMAT (/, '***** MODULE INTERNAL EFFECTS CALCULATION *****',//, 

$ ' CALCULATION OPTION SELECTED . . . ■ ) 

IF (IOPT. EQ. 1) WRITE (2,1871) 

1871 FORMAT (/,3X, 'SHOCK WAVE PASSAGE ONLY') 

IF (IOPT.EQ.2) WRITE (2,1872) 

1872 FORMAT (/, 3X, ' INTERNAL DEBRIS CLOUD PASSAGE ONLY') 

IF (IOPT.EQ.3) WRITE (2,1874) 

1874 FORMAT (/, 3X, 'COMBINED SHOCK WAVE AND INTERNAL DEB CLD PASSAGE') 

C 

IF (IOPT. EQ.l. OR. IOPT.EQ.3) THEN 
ESDCA=0.5*MDC* (V*1000. 0) * (V*1000. 0) 

ESDCE=0.5*MDC* (VEXP*1000 . 0 ) * (VEXP*1000. 0) 

EL1=EXTDC* (MDC*AK2) 

EL2=EXTPW* (MW*AK2 ) 

ESW=ESDCA+ESDCE- (EL1+EL2) 

SWR1=EL1/ ( ES DCA+ES DCE ) 

SWR2=EL2/ (ESDCA+ESDCE) 

SWR3=ESW/ (ESDCA+ESDCE) 

WRITE (2,306) MDC*1000 . 0, V, ESDCA, ESDCE, ESDCA+ESDCE, EXTDC, 

$ AK2*MDC* 1000 . 0, AK2 , ELI, SWR1*100 . 0, EXTPW, 

$ AK2*MW* 1000 . 0, AK2, EL2, SWR2*100 . 0 

306 FORMAT (/,3X, 'SECONDARY DEBRIS CLOUD MASS ’,F10.3, 

$' GMS ',/, 3X, ' SECONDARY DEBRIS CLOUD VELOCITY \F10.3* 

$’ KM/S' ,/,3X, 'AXIAL ENERGY OF SECONDARY DEBRIS CLOUD ... ',E10.4, 

$' JOULES' ,/,3X, 'EXPAN ENERGY OF SECONDARY DEBRIS CLOUD ... ',E10.4, 
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$' JOULES 3X, 'TOTAL ENERGY OF SECONDARY DEBRIS CLOUD ... ',E10.4, 

S' JOULES 3X, ' SDC SHOCK HEATING WASTE HEAT ',E10.4 

' JOULES/KG' , /, 3X, ' SDC MASS SUBJECTED TO SHOCK HEATING 

$F10 . 3, GMS (’, F4 .2, ’MDC) 3X, 'ENERGY LOST TO SDC SHOCK HEATING 

5 ' , E10 . 4 , ' JOULES (', F5. 2, '%) ',/,3X, ' PWALL SHOCK HEATING 

$WASTE HEAT ’,E10.4,' JOULES / KG ',/, 3X, ' PWALL MASS SUBJE 

$CTED TO SHOCK HEATING .... ',F10.3,' GMS ( ' , F4 . 2 , ’MW) ' , / , 3X, ' ENERG 

$Y LOST TO PWALL SHOCK HEATING ',E10.4,' JOULES (’,F5.2,'%) 

$' ) 


IF (IOPT.EQ.l) WRITE (2,3601) ESW, SWR3*100 . 0 

3601 FORMAT (3X, 'ENERGY TRANSFORMED INTO SW MOTION ' ,E10.4, 

$ ' JOULES ( ' , F5 . 2 , ' % ) ’ ) 

IF (ESW. LE. 0.0) THEN 
WRITE (2,990) 

WRITE (*,990) 

990 FORMAT ('*** FATAL ERROR 003: INSUFFICIENT ENERGY AVAILABLE FOR SHO 
$CK WAVE FORMATION ***') 

STOP 

ENDIF 

ENDIF 


C 


IF ( IOPT . EQ. 2 . OR. IOPT . EQ. 3 ) THEN 
ESDCA=0 . 5*MDC* (V*1000. 0) * (V*1000. 0) 

ESDCE=0 . 5*MDC* (VEXP*1000 . 0) * (VEXP*1000 . 0) 

EL1=EXTDC* (MDC*AK2 ) 

EL2=EXTPW* (MW*AK2 ) 

DDCM=MDC+MWH 
UMC=MDC* V/ DDCM 

EL3=0 . 5*DDCM* (UMC*1000 . 0 ) * (UMC*1000.0) 

DELE=ESDCA+ESDCE- (EL1+EL2+EL3 ) 

IF (DELE. LE. 0.0) THEN 
WRITE (2,992) 

WRITE (*,992) 

992 FORMAT)'*** FATAL ERROR 004: INSUFFICIENT ENERGY AVAILABLE FOR INT 
$ DEB CLD EXPANSION ***•) 

STOP 

ENDIF 

UEXP=DSQRT (DELE/ (DDCM/2.0) ) /1000. 0 
ENDIF 
C 

c READ COORDINATES OF THE POINT OF INTEREST 

C 

READ (4,1854) NPOI 

1854 FORMAT (12) 

DO 9099 IPOI=l , NPOI 
C 

c INITIALIZE TOTAL PRESSURE AND TEMPERATURE INCREASES 

C 

PINCTOT=0 . 0 
TINCTOT=0 . 0 
C 

READ (4,1855) XPOI , YPOI , ZPOI 

1855 FORMAT (3F10. 5) 

C 

WRITE (2,1856) IPOI,NPOI 

1856 FORMAT (5X, /,' ***** NOW ANALYZING POINT-OF- INTEREST NO.',I2,'/', 

$ 12 , ' *****',/) 

C 

IF (IOPT.EQ.l) CONTINUE 
IF (IOPT.EQ.2) GOTO 1000 


107 



IF (I0PT.EQ.3) THEN 

UAPOLD=UAP 

UASOLD=UAS 

ESWOLD=ESW 

EON=ESW 

EOD= (MDC+MWH) * (UAP*1000.0) * (UAP*1000.0) /2.0 
EO=EON/EOD 

UMC2= (UMC*1000. 0) * (UMC*1000. 0) 

UDE2=DELE/ (0.5* (MDC+MWH) ) 

SQ1=UMC2 

SQ2= (EO+1 ) * (UMC2-UDE2 ) 

SQ=SQ1-SQ2 
IF (SQ.LE.0.0) THEN 
WRITE (2,994) 

WRITE (*,994) 

994 FORMAT ('*** FATAL ERROR 005: INAPPROPRIATE SHOCK WAVE/DEBRIS CLOUD 
$ INTERACTION MODELING ***•) 

STOP 

ENDIF 

UAP= (UMC*1000 . 0+DSQRT ( SQ) ) / (E0+1 )/ 1000 . 0 
ESW=ESWOLD* (UAP/UAPOLD) * (UAP/UAPOLD) 

SWR3=ESW/ ( ESDCA+ESDCE) 

IF (IOPT.EQ.3) WRITE (2,3611) ESW, SWR3*100. 0 
3611 FORMAT (3X, ’ENERGY TRANSFORMED INTO SW MOTION .. ' E10 4 

$' JOULES (',F5.2, '%) ') 

UAS=COA+KA*UAP 

PHA=RA*UAS*UAP 

PHAB= ( PHA/101325 . 0 ) *1 . 0E+O9 

WRITE (2,1751) UAP,UAS, PHA, PHAB 

1751 FORMAT (/, 'MODIFIED INTERFACE VALUES ...',/, 3X, ' INTERFACE PARTICLE 
$VEL . . . ' , F8 . 5 , ' KM/S’ ,/, 3X, 'TRANS SHOCKWAVE SPEED ...',F8.5, 

$’ KM/S ',/, 3X, ' INTERFACE PRESSURE ',E12.5,' GPA ( ' F5 1 

$ • ATM) ' ) > • t 

ENDIF 

UAPT=UAP+UPW 

UAST=COA+ KA* UAPT 

PHAT=RA*UAST*UAPT 

PHATB= ( PHAT/101325 . 0 ) *1.0E+09 

WRITE (2,1752) UAPT, UAST, PHAT, PHATB 

1752 FORMAT (3X, 'TOT TRANS PARTICLE VELOCITY ',F8.5,' KM/S' / 3X 

$ ’TOT TRANS SHOCK WAVE SPEED ',F8.5,' KM/S *,/, 3X, ' TOT TRANS 

^PRESSURE ' , E12 . 5, ' GPA (',F5.1,' ATM)') 

WRITE (*,184) 

184 FORMAT (/,' BEGINNING SHOCK WAVE PASSAGE CALCULATIONS ...') 

CALL INTEFF(UAST, COA, KA, RA, E0A, PA, CPA, SWUS, SWPHB, DELT, RHA) 

WRITE (2,300) UAST, SWUS 

300 FORMAT (/,' INITIAL SHOCK WAVE PASSAGE ...',/, 3X, 'TOTAL SHOCK WAVE V 

VELOCITY ' , F10 . 3, ' KM/S ',/, 3X, ' INTERNAL SHOCK WAVE V 

$ELOCITY ' , F10 . 3, ' KM/S') 


IF (SWUS. LT. COA) THEN 
DPPOI=0 . 0 
DELTPOI=0 . 0 
WRITE (2,303) 

303 FORMAT (/, 3X, ' *** SHOCK WAVE SPEED LESS THAN SPEED OF SOUND ***') 
GOTO 309 
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ENDIF 

C 

SWPINC=SWPHB-POA 

SWTINC=DELT 

WRITE (2,3061) SWPINC, SWTINC 

3061 FORMAT <3X, ' INITAL PRESSURE INCREASE ’,F10.3, 

$’ ATM' ,/, 3X, ' INITIAL TEMPERATURE INCREASE ',F10.3, 

$ ' DEG-C ' ) 

C 

IF (IOPT.EQ. 1) JOPT=l 
IF (IOPT.EQ. 3) JOPT=31 

CALL ATTEN (AMEXP , SWPINC, COA, KA, SWUS , RA, PA, EOA, CPA, ESW, XPOI , YPOI , 

$ ZPOI , JOPT, DPPOI , DELTPOI ) 

C 

309 WRITE (*,310) 

310 FORMAT (' SHOCK WAVE PASSAGE CALCULATIONS COMPLETE.’) 

C 

IF (IOPT.EQ. 1) GOTO 9099 
C 

IF (IOPT.EQ. 3) THEN 
PINCTOT=PINCTOT+DPPOI 
TINCTOT=TINCTOT+DELTPOI 
GOTO 1000 
ENDIF 
C 

1000 WRITE (*,311) 

311 FORMAT {/,' BEGINNING DEBRIS CLOUD PASSAGE CALCULATIONS ...') 

C 

UEXP=UAP-UMC 
IF (UEXP.LE.0.0) THEN 
WRITE (2,9922) 

WRITE (*,9922) 

9922 FORMAT (//, '*** FATAL ERROR 006: INSUFFICIENT ENERGY AVAILABLE FOR 
$INT DEB CLD EXPANSION ***') 

WRITE (2,4060) MWH*1000. 0,AK3,AK2, DDCM*1000. 0, UMC, UAP, UEXP 

4060 FORMAT (3X, ' PR WALL MASS CONTRIB TO DIFF DEB CLD ',F10.3, 

$' GMS (',F4.2, '*',F4. 2, ’MW1) ',/,3X, 'DIFFUSE DEBRIS CLOUD MASS 

$ ' , F10 . 3, ' GMS' ,/, 3X, 'DIFF DEB CLD AXIAL CHARACT VELOC 

$ITY ' , F10 . 3, ' KM/S', /,3X, 'DIFF DEB CLD LEADING EDGE VELOCIT 

$ Y ' , F10 . 3, ' KM/S', /,3X, 'DIFF DEB CLD EXPANSION VELOCITY .. 

$ ’ ,F10.3, ' KM/S' ) 

STOP 

ENDIF 

C 

CALL INTEFF (UAP, COA, KA, RA, EOA, PA, CPA, DCUS, DCPHB, DELT, RHA) 

C 

WRITE (2,400) UAP, DCUS 

400 FORMAT (/, 'DEBRIS CLOUD PASSAGE ...',/, 3X, ' INT DEB CLD LEADING EDGE 

$ VELOCITY * , F10. 3, ' KM/S ' , / , 3X, ' DEB CLD INDUCED SHOCK WAV 

$E VELOCITY ' , F10 .3, ' KM/S') 

C 

IF (DCUS. LT. COA) THEN 
DPPOI=0. 0 
DELTPOI=0 . 0 
WRITE (2,403) 

403 FORMAT (/, 3X, ' *** SHOCK WAVE SPEED LESS THAN SPEED OF SOUND ***') 
GOTO 409 
ENDIF 
C 

DCPINC=DCPHB-P0A 
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DCTINC=DELT 


EDDCA=0 . 5 * DDCM* (UMC*1000.0) * (UMC* 1000.0) 

EDDCE=0 . 5*DDCM* (UEXP* 1000 . 0 ) * (UEXP*1000. 0) 

EDDC=EDDCA+EDDCE 
DDCR3=EDDCA/ (ESDCA+ESDCE) 

DDCR4=EDDCE/ (ESDCA+ESDCE) 

DDCR5=EDDC/ (ESDCA+ESDCE) 

WRITE (2,4061) DCPINC, DCTINC, MWH* 1000 . 0,AK3,AK2, DDCM* 1000 . 0, UMC, 
$ EDDCA, DDCR3* 100 . 0, UEXP, EDDCE, DDCR4* 100 . 0, EDDC, 

$ DDCR5*100 . 0 


C 

C 

C 

C 


4061 FORMAT (3X, ' INIT PRESS INCR BEH DEB CLD INDUCED SW ',F10.3, 

$' ATM' ,/, 3X, ' INIT TEMP INCR BEH DEB CLD INDUCED SW .... ^FIO*. 3, 

$' DEG-C' , /, 3X, ' PR WALL MASS CONTRIB TO DIFF DEB CLD ' ,Flo!3 

$’ GMS (', F4 .2, F4. 2, 'MWl) ',/, 3X, 'DIFFUSE DEBRIS CLOUD MASS 

f ' , F10 . 3, ' GMS ',/,3X, 'DIFF DEB CLD AXIAL CHARACT VELOC 

$ ITY ' , FI 0 . 3 , ' KM/S' ,/,3X, 'DIFFUSE DEB CLD AXIAL KINETIC ENE 

$ RGY ' , E10 . 4 , ' JOULES ( ' , F5.2, * %) ' ,/,3X, 'DIFF DEB CLD EXPAN C 

$HARACT VELOCITY ' ,F10.3,' KM/S ’,/, 3X, ’ DIFFUSE DEB CLD EXPAN 

$ KINETIC ENERGY ',E10.4,' JOULES (', F5. 2, '%)’,/, 3X, ' DIFFUSE 

$DEB CLD TOTAL KINETIC ENERGY ',E10.4,' JOULES ( ’ , F5 . 2 , ' % ) ' ) 


CALCULATE INTERNAL DEBRIS CLOUD SPREAD, ATTENUATION INITIATION 
DISTANCE, AND DEBRIS CLOUD RADIUS AT ATTENUATION ONSET 


TPHDC=UEXP/UMC 

ANG=( 180. 0/PI) *ATAN (TPHDC) 

SPHDC=DS IN ( PI *ANG/ 180.0) 

FCNPHI=SPHDC/ ( 1 . 0+SPHDC) 

DEN=PI*CD* (RA*1000. 0) *FCNPHI*FCNPHI* (UMC*1000. 0) 
$ * (UMC* 1000 . 0 ) /2 . 0 

H=100 . 0* (EDDCA/ DEN) **0.3333333333 
IF (H.GT.HMOD) H=HMOD 
RDC=H*FCNPHI 
DEL=0 . 8 * RDC* RA/ RHA 
AM=DCUS/COA 
AN=DEL+RDC 

AD=AM/DSQRT (AM*AM-1 . 0) -1 . 0 
A=AN/AD 

B=A/DSQRT (AM*AM-1 .0) 


CHECK IF POINT OF INTEREST WITHIN DEBRIS CLOUD PATH 


IF (XPOI.GE.H) GOTO 4023 
RXPOI=XPOI*TPHDC 


QCHKl=YPOI*YPOI+ZPOI*ZPOI 
QCHK2 =RX POI * RXPOI 

WRITE (2,4011) UMC, UEXP, ANG, RDC, H, AM, DEL, A, B 

4011 FORMAT (/,' INTERNAL DEBRIS CLOUD CHARACTERISTICS ’,/,3X, 'UCOM = 

$ ' , F6 . 3 , KM/S ’,/, 3X, 'UEXP = ',F6.3,' KM/S ' , /, 3X, ' 1/2-ANG = ',F6.2, 
$' DEG' , /, 3X, 'MAX RDC = ’,F6.2,' CM' , / , 3X, ' H-DIST = *,F6.2,’ CM', 
$//,'INT DEB CLD SHOCK WAVE CHARACTERISTICS ...',/, 3X, 'M# = 

$,/,3X,'DEL= ' , F6 . 2, ' CM',/,3X,’A = ',F10.2,' CM',/,3X,'B 

$F10 . 2 , ' CM' ) 


, F6.2 

t 


I F ( QCHK1 . LE . QCHK2 ) THEN 
WRITE (2,4012) XPOI, YPOI, ZPOI 

4012 FORMAT!/, '*** WARNING: PT OF INT ( ' , F6. 2, • , • , F6 . 2, • , • , F6 . 2, • ) LIES 
$ WITHIN DEB CLD CONE ***■) 

WRITE (*,410) 

GOTO 9099 
ENDIF 
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c 

c IF POINT of interest is outside debris cloud cone angle, check if 

c POINT OF INTEREST IS OVER- RUN BY FULL STRENGTH BOW SHOCK 

C 

IF (XPOI.GE.H) GOTO 4023 
Tl=XPOI- (H+DEL+A) 

RXPOI2= (Tl/A) * (Tl/A) -1 . 0 
RXPOI=DSQRT ( RXPOI2 ) 

QCHK1= (YPOI/B) * (YPOI/B) + (ZPOI/B) * (ZPOI/B) 

QCHK2=RXPOI*RXPOI 
Tll=- (H+DEL+A) 

RPWP2= (Tll/A) * (Tll/A) -1 . 0 
RPWP=DSQRT (RPWP2) 

WRITE (2,4021) B*RPWP 

4021 FORMAT (3X, ’ RPW-MIN = ',F6.2,' CM') 

IF ( QCHK1 . LE . QCHK2 ) THEN 

WRITE (2,4022) XPOI , YPOI, ZPOI 

4022 FORMAT (/,’*** WARNING: PTOFINT ( ’ , F6.2, ' , ’ , F6.2, ' , ' , F6.2, ' ) LIES 
$ OUTSIDE DEB CLD CONE ',/, 13X, ' BUT WITHIN FULL STRENGTH BOW SHOCK * 
$**') 

DPPOI=DCPHB-P0A 
DELTPOI=DCTINC 
GOTO 409 
ENDIF 
C 

c IF POINT OF INTEREST OUTSIDE DEBRIS CLOUD CONE AND OUTSIDE BOW 

c SHOCK, PROCEED WITH ATTENUATION OF PEAK PRESSURE AND TEMPERATURE 

C 

4023 IF (IOPT.EQ.2) JOPT=2 
IF (IOPT.EQ.3) JOPT=32 

CALL ATTEN (AMEXP, DCPINC, COA, KA, DCUS , RA, PA, EOA, CPA, EDDC, XPOI , YPOI 
$ , ZPOI, JOPT, DPPOI, DELTPOI ) 

C 

409 WRITE (*,410) 

410 FORMAT ( 1 DEBRIS CLOUD PASSAGE CALCULATIONS COMPLETE.') 

C 

IF (IOPT.EQ.2) GOTO 9099 

C 

IF (IOPT.EQ.3) THEN 

PINCTOT=PINCTOT+DPPOI 

TINCTOT=TINCTOT+DELTPOI 

WRITE (2,421) PINCTOT, PINCTOT* 14 . 7 , TINCTOT 

421 FORMAT (/, 'CUMULATIVE EFFECT OF SHOCK WAVE AND DEBRIS CLOUD PASSAGE 

$ ' / / , 3X , ' TOT PRESSURE INCREASE ...',F7.2,' ATM C,F9.3, 

$* PSI) ’ ,/,3X, 'TOT TEMPERATURE INCREASE ...',F7.2,' DEG-C ' ) 

PRT (IPOI, 1) =PINCTOT*14 . 7 
PRT (IPOI , 2 ) =TINCTOT 
GOTO 9099 
ENDIF 
C 

9099 CONTINUE 
C 

DO 9199 IPOI=l, NPOI 

WRITE (*,422) IPOI, PRT (IPOI, 1) , PRT(IPOI,2) 

422 FORMAT ( ' POI #',I2,5X,' TOT PRES INCR ...*,F7.3,' PSI', 5X, 'TOT TEM 
$P INCR . . . ' , F7.2, ' C ' ) 

9199 CONTINUE 
C 

CLOSE (1) 

CLOSE (2) 
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CLOSE (4) 

C 

STOP 

END 

C 

SUBROUTINE RELS (CO, K, R,G,VO,Vl, PHO, EX, UP, PHA,VF) 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DOUBLE PRECISION K, PH ( 201 ) , EH (201 ) , V(201) , P (201 ) , E (201 ) 


C 

c THIS SUBROUTINE CALCULATES THE RELEASE OF A SHOCKED MATERIAL 

c USING THE MIE-GRUNHEISEN EQUATION OF STATE. INCLUDED IS A 

C CALCULATION OF THE FINAL SPECIFIC VOLUME AND THE WASTE HEAT 

C GENERATED BY THE RELEASE PROCESS. WHEN THE PRESSURE ALONG THE 

c ISENTROPE DROPS BELOW THE REFLECTED PRESSURE ARE CALCULATED BY 

c THE IMPEDANCE MATCH PROCESS, THE RELEASE PROCESS IS TERMINATED 

C 


V(1)=V1 
PH ( 1 ) =PHO 

EH(1) =0. 5*PH (1) * (VO- VI )/ 1000. 0 
DV= ( VO- VI ) / 5 0 . 0 
DE=0 . 0 

DV2=DV/ 1000 . 0 
E ( 1 } =EH ( 1 ) 

P ( 1) =PH (1) 

DEN1=1 . 0+G*DV2*0 . 5 
C 

11=0 
UR=0 . 0 

DO 10 1=2,201 
V(I)=V(I-1) +DV 

PH (I) =CO**2*R*1000 . 0* ( 1 . 0-V(I) /VO) / ( 1. 0-K* (1. 0-V(I) /VO) ) **2 

PH (I ) =PH (I ) *1 . 0E06 

EH ( I ) =0 . 5*PH ( I ) * (VO-V ( I ) ) /1000. 0 

P(I) = (PH (I) +G* ( E ( I — 1 ) -EH (I) -0. 5*P ( I — 1 ) *DV2) ) /DENI 
E (I ) =E (I-l)-O.S* (P (I ) +P ( 1-1 ) ) *DV2 
DP=P (I) -P ( I— 1 ) 

DUR=DSQRT (-DP* (DV/1000. 0) ) 

UR=UR+DUR/ 1000.0 
11=11+1 

IF (P (I) .GE.0. 0) DE=DE+0 . 5*DV2* ( P ( I ) +P ( 1-1) ) 

IF (P(I) .LT.0.0) GOTO 15 
IF (P(I) .LE.PHA) GOTO 15 
10 CONTINUE 
C 

15 Q=P(II)/(P(II)-P(II+1) ) 

DE=DE+0 . 5*Q*DV2 * P ( 1 1 ) 

EX=EH ( 1 ) -DE 

VF=V ( II ) +Q* (V ( II+l ) -V ( II ) ) 

UFS1=UP+UR 
UFS2=2 . 0*UP 
C 

WRITE (2,20) VF , EH ( 1 ) , DE , EX 

20 FORMAT {/, 'SPECIFIC VOL AFTER RELEASE VF =',F5.3, 

$' CU. CM. /GM* ,/, 'ENERGY DUE TO DEB CLD IMPACT .... ',E10.4, 

$' JOULES/KG' ,/, 'ENERGY RECOVERED BY RELEASE ',E10.4, 

$' JOULES / KG ',/, 'WASTE HEAT GENERATED \E10.4, 

$' JOULES/KG’) 

C 

RETURN 

END 
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c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 


c 

c. 

c. 

c. 

c. 

c. 

c. 

c 


SUBROUTINE TINC ( SHS, SHL, TM, TV, HF, HV, EXH) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DOUBLE PRECISION IME,IVE 

THIS SUBROUTINE CALCULATES THE RESIDUAL TEMPERATURE INCREASE 

IN A MATERIAL THAT HAS BEEN RELEASED FROM THE SHOCKED STATE 

. ... ESTIMATES THE PERCENTAGE OF VAPORIZED, MELTED, AND SOLID 
. ... MATERIAL DUE TO THE RELEASE PROCESS 

SHS=SHSM186 . 0 
SHL=SHL*4 186.0 
HF=HF*4186 . 0 
HV=HV* 4186.0 

CALCULATE ENERGIES REQUIRED TO INITIATE MATERIAL MELT AND 

VAPORIZATION. 

IME=TM*SHS 

IVE=IME+HF+ (TV-TM) *SHL 

IF WASTE HEAT IS LESS THAN THE ENERGY REQ * D TO START MELT, 

, ... CALCULATE TEMPERATURE RISE USING W . H . =S . H . * (TEMP . INCR. ) 

IF (EXH.LT.IME) THEN 

DT=EXH/SHS 

TR=DT 

DEL=0 . 0 

WRITE (2, 50) IME, DEL, EXH 

50 FORMAT (' ENERGY REQ, INCIPIENT MELT ... *,£10. 4,' JOULE S/KG 9 ,/, 

$' ENERGY AVAILABLE FOR MELT .... 1 ,E10.4,' JOULES/KG',/, 

$ 1 EXCESS ENERGY AVAILABLE ',E10.4,' JOULES / KG' ) 

PV=0 . 0 
PL=0 . 0 
PS=100.0 
GOTO 100 
ENDIF 

IF WASTE HEAT EXCEEDS THE ENERGY REQ'D TO START MELT, BUT IS 
LESS THAN THAT REQ'D TO COMPLETE MELT, RESET THE VALUE OF THE 
ENERGY AVAILABLE FROM THE WASTE HEAT VALUE TO THE VALUE REQ'D 
... TO START MELT. THIS IMPLIES THAT SOME ENERGY IS AVAILABLE FOR 
. . . MELTING A PORTION OF THE MATERIAL. NOTE: THE TEMPERATURE RISE 
... EQUALS THE MELT TEMPERATURE OF THE MATERIAL. 

IF (EXH. GE. IME. AND. EXH. LT . IME+HF) THEN 
TR=TM 

DEL=EXH-IME 
REQM= IME+HF 

WRITE (2, 60) IME, REQM, DEL 

60 FORMAT (' ENERGY REQ, INCIPIENT MELT ... ',E10.4,' JOULES/KG *,/ , 

$' ENERGY REQ, COMPLETE MELT .... ' ,E10.4,* JOULES/ KG ',/ , 

$'ENERGY AVAILABLE FOR MELT .... ',E10.4,' JOULES/KG') 

PV=0 .0 

PL=100.0*DEL/HF 
PS=100. 0-PL 
GOTO 100 
ENDIF 


C 


IF THE WASTE HEAT EXCEEDS THE ENERGY REQ'D TO COMPLETELY MELT 
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c the material, but is less than THAT REQ'D to start vaporiza- 

c TION, compute the temperature increase caused by the excess 

c energy and add it to the melt temperature of the material. 

c 


IF (EXH.GE.IME+HF.AND.EXH.LT.IVE) then 

DEL=EXH-IME-HF 

DT=DEL/SHL 

TR=TM+DT 

REQM=IME+HF 

WRITE (2,70) IME, REQM, DEL 

70 FORMAT { 1 ENERGY REQ, INCIPIENT MELT ... ',E10.4, 1 JOULES/KG , 

$’ ENERGY REQ, COMPLETE MELT .... \ El 0.4,' JOULES/KG 1 ,/ , 

$ * EXCESS ENERGY AVAILABLE ^ElO^, 1 JOULES / KG 1 ) 

PV=0 . 0 
PL=100 . 0 
PS=100 . 0-PL 
GOTO 100 
ENDIF 
C 

IF (EXH. GE.IVE.AND.EXH.LT. IVE+HV) THEN 

DEL=EXH-IVE 

REQV= IVE+HV 

TR=TV 

WRITE (2, 80) IVE,REQV, DEL 

80 FORMAT ( * ENERGY REQ, INCIPIENT VAP \E10.4,' JOULES/KG ' , / , 

$' ENERGY REQ, COMPLETE VA? ’,E10.4,' JOULES/KG ',/ , 

$ 1 EXCESS ENERGY AVAILABLE ',£10.4,* JOULES/KG') 

PV=100. 0*DEL/HV 
PL=100.0-PV 
PS=100 . 0-PL 
GOTO 100 
ENDIF 
C 

IF (EXH. GE. IVE+HV) THEN 

ECVAP= IVE+HV 

PV=100. 0 

PL=0 . 0 

PS=0. 0 

WRITE (2,90) ECVAP 

90 FORMAT (• ENERGY REQ, COMPLETE VAP -,£10.4,- JOULES/KG 1 ,/ , 

$'+* + THE MATERIAL IS COMPLETELY VAPORIZED ***•) 

GOTO 120 
ENDIF 
C 

100 WRITE (2 , 110) TR, PS, PL, PV 

110 FORMAT (' RESIDUAL MATERIAL TEMP 1 , F10.3,' DEG-C ',//,* PERCEN 

$T SHKD AND REL PRESS WALL MATERIAL 3X, 1 IN SOLID STATE ... 1 

$, F6 . 2, * % * , /, 3X, 1 IN MOLTEN FORM ... ' , F6 . 2 , 1 % 1 , / , 3X, 1 IN VAPOR FORM 

$ 1 , F6 . 2 , ' % ’ ) 

C 

120 RETURN 
END 
C 

SUBROUTINE IMPMTCH ( RMA, RMB , KMA, KMB, COMA, COMB, VI , P4 , UPMB ) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DOUBLE PRECISION KMA, KMB 
C 

A=RMA* KMA- RMB* KMB 

B^RMA* COMA+RMB* COMB+4 . 0*RMA*KMA*V1 
Pl=2 . 0* RMA* VI* ( COMA+2 . 0* KMA* VI) 
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c 


c 


c 


c 


c 


c 

c 

c 


c 


P1=P1-P4/1 . 0E+09 

DISC=B*B-4 . 0*A*P1 

UPMB= (B-DSQRT (DISC) )/ (2.0*A) 

RETURN 

END 

SUBROUTINE INTEFF (V, CO, K, RO, EO, P, CP, US, PHB, DT, R) 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DOUBLE PRECISION K 


US=V 

UP=(US-CO) /K 
PH=P/1 . 0E+O9+RO*US*UP 
PHB= ( PH/ 101325 . 0 ) * 1 . 0E+09 
R=US*R0/ (US-UP) 

EH=E0+0 . 5* ( PH* 1 . 0E+09+P) * ( 1 . 0/R0- 1 . 0/R) / 1000 . 0 
DT= (EH-EO ) / (CPM186.0) 

RETURN 

END 

SUBROUTINE ATTEN (AM, DPI, CO, K,US, RO, PO, EO, CP, ES,XP, YP, ZP, JOPT, 

$ DPPOI , DELTPOI ) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DOUBLE PRECISION K 
COMMON/HYPERB/A, B, DEL, H 

PI=3. 141592 

IF (JOPT.EQ.l.OR. JOPT.EQ.31) RPOI=DSQRT (XP*XP+YP*YP+ZP*ZP) 

WRITE (2,5) 

5 FORMAT (/, 'PEAK SHOCK VALUES 2X, ' DIST (CM)’,3X, ' SW VEL (KM/S 

$)’,3X, 'SH DEN (GM/CC) ’ , 4X, ' PR INC (ATM) ' , 3X, ' T INC (DEG-C)') 

PB1=3. 803/ (0.01**1.2627) 

DO 100 1=1,2000 
RB=I/ 100 . 0 

IF (RB.LE.0.1) PB=3. 803/ (RB**1. 2627) 

IF ( RB . GT .0.1. AND . RB . LE . 1 . 0 ) PB=0 . 476/ ( RB**2 . 2416) 

IF (RB. GT. 1.0. AND. RB.LE. 10.0) PB=0 . 44 1/ (RB**1 . 3821 ) 

IF (P.B.GT.10.0) PB=0. 241/ (RB**1. 1120) 

RATIO=PB/PBl 

PH=(DPl*RATIO+l. 0) *101325.0 
DP=PH/101325. 0-1.0 

IF (DP. LT. 0.01. AND. JOPT. EQ.l) GOTO 200 

IF (DP. LT. 0.01. AND. JOPT.EQ.31) GOTO 200 

IF (DP. LT. 0.01. AND. JOPT. EQ. 2) GOTO 300 

IF (DP. LT. 0.01. AND. JOPT. EQ. 32) GOTO 300 

IF (I. EQ.l) THEN 

UP=(US-CO) /K 

GOTO 7 

ENDIF 

AA=R0*K 

BB=R0*CO 

CC= ( P0-PH) /I . 0E+09 

UP= (-BB+DSQRT (BB*BB-4 . 0*AA*CC) )/ (2.0*AA) 

US=CO+K*UP 
7 R=R0*US/ (US-UP) 
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E=E0+0 . 5* (PH+PO) * (l.O/RO-1. O/R) / 1000.0 
RA1=RB* ( (ES/PO) **0.3333333333) *100.0 
IF ( JOPT.EQ. l.OR. JOPT.EQ.31) AK4=1.0 
IF (JOPT.EQ. 2. OR. JOPT.EQ. 32) THEN 
IF (XP.LE.H) AK4=1 . 0 

IF (XP.GT.H) AK4=1 . 0+AM* ( (RB/0 . 01) **AM) 

ENDIF 

RA=AK4*RA1 

DELT= (E-E0) / (CP*4186. 0) 

IF (I.EQ.l) WRITE (2,10) RA, US, R, DP, DELT 
10 FORMAT ( 3X, F6 . 2 , 9X, F5 . 3, 9X, E9 . 3 , 10X, F7 . 3, 6X, F9 . 
C 

IF (JOPT.EQ. l.OR. JOPT.EQ.31) THEN 

IF (RPOI.GT.RA) GOTO 90 

IF (RPOI.EQ.RA) THEN 

DPPOI=DP 

DELTPOI=DELT 

GOTO 150 

ENDIF 

IF (RPOI.LT.RA) THEN 
Tl= ( RPOI-RAF) / (RA-RAF) 

DPPOI=DPF+ (DP-DPF) *T1 
DELTPOI=DELTF+ (DELT-DELTF) *T1 
GOTO 150 
ENDIF 
ENDIF 
C 

IF (JOPT.EQ. 2. OR. JOPT.EQ. 32) THEN 
T1=XP- (H+DEL+A+RA) 

RXP2= (Tl/A) * (Tl/A) -1 . 0 

QCHK1= (YP/B) * (YP/B) + (ZP/B) * (ZP/B) 

QCHK2=RXP2 

XFS=H+DEL+A-A*DSQRT (QCHK1+1. 0) 

IF ( QCHK1 . GT . QCHK2 ) GOTO 90 

IF ( QCHK1 . EQ . QCHK2 ) THEN 

DPPOI=DP 

DELTPOI=DELT 

GOTO 150 

ENDIF 

IF ( QCHK1 . LT . QCHK2 ) THEN 
Tl= (XP- (XFS+RAF) )/ (RA-RAF) 

DPPOI=DPF+ (DP-DPF) *T1 
DELTPOI=DELTF+ (DELT-DELTF) *T1 
GOTO 150 
ENDIF 
ENDIF 
C 

90 RAF=RA 
USF=US 
RF=R 
DPF=DP 
DELTF=DELT 
C 

100 CONTINUE 
C 

150 WRITE (2,15) XP,YP,ZP,AM 
15 FORMAT (/, ' POINT-OF-INTEREST COORDINATES 

$F6 . 2 , ' CM' , / , 3X, ' Y-COORD = ',F6.2,' CM',/,3X,* 
$' CM' ,//,3X, 'M-CONST = ’,F6.2) 

IF (JOPT.EQ. l.OR. JOPT.EQ.31) WRITE (2,16) RPOI 


3X, 'X-COORD = ' , 
-COORD = ' , F6.2, 
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16 FORMAT (3X, ' R-DI ST = ',F6.2,' CM') 

IF ( J0PT.EQ.2.0R. JOPT.EQ.32) WRITE (2,161) XP-XFS 
161 FORMAT (3X, 'ETA- DST = ',F6.2,' CM') 

WRITE (2,17) US, R, DPPOI, DPPOI*14 . 7, DELTPOI 

17 FORMAT (/, 'ATTENUATED SHOCK VALUES 3X, ' SH WAVE VEL = ',F9.3, 

$' KM/S' ,/,3X, 'SH DENSITY = ',E13.3,' GM/CC 3X, * PRESS INCR = ' 
$ , F9 . 3 , ' ATM ( ' , F9 . 3 , ' PSI ) ' , /, 3X ' TEMP INCR = ',F9.3,' DEG-C') 

RETURN 

200 WRITE (2,10) RAF, USF, RF, DPF, DELTF 
DPPOI=0. 0 
DELTPOI=0 . 0 

WRITE (2,205) XP,YP,ZP 

205 FORMAT (/, ' POINT-OF-INTEREST COORDINATES 3X, 'X-COORD = ', 

$F6. 2, ' CM' ,/, 3X, 'Y-COORD = ' , F6 . 2 , ' CM' , / , 3X, ' Z-COORD = ',F6.2, 

$’ CM') 

WRITE (2,210) RPOI , RAF 

210 FORMAT('*** NOTE: DIST TO PT OF INT (',F6.2, ' CM) BEYOND DP=1%ATM 
$DIST ( ' , F6.2, 'CM) *** ' ) 

RETURN 

300 WRITE (2,10) RAF, USF, RF, DPF, DELTF 
DPPOI=0 . 0 
DELTPOI=0 . 0 

WRITE (2,305) XP,YP,ZP 

305 FORMAT (/,' POINT-OF-INTEREST COORDINATES 3X, 'X-COORD = ', 

$F7 . 2 , ' CM ',/,3X, 'Y-COORD = ',F7.2,' CM' ,/, 3X, ' Z-COORD = ',F7.2, 

$' CM' ) 

AP=A+RAF 

QCHK22= (H+DEL+AP) * (H+DEL+AP) / (A*A) -1.0 
QCHK2=DSQRT (QCHK22) 

QCHK12= ( YP/B) * (YP/B) + (ZP/B) * (ZP/B) 

QCHK1=DSQRT (QCHK12 ) 

RPW=B*QCHK2 

IF ( QCHK1 . GE . QCHK2 ) WRITE (2,310) RPW 
310 FORMAT (’*** NOTE: PT OF INT OUTSIDE MAX SHOCK WAVE CYL (RPW-MAX = 

$ ' , F6. 2, ' CM) **+' ) 

IF ( QCHK1 . LT . QCHK2 ) WRITE (2,315) XP-XFS, XP,XFS, RAF 
315 FORMAT ( ' *** NOTE: X-DIR DIST (',F6.2,' CM) TO PT OF INT (X-COORD = 
$ ' , F6.2, ' CM) FROM' ,/, 10X, 'FULL STRENGTH SHOCK (X-COORD = ',F6.2, 
$’ CM) BEYOND DP=1%ATM' , /, 10X, ' DIST (',F6.2,' CM) ***') 

RETURN 

END 

SUBROUTINE ROOTS (AK1, R2 ,MIB,MDC, MDC1, Cl, C2 ) 

IMPLICIT DOUBLE PRECISION (A-H, O-Z ) 

DOUBLE PRECISION MIB,MDC,MDC1 

R1=MIB/MDC1 

T0= (Rl+1 . 0 ) *AK1*AK1 

DD=(MDC*R2) * (MDC*R2 ) *R1* (TO— 1 . 0) 

T11=MDC*R2*R1 

T12=DSQRT(DD) 

C2N=T11-T12 
C2D=MIB* ( Rl+1 . 0 ) 

C2=C2N/C2D 

Cl= (MDC*R2-MIB*C2 ) /MDC1 

RETURN 

END 
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APPENDIX B - SAMPLE INPUT FILES 
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INPUT FILE: IMPDAT 



MAT 1 L — | 

— EL. MOD . - | 
— T .MELT — | 

CO 

NU 

--T.VAP-- 

1 K | 

1 —ALFA 1 

1— H.FUS— | 

RHO | 

CPS | 

— H . VAP — | 

— GAMO | 

CPL | 

— 

AL 

ALUMINUM 
0. 103E+08 
660.0 

5.380 

0.35 

2450.0 

1.340 
0 . 240E-04 
95.0 

2.712 

0.235 

2450.0 

2.130 

0.255 

— 

A1 

2 XXX ALUM 
0. 106E+08 
640.0 

5.350 

0.33 

2450.0 

1.340 
0 . 209E-Q4 
85.0 

2.800 

0.212 

2450.0 

2.000 

0.242 


A2 

5 XXX ALUM 
0 . 101E+08 
641.0 

5.310 

0.33 

2450.0 

1.340 
0 . 225E-04 
85.0 

2.670 

0.215 

2450.0 

2.000 

0.245 

— 

A3 

6XXX ALUM 
0 . 100E+08 
652.0 

5.380 

0.33 

2450.0 

1.340 
0 . 233E-04 
85.0 

2.700 

0.212 

2450.0 

2.000 

0.242 

— 

A4 

7 XXX ALUM 
0. 103E+08 
636.0 

5.290 

0.33 

2450.0 

1.340 
0 . 221E-04 
85.0 

2.810 

0.217 

2450.0 

2.000 

0.245 


BE 

BERYLLIUM 
0 . 419E+08 
1281.0 

7.975 

0.08 

2884.0 

1.124 
0. 140E-04 
260.0 

1.820 

0.570 

8195.0 

1.160 

0.832 

— 

CD 

CADMIUM 
0 . 672E+07 
321.0 

2.307 

0.33 

765.0 

1.640 

0.343E-04 

13.5 

8.640 

0.058 

212.0 

2.270 

0.063 

— 

cu 

COPPER 
0 . 190E+08 
1083.0 

3.940 

0.34 

2590.0 

1.489 
0. 170E-04 
49.0 

8.930 

0.097 

1150.0 

2.000 

0.114 


EP 

EPOXY 
0. 650E+06 
350.0 

3.020 

0.50 

-1.0 

1.520 
0. 500E-04 
-1.0 

1.180 

0.250 

-1.0 

0.800 

0.285 

— 

FE 

IRON 

0 . 290E+08 
1539.0 

4.580 

0.30 

3035.0 

1.490 
0 . 120E-04 
65.0 

7.870 

0.120 

1591.0 

1.570 

0.150 
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PB 


LEAD 

0.200E+07 

327.0 

2.030 

0.45 

1740.0 

1.470 
0 , 293E-04 
6.0 

11.340 

0.031 

210.0 

2.770 

0.036 

LX 

LEXAN 
0 . 345E+06 
225.0 

2.750 

0.50 

-1.0 

1.480 
0. 650E-04 
-1.0 

1.180 

0.290 

-1.0 

0.860 

0.315 

MO 

MOLYBDENUM 
0 . 460E+08 
2610.0 

5.173 

0.31 

5555.0 

1.220 
0. 061E-04 
70.0 

10.200 

0.079 

1242.0 

1.520 

0.104 

NI 

NICKEL 
0 . 330E+08 
1454.0 

4.667 

0.30 

2865.0 

1.530 
0. 143E-04 
74.0 

8.860 

0.130 

1523.0 

1.800 

0.157 

PT 

PLATINUM 
0 * 277E+08 
1769.0 

3.680 

0.39 

4349.0 

1.500 
0. 110E-04 
26.0 

21.370 

0.037 

632.0 

2.940 

0.042 

SI 

304 STEEL 
0.284E+08 
1425.0 

4.590 

0.28 

3035.0 

1.550 
0. 112E-04 
65.0 

7.910 

0.110 

1590.0 

1.670 

0.125 

S2 

430 STEEL 
0 . 299E+08 
1470.0 

4.680 

0.29 

3035.0 

1.550 
0. 104E-04 
65.0 

7.830 

0.110 

1590.0 

1.670 

0.125 

S3 

4340 STEEL 
0 . 290E+08 
1510.0 

4.570 

0.30 

3070.0 

1.550 
0. 112E-04 
65.0 

7.830 

0.110 

1590.0 

1.670 

0.125 

TA 

TANTALUM 
0 . 260E+08 
2996.0 

3.374 

0.35 

5425.0 

1.201 
0. 065E-04 
38.0 

16.650 

0.033 

1007.0 

1.690 

0.039 

SN 

TIN 

0 . 603E+07 
235.0 

2.560 

0.33 

2450.0 

1.520 
0 . 269E-04 
14.0 

7.280 

0.058 

580.0 

1.850 

0.062 

TI 

TITANIUM 
0 . 180E+08 
1676.0 

4.786 

0.30 

3260.0 

1.049 
0. 100E-04 
99.0 

4.512 

0.150 

2182.0 

1.100 

0.167 

W 

TUNGSTEN 
0 . 590E+08 
3410.0 

4.150 

0.30 

5900.0 

1.237 
0 . 040E-04 
53.0 

19.170 

0.035 

1054.0 

1.480 

0.046 
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ZINC 

0 . 108E+08 
420.0 

3.042 

0.33 

907.0 

1.500 
0 . 274E-04 
25.0 

7.140 

0.100 

420.0 

2.150 

0.115 

AU 

GOLD 

0 . 124E+08 
1063.0 

3.060 

0.42 

2960.0 

1.570 
0. 161E-04 
16.0 

19.240 

0.034 

413.0 

3.100 

0.038 

AG 

SILVER 
0 . 120E+08 
961.0 

3.230 

0.37 

2210.0 

2.500 
0 . 211E-04 
25.0 

10.490 

0.062 

554.0 

2.500 

0.071 

MG 

MAGNESIUM 
0 . 640E+07 
650.0 

4.490 

0.29 

1110.0 

1.240 

0.300E-04 

88.0 

1.740 

0.295 

1326.0 

1.500 

0.336 

XX 

INPUT FILE: 

GPARAM 



AL 

AL 

AL 

AA 

1.27 

0.19 

0.16 

0.30 

0.48 

0.50 

11.43 

0.60 

5.72 


1 

3 

4 


18.0 

61.0 

0.0 

119.0 

61.0 

0.0 

203.0 

61.0 

0.0 

203.0 

0.0 

0.0 


0.20 
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APPENDIX C - SAMPLE OUTPUT FILE 
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NORMAL IMPACT OF A ALUMINUM PROJECTILE ON A DUAL-WALL SYSTEM WITH A 
ALUMINUM BUMPER, A MLI-BLNKT INNER BUMPER, AND A ALUMINUM PRESSURE WALL 

PROJECTILE PROPERTIES . . . 


MAT 

as 

ALUMINUM 

CO 

= 

5.380 

KM/S 

K 

= 

1.340 


RHO 

= 

2.712 

GM/ CU.CM. 

DP 

= 

1.270 

CM 

MP 

- 

2.909 

GMS 

VP 

= 

6.700 

KM/S 

BUMPER 

PROPERTIES . . . 

MAT 

= 

ALUMINUM 

CO 

= 

5.380 

KM/S 

K 

= 

1.340 


RHO 

= 

2.712 

GM/ CU.CM. 

TS 

= 

. 160 

CM 

DH 


2.257 

CM 

MB 

= 

1.736 

GMS 

**** PROJECTILE 

AND BUMPER SHOCK LOADING RESPONSE AND RELEASE CALCULATIONS 




PROJECTILE IMPACT VELOCITY 

PROJ MATL PARTICLE VELOCITY . . . 
PROJ MATL SHOCK WAVE SPEED 

HUGONIOT IMPACT PRESSURE 

BMPR MATL PARTICLE VELOCITY . . . 

BMPR MATL SHOCK WAVE SPEED 

HUGONIOT IMPACT PRESSURE 


VP 

= 

6.700 

KM/S 

UP 

= 

3.350 

KM/S 

US 

as 

9.869 

KM/S 

PH 

= 

89.662 

GPA 

UP 

= 

3.350 

KM/S 

US 

= 

9.869 

KM/S 

PH 

= 

89.662 

GPA 


PARAMETERS REQUIRED FOR CALCULATING PROJECTILE MATERIAL RESPONSE AND 
RELEASE FROM SHOCKED STATE USING THE MIE-GRUNEISEN E-O-S: 


ELASTIC MODULUS 

. . E 

= 

. 7102E+11 

N/SQ.M. 

POISSON RATIO 

. . NU 

= 

.350 


BULK MODULUS 

. . K 

= 

. 7891E+11 

N/SQ.M. 

LIN. COEF. OF THERM. EXP. 

. . ALFA 

= 

. 2400E-04 

/DEG-C 

SP HEAT (SOLID) 

. . CPS 

= 

.235 

CAL/GM/DEG-C 

SP HEAT (LIQUID) 

. . CPL 

= 

.255 

CAL/GM/DEG-C 

HUGON IMP PRESS ( PA,MBAR) 

. . PH 

= 

. 8966E+11, 

.897 

SP VOL AT REST 

. . VO 

= 

.369 

CU.CM./GM 

SP VOL AT IMPACT 

. . VI 

= 

.244 

CU.CM./GM 

AMB M-GRUN COEF (CAL,INP) 

. . GAM0 

= 

2 . 129, 

2.130 

MELT TEMPERATURE 

. . TM 

= 

660.00 

DEG-C 

VAPOR TEMPERATURE 

. . TV 

= 

2450.00 

DEG-C 

HEAT OF FUSION 

. . HF 

= 

95.00 

CAL/GM 

HEAT OF VAPORIZATION 

. . HV 

= 

2450.00 

CAL/GM 


RELEASE OF SHOCKED PROJECTILE MATERIAL . . . 


SPECIFIC VOL AFTER RELEASE 

ENERGY DUE TO DEB CLD IMPACT . , . 

ENERGY RECOVERED BY RELEASE 

WASTE HEAT GENERATED 

ENERGY REQ, INCIPIENT MELT . . . 

ENERGY REQ, COMPLETE MELT 

ENERGY AVAILABLE FOR MELT 

RESIDUAL MATERIAL TEMP 


VF = .399 CU.CM./GM 
. 5611E+07 JOULES/KG 
. 4616E+07 JOULES/KG 
. 9957E+06 JOULES/KG 
6492E+06 JOULES/KG 
1047E+07 JOULES/KG 
3464E+06 JOULES/KG 
660.000 DEG-C 
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PERCENT SHKD AND REL PRESS WALL MATERIAL 
IN SOLID STATE ... 12.89% 

IN MOLTEN FORM ... 87.11% 

IN VAPOR FORM .... . 00% 


PARAMETERS REQUIRED FOR CALCULATING BUMPER MATERIAL RESPONSE AND 


RELEASE FROM SHOCKED STATE USING 

ELASTIC MODULUS 

POISSON RATIO 

BULK MODULUS 

LIN. COEF. OF THERM. EXP. ... 

SP HEAT (SOLID) 

SP HEAT (LIQUID) 

HUGON IMP PRESS ( PA, MBAR ) ... 

SP VOL AT REST 

SP VOL AT IMPACT 

AMB M-GRUN COEF (CAL,INP) ... 

MELT TEMPERATURE 

VAPOR TEMPERATURE 

HEAT OF FUSION 

HEAT OF VAPORIZATION 


THE MIE-GRUNEISEN E-O-S: 

E = . 7102E+11 N/SQ.M. 

NU = .350 

K = . 7891E+11 N/SQ.M. 

ALFA = . 2400E-04 /DEG-C 
CPS = .235 CAL/GM/DEG-C 

CPL = .255 CAL/GM/DEG-C 

PH = . 8966E+11, .897 
VO = .369 CU.CM. /GM 

VI = .244 CU.CM. /GM 

GAMO = 2.129,2.130 

TM = 660.00 DEG-C 

TV = 2450.00 DEG-C 

HF = 95.00 CAL/GM 

HV = 2450.00 CAL/GM 


RELEASE OF SHOCKED BUMPER MATERIAL . . . 


SPECIFIC VOL AFTER RELEASE VF = .399 CU.CM./GM 

ENERGY DUE TO DEB CLD IMPACT 5611E+07 JOULES/KG 

ENERGY RECOVERED BY RELEASE 4616E+07 JOULES/KG 

WASTE HEAT GENERATED 9957E+06 JOULES/KG 

ENERGY REQ, INCIPIENT MELT ... .6492E+06 JOULES/KG 

ENERGY REQ, COMPLETE MELT 1047E+07 JOULES/KG 

ENERGY AVAILABLE FOR MELT 3464E+06 JOULES/KG 

RESIDUAL MATERIAL TEMP 660.000 DEG-C 

PERCENT SHKD AND REL PRESS WALL MATERIAL . . . 

IN SOLID STATE ... 12.89% 

IN MOLTEN FORM ... 87.11% 

IN VAPOR FORM .00% 


PRIMARY DEBRIS CLOUD CHARACTERISTIC VELOCITIES . . . 
VEXP = 2.918 KM/S 

VCOM = 4.196 KM/S 

VLE = 6.700 KM/S 

VPW = .118 KM/S 

SELECTED PRIMARY DEBRIS CLOUD IMPACT VELOCITY 
VDC = 4.196 KM/S 


PRIMARY DEBRIS CLOUD 
1/2-ANG = 34.821 


RATIO- P = .626 
RPSR = 2.504 
RATIO-B = .374 
RBSR = 2.504 
CO = 5.380 
K = 1.340 
RHO = 2.504 
MDC = 4.645 


PROPERTIES . . . 
DEG 

GM/ CU.CM. 

GM/ CU.CM. 

KM/S (AVG) 

(AVG) 

GM/ CU.CM. (AVG) 
GMS 


INNER BUMPER PROPERTIES . . . 
RHO-A = .072 GM/SQ.CM. 

DH = 6.293 CM 
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S2 = 5.720 CM 

SECONDARY DEBRIS CLOUD PROPERTIES (GENERAL) ... 

SEC DEB CLD AVG AXL VEL 2.885 KM/S 

PRIM DEB CLD AXL ENERGY 4088E+05 J 

PRIM DEB CLD EXP ENERGY 1978E+05 J 

SEC DEB CLD AXL ENERGY 2812E+05 J 

PRIM DEB CLD AXL EN LOST ... . 1147E+04 J 

PRIM DEB CLD EXP EN LOST ... .5548E+03 J 

EN AVAIL SEC DEB CLD EXP ... .3084E+05 J 

SEC DEB CLD AVG EXP VEL 3.022 KM/S 

ADDL PRESS WALL VEL 308 KM/S 

SECONDARY DEBRIS CLOUD PROPERTIES (1ST COMP) ... 

1/2-ANG = 34.821 DEG 
RHO = 2.504 GM/ CU.CM. 

N = .500 

K1 = .972 (> .788) 

MDC1 = 4.515 GMS 

C#1 = 1.453 

C#2 = .965 

VMC = 4.192 KM/S 

VEXP = 2.916 KM/S 

SECONDARY DEBRIS CLOUD PROPERTIES (2ND COMP) ... 

1/2-ANG = 34.821 DEG 
RHO-A = .072 GM/ CU.CM. 

MIB = 2.239 GMS 

C#1 = .087 

C#2 = .058 

VMC = .252 KM/S 

VEXP = .175 KM/S 

PRESSURE WALL PROPERTIES . . . 

MAT = ALUMINUM 

CO = 5.380 KM/S 

K = 1.340 

RHO = 2.712 GM/ CU.CM. 

TW = .480 CM 

S = 11.430 CM 

MODULE AIR PROPERTIES . . . 

CO = .330 KM/S 

K = 1.070 

RHO = . 122E-02 GM/ CU.CM. 

CP = .240 CAL/GM/DEG-C 

TO = 24.000 DEG-C 
P0 = 1.000 ATM 

MODULE DIAMETER . . . 

HMOD = 400.00 CM 

*+** PRESSURE WALL RESPONSE AND IMPEDANCE MATCH CALCULATIONS **** 

ASSIGNED PRESS WALL IMP VEL VW = 4.192 KM/S 

PRESSURE WALL PARTICLE VELOCITY ... UP = 2.096 KM/S 

PRESSURE WALL SHOCK WAVE SPEED US = 8.188 KM/S 

HUGONIOT IMPACT PRESSURE PH = 46.543 GPA 

PARAMETERS REQUIRED FOR CALCULATING PRESSURE WALL MATERIAL RESPONSE AND 
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RELEASE 

FROM SHOCKED STATE USING 

THE MIE-GRUNEISEN E-O- 

S: 

MATL 

ELASTIC MODULUS 

E 

= 

. 7102E+11 

N/SQ.M. 

MATL 

POISSON RATIO 

NU 

= 

.350 


MATL 

BULK MODULUS 

K 

= 

. 7891E+11 

N/SQ.M. 

MATL 

LIN. COEF. OF THERM. EXP. 

. , * ALFA 

= 

.2400E-04 

/DEG-C 

MATL 

SP HEAT (SOLID) 

CPS 

= 

.235 

CAL/ GM/ DEG 

MATL 

SP HEAT (LIQUID) 

CPL 

= 

.255 

CAL/ GM/ DEG 

MATL 

HUGON IMP PRESS ( PA, MBAR) 

... PH 

= 

. 4654E+11, 

.465 

MATL 

SP VOL AT REST 

VO 

= 

.369 

CU.CM./GM 

MATL 

SP VOL AT IMPACT 

VI 

= 

.274 

CU.CM./GM 

MATL 

AMB M-GRUN COEF (CAL,INP) 

. . . GAM0 

= 

2.129 f 

2.130 

MATL 

MELT TEMPERATURE 

TM 

= 

660.00 

DEG-C 

MATL 

VAPOR TEMPERATURE 

TV 

= 

2450.00 

DEG-C 

MATL 

HEAT OF FUSION 

HF 

= 

95.00 

CAL/GM 

MATL 

HEAT OF VAPORIZATION 

HV 

=: 

2450.00 

CAL/GM 

RESULTS 

OF IMPEDANCE MATCH CALCULATIONS . . . 





INTERFACE ( REFL & TRANS) PARTICLE VEL ... 4.19006 KM/S 

TRANSMITTED SHOCK WAVE SPEED 4.81337 KM/S 

INTERFACE (REFL & TRANS) PRESSURE 24605E-01 GPA (242.8 ATM) 

RELEASE OF SHOCKED SECONDARY DEBRIS CLOUD MATERIAL . . . 

SPECIFIC VOL AFTER RELEASE VF = .414 CU.CM./GM 

ENERGY DUE TO DEB CLD IMPACT 2379E+07 JOULES/KG 

ENERGY RECOVERED BY RELEASE 1950E+07 JOULES/KG 

WASTE HEAT GENERATED 4289E+06 JOULES/KG 

RELEASE OF SHOCKED PRESSURE WALL MATERIAL . . . 


SPECIFIC VOL AFTER RELEASE . . 
ENERGY DUE TO DEB CLD IMPACT 
ENERGY RECOVERED BY RELEASE 

WASTE HEAT GENERATED 

ENERGY REQ, INCIPIENT MELT . 
ENERGY AVAILABLE FOR MELT . . 

EXCESS ENERGY AVAILABLE 

RESIDUAL MATERIAL TEMP 


VF = .379 CU.CM./GM 
. 2196E+07 JOULES/KG 
. 1847E+07 JOULES/KG 
. 3496E+06 JOULES/KG 
6492E+06 JOULES/KG 
0000E+00 JOULES/KG 
3496E+06 JOULES/KG 
355.374 DEG-C 


PERCENT SHKD AND REL PRESS WALL MATERIAL ... 

IN SOLID STATE ... 100.00% 

IN MOLTEN FORM ... .00% 

IN VAPOR FORM 00% 

***** MODULE INTERNAL EFFECTS CALCULATION ***** 

CALCULATION OPTION SELECTED . . . 


COMBINED SHOCK WAVE AND INTERNAL DEB CLD PASSAGE 


SECONDARY DEBRIS CLOUD MASS 

SECONDARY DEBRIS CLOUD VELOCITY 

AXIAL ENERGY OF SECONDARY DEBRIS CLOUD . . . 
EXPAN ENERGY OF SECONDARY DEBRIS CLOUD . . . 
TOTAL ENERGY OF SECONDARY DEBRIS CLOUD . . . 

SDC SHOCK HEATING WASTE HEAT 

SDC MASS SUBJECTED TO SHOCK HEATING 

ENERGY LOST TO SDC SHOCK HEATING 

PWALL SHOCK HEATING WASTE HEAT 

PWALL MASS SUBJECTED TO SHOCK HEATING 


4.515 GMS 
4.192 KM/S 
. 3966E+05 JOULES 
. 1919E+05 JOULES 
. 5085E+O5 JOULES 
. 4289E+06 JOULES/KG 

.858 GMS ( . 19MDC) 

. 3679E+03 JOULES ( .63%) 

. 3496E+06 JOULES/KG 
9.390 GMS ( . 19MW) 
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ENERGY LOST TO PWALL SHOCK HEATING 


. 3283E+04 JOULES ( 5.58%) 


***** NOW ANALYZING POINT-OF-INTEREST NO. 1/4 ***** 

ENERGY TRANSFORMED INTO SW MOTION 2986E+05 JOULES (50.74%) 

MODIFIED INTERFACE VALUES . . . 

INTERFACE PARTICLE VEL ... 3.08170 KM/S 
TRANS SHOCK WAVE SPEED ... 3.62742 KM/S 

INTERFACE PRESSURE 13638E-01 GPA (134.6 ATM) 

TOT TRANS PARTICLE VELOCITY 3.38979 KM/S 

TOT TRANS SHOCK WAVE SPEED 3.95708 KM/S 

TOT TRANS PRESSURE 16365E-01 GPA (161.5 ATM) 


INITIAL SHOCK WAVE PASSAGE . . . 

TOTAL SHOCK WAVE VELOCITY 3.957 KM/S 

INTERNAL SHOCK WAVE VELOCITY 3^957 KM/S 

INITAL PRESSURE INCREASE 161.507 ATM 

INITIAL TEMPERATURE INCREASE 5789.626 DEG-C 


PEAK SHOCK VALUES . . . 

DIST (CM) SW VEL (KM/S) SH DEN (GM/CC) 
•67 3.957 . 851E-02 


PR INC (ATM) T INC (DEG-C) 
161.507 5789.626 


POINT-OF-INTEREST COORDINATES . . . 
X-COORD = 18.00 CM 

Y-COORD = 61.00 CM 

Z -COORD = .00 CM 


M-CONST = . 60 

R-DIST = 63.60 CM 


ATTENUATED SHOCK 
SH WAVE VEL = 
SH DENSITY = 
PRESS INCR = 
TEMP INCR 


VALUES . . . 

.347 KM/S 
. 128E-02 GM/CC 
.067 ATM ( 
3.931 DEG-C 


.981 PSI) 


DEBRIS CLOUD PASSAGE . . . 

INT DEB CLD LEADING EDGE VELOCITY 

DEB CLD INDUCED SHOCK WAVE VELOCITY 

INIT PRESS INCR BEH DEB CLD INDUCED SW . . . 

INIT TEMP INCR BEH DEB CLD INDUCED SW 

PR WALL MASS CONTRIB TO DIFF DEB CLD 

DIFFUSE DEBRIS CLOUD MASS 

DIFF DEB CLD AXIAL CHARACT VELOCITY 

DIFFUSE DEB CLD AXIAL KINETIC ENERGY 

DIFF DEB CLD EXPAN CHARACT VELOCITY 

DIFFUSE DEB CLD EXPAN KINETIC ENERGY 

DIFFUSE DEB CLD TOTAL KINETIC ENERGY 


3.082 

KM/S 


3.082 

KM/S 


95.422 

ATM 


3360.483 

DEG-C 


2.817 

GMS ( . 

30* . 19MW1 ) 

7.332 

GMS 


2.581 

KM/S 


2442E+05 

JOULES 

(41.50%) 

.501 

KM/S 


9183E+03 

JOULES 

( 1.56%) 

2534E+05 

JOULES 

(43.06%) 


INTERNAL DEBRIS CLOUD CHARACTERISTICS 
UCOM = 2.581 KM/S 

UEXP = .501 KM/S 

1/2-ANG = 10.97 DEG 

MAX RDC = 11.52 CM 

H-DIST = 72.05 CM 


INT DEB CLD SHOCK WAVE CHARACTERISTICS . . . 
M# = 9.34 
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DEL = 1.53 CM 

A = 2256.06 CM 

B = 242.98 CM 

RPW-MIN = 62.56 CM 

PEAK SHOCK VALUES ... 

DIST (CM) SW VEL (KM/S) SH DEN (GM/CC) PR INC (ATM) T INC (DEG-C) 
•63 3.082 . 737E-Q2 95.422 3360.483 

POINT-OF- INTEREST COORDINATES . . . 

X-COORD = 18.00 CM 

Y-COORD = 61.00 CM 

Z-COORD = .00 CM 

M-CONST = .60 

ETA-DST = 14.43 CM 

ATTENUATED SHOCK VALUES . . . 

SH WAVE VEL = .501 KM/S 

SH DENSITY = . 179E-02 GM/CC 

PRESS INCR = .970 ATM ( 14.258 PSI) 

TEMP INCR = 39.274 DEG-C 

CUMULATIVE EFFECT OF SHOCK WAVE AND DEBRIS CLOUD PASSAGE 

TOT PRESSURE INCREASE ... 1.04 ATM ( 15.239 PSI) 

TOT TEMPERATURE INCREASE ... 43.20 DEG-C 

***** NOW ANALYZING POINT-OF-INTEREST NO. 2/ 4 ***** 

ENERGY TRANSFORMED INTO SW MOTION 2986E+05 JOULES (50.74%) 

MODIFIED INTERFACE VALUES . . . 

INTERFACE PARTICLE VEL ... 3.08170 KM/S 
TRANS SHOCK WAVE SPEED ... 3.62742 KM/S 

INTERFACE PRESSURE 13638E-01 GPA (134.6 ATM) 

TOT TRANS PARTICLE VELOCITY 3.38979 KM/S 

TOT TRANS SHOCK WAVE SPEED 3.95708 KM/S 

TOT TRANS PRESSURE 16365E-01 GPA (161.5 ATM) 

INITIAL SHOCK WAVE PASSAGE . . . 

TOTAL SHOCK WAVE VELOCITY . . 

INTERNAL SHOCK WAVE VELOCITY 
INITAL PRESSURE INCREASE . . . 

INITIAL TEMPERATURE INCREASE 

PEAK SHOCK VALUES . . . 

DIST (CM) SW VEL (KM/S) SH DEN (GM/CC) PR INC (ATM) T INC (DEG-C) 
•67 3.957 . 851E-02 161.507 5789.626 

POINT-OF-INTEREST COORDINATES . . . 

X-COORD = 119.00 CM 
Y-COORD = 61.00 CM 

Z -COORD = .00 CM 

M-CONST = .60 

R-DIST = 133.72 CM 

ATTENUATED SHOCK VALUES . . . 

SH WAVE VEL = .336 KM/S 

SH DENSITY = .124E-02 GM/CC 


3.957 KM/S 
3.957 KM/S 
161.507 ATM 
5789.626 DEG-C 
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PRESS INCR 
TEMP INCR 


.021 ATM ( 
1.312 DEG-C 


.313 PSI) 


DEBRIS CLOUD PASSAGE . . . 

INT DEB CLD LEADING EDGE VELOCITY 

DEB CLD INDUCED SHOCK WAVE VELOCITY ...... 

INIT PRESS INCR BEH DEB CLD INDUCED SW . . . 

INIT TEMP INCR BEH DEB CLD INDUCED SW 

PR WALL MASS CONTRIB TO DIFF DEB CLD 

DIFFUSE DEBRIS CLOUD MASS 

DIFF DEB CLD AXIAL CHARACT VELOCITY 

DIFFUSE DEB CLD AXIAL KINETIC ENERGY 

DIFF DEB CLD EXPAN CHARACT VELOCITY 

DIFFUSE DEB CLD EXPAN KINETIC ENERGY 

DIFFUSE DEB CLD TOTAL KINETIC ENERGY 

PEAK SHOCK VALUES . . . 

DIST (CM) SW VEL (KM/S) SH DEN (GM/CC) 
1.01 3.082 . 737E-02 

POINT-OF- INTEREST COORDINATES . . . 

X- COORD = 119.00 CM 
Y-COORD = 61.00 CM 

Z-COORD = .00 CM 


3.082 
3.082 
95.422 
3360.483 
2.817 
7.332 
2.581 
.2442E+05 
.501 
. 9183E+03 
.2534E+05 


KM/S 

KM/S 

ATM 

DEG-C 

GMS ( , 

GMS 

KM/S 

JOULES 

KM/S 

JOULES 

JOULES 


30* . 19MW1 ) 


(41.50%) 

( 1.56%) 
(43.06%) 


PR INC (ATM) T INC (DEG-C) 
95.422 3360.483 


M-CONST = . 60 

ETA-DST = 115.43 CM 


ATTENUATED SHOCK 
SH WAVE VEL = 
SH DENSITY = 
PRESS INCR = 
TEMP INCR 


VALUES . . . 

.426 KM/S 
. 154E-02 GM/CC 
.468 ATM ( 
21.694 DEG-C 


6.872 PSI) 


CUMULATIVE EFFECT OF SHOCK WAVE AND DEBRIS CLOUD PASSAGE 

TOT PRESSURE INCREASE ... .49 ATM ( 7.185 PSI) 

TOT TEMPERATURE INCREASE ... 23.01 DEG-C 

***** NOW ANALYZING POINT-OF-INTEREST NO. 3/4 ***** 

ENERGY TRANSFORMED INTO SW MOTION 2986E+05 JOULES (50.74%) 

MODIFIED INTERFACE VALUES . . . 

INTERFACE PARTICLE VEL ... 3.08170 KM/S 
TRANS SHOCK WAVE SPEED ... 3.62742 KM/S 

INTERFACE PRESSURE 13638E-01 GPA (134.6 ATM) 

TOT TRANS PARTICLE VELOCITY 3.38979 KM/S 

TOT TRANS SHOCK WAVE SPEED 3.95708 KM/S 

TOT TRANS PRESSURE 16365E-01 GPA (161.5 ATM) 

INITIAL SHOCK WAVE PASSAGE . . . 

TOTAL SHOCK WAVE VELOCITY . . 

INTERNAL SHOCK WAVE VELOCITY 
INITAL PRESSURE INCREASE . . . 

INITIAL TEMPERATURE INCREASE 

PEAK SHOCK VALUES . . . 

DIST (CM) SW VEL (KM/S) SI 
.67 3.957 


3.957 KM/S 
3.957 KM/S 
161.507 ATM 
5789.626 DEG-C 


DEN (GM/CC) PR INC (ATM) T INC (DEG-C) 
. 851E-02 161.507 5789.626 
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POINT-OF- INTEREST COORDINATES . . . 

X-COORD = 203.00 CM 
Y-COORD = 61.00 CM 

Z-COORD = .00 CM 

M-CONST = .60 

R-DIST = 211.97 CM 

ATTENUATED SHOCK VALUES . . . 

SH WAVE VEL = .333 KM/S 

SH DENSITY = .123E-02 GM/CC 

PRESS INCR = .011 ATM ( .166 PSI) 

TEMP INCR = .701 DEG-C 

DEBRIS CLOUD PASSAGE . . . 

INT DEB CLD LEADING EDGE VELOCITY 

DEB CLD INDUCED SHOCK WAVE VELOCITY 

INIT PRESS INCR BEH DEB CLD INDUCED SW . . . 

INIT TEMP INCR BEH DEB CLD INDUCED SW 

PR WALL MASS CONTRIB TO DIFF DEB CLD 

DIFFUSE DEBRIS CLOUD MASS 

DIFF DEB CLD AXIAL CHARACT VELOCITY 

DIFFUSE DEB CLD AXIAL KINETIC ENERGY 

DIFF DEB CLD EXPAN CHARACT VELOCITY 

DIFFUSE DEB CLD EXPAN KINETIC ENERGY 

DIFFUSE DEB CLD TOTAL KINETIC ENERGY 

PEAK SHOCK VALUES . . . 

DIST (CM) SW VEL (KM/S) SH DEN (GM/CC) PR INC (ATM) T INC (DEG-C) 
1.01 3.082 . 737E-02 95.422 3360.483 

POINT-OF- INTEREST COORDINATES . . . 

X-COORD = 203.00 CM 
Y-COORD = 61.00 CM 

Z-COORD = .00 CM 

M-CONST = .60 

ETA-DST = 199.43 CM 

ATTENUATED SHOCK VALUES . . . 

SH WAVE VEL = .378 KM/S 

SH DENSITY = .138E-02 GM/CC 

PRESS INCR = .207 ATM ( 3.043 PSI) 

TEMP INCR = 10.941 DEG-C 

CUMULATIVE EFFECT OF SHOCK WAVE AND DEBRIS CLOUD PASSAGE 

TOT PRESSURE INCREASE ... .22 ATM ( 3.208 PSI) 

TOT TEMPERATURE INCREASE ... 11.64 DEG-C 

***** NOW ANALYZING POINT-OF-INTEREST NO. 4/4 ***** 

ENERGY TRANSFORMED INTO SW MOTION 2986E+05 JOULES (50.74%) 

MODIFIED INTERFACE VALUES . . . 

INTERFACE PARTICLE VEL ... 3.08170 KM/S 
TRANS SHOCK WAVE SPEED ... 3.62742 KM/S 

INTERFACE PRESSURE 13638E-01 GPA (134.6 ATM) 

TOT TRANS PARTICLE VELOCITY 3.38979 KM/S 

TOT TRANS SHOCK WAVE SPEED 3.95708 KM/S 

TOT TRANS PRESSURE 16365E-01 GPA (161.5 ATM) 


3.082 KM/S 
3.082 KM/S 
95.422 ATM 
3360.483 DEG-C 

2.817 GMS ( .30* .19MW1) 
7.332 GMS 
2.581 KM/S 

. 2442E+05 JOULES (41.50%) 
.501 KM/S 

. 9183E+03 JOULES ( 1.56%) 

. 2534E+05 JOULES (43.06%) 
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INITIAL SHOCK WAVE PASSAGE . . . 

TOTAL SHOCK WAVE VELOCITY 3.957 KM/S 

INTERNAL SHOCK WAVE VELOCITY 3.957 KM/S 

INITAL PRESSURE INCREASE 161.507 ATM 

INITIAL TEMPERATURE INCREASE 5789.626 DEG-C 

PEAK SHOCK VALUES . . . 

DIST (CM) SW VEL (KM/S) SH DEN (GM/CC) PR INC (ATM) 
• 67 3.957 . 851E-02 161.507 


POINT-OF- INTEREST COORDINATES . . . 
X-COORD = 203.00 CM 
Y-COORD = .00 CM 

Z -COORD = .00 CM 

M-CONST = . 60 

R-DIST = 203.00 CM 


ATTENUATED SHOCK 
SH WAVE VEL = 
SH DENSITY = 
PRESS INCR = 
TEMP INCR 


VALUES . . . 

.333 KM/S 
. 123E-02 GM/CC 
.012 ATM ( 

.744 DEG-C 


.176 PSI) 


DEBRIS CLOUD PASSAGE . . . 

INT DEB CLD LEADING EDGE VELOCITY 3.082 KM/S 

DEB CLD INDUCED SHOCK WAVE VELOCITY 3.082 KM/S 

INIT PRESS INCR BEH DEB CLD INDUCED SW ... 95.422 ATM 

INIT TEMP INCR BEH DEB CLD INDUCED SW 3360.483 DEG-C 

PR WALL MASS CONTRIB TO DIFF DEB CLD 2.817 GMS ( 

DIFFUSE DEBRIS CLOUD MASS 7.332 GMS 

DIFF DEB CLD AXIAL CHARACT VELOCITY 2.581 KM/S 

DIFFUSE DEB CLD AXIAL KINETIC ENERGY 2442E+05 JOULES 

DIFF DEB CLD EXPAN CHARACT VELOCITY .501 KM/S 

DIFFUSE DEB CLD EXPAN KINETIC ENERGY 9183E+03 JOULES 

DIFFUSE DEB CLD TOTAL KINETIC ENERGY 2534E+05 JOULES 


PEAK SHOCK VALUES . . . 

DIST (CM) SW VEL (KM/S) SH DEN (GM/CC) 
1.01 3.082 . 737E-02 


PR INC (ATM) 
95.422 


POINT-OF- INTEREST COORDINATES . . . 
X-COORD = 203.00 CM 
Y-COORD = .00 CM 

Z -COORD = . 00 CM 

M-CONST = .60 

ETA-DST = 129.43 CM 


ATTENUATED SHOCK VALUES . . . 


SH WAVE VEL = 
SH DENSITY = 
PRESS INCR = 
TEMP INCR 


.411 KM/S 
. 150E-02 GM/CC 
.394 ATM ( 
18.842 DEG-C 


5.790 PSI) 


CUMULATIVE EFFECT OF SHOCK WAVE AND DEBRIS CLOUD PASSAGE 
TOT PRESSURE INCREASE ... .41 ATM ( 5.966 PSI) 

TOT TEMPERATURE INCREASE ... 19.59 DEG-C 


INC (DEG-C) 
5789.626 


30* . 19MW1 ) 


(41.50%) 

( 1.56%) 
(43.06%) 


INC (DEG-C) 
3360.483 
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